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ABSTRACT 


This  report  covers  research  in  the  fields  of  radio  physics  and  astronomy 
at  Lincoln  Laboratory  for  the  period  from  1  December  1965  through  31  May 
1966.  In  a  few  cases,  material  from  earlier  work  is  introduced  to  clarify 
the  discussion  of  a  current  program,  e.g.,  " Experimental  Evaluation  of  a 
1000-Wavelength  Antenna."  Twenty-three  papers  presenting  work  by  the 
Radio  Physics,  Communications,  Solid  State  and  Radar  Divisions  of  Lincoln 
Laboratory  and  by  the  M.I.T.  Research  Laboratory  of  Electronics  with  the 
cooperation  of  Lincoln  Laboratory  are  included.  Twenty  of  the  papers  dis¬ 
cuss  work  involving  the  Millstone  Hill  Field  Station  which  is  the  Labora¬ 
tory's  primary  facility  for  experimental  research  in  radio  physics  and  as¬ 
tronomy.  Two  papers  describe  the  28-ft  millimeter  research  facility 
(1000  wavelength)  at  Lincoln  Laboratory  in  Lexington  and  one  is  a  theoret¬ 
ical  paper  on  scattering  from  objects  with  irregular  surfaces. 
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FOREWORD 


This  report  is  the  second  in  a  series  on  research  in  the  field  of  radio  physics 
and  astronomy  being  carried  out  at  Lincoln  Laboratory.  It  covers  activities 
for  a  period  of  approximately  six  months.  Our  objective  for  this  series  is  to 
provide  a  medium  for  the  presentation  of  scientific  results,  technical  and  en¬ 
gineering  developments,  studies,  and  background  reviews;  in  part,  this  re¬ 
port  serves  as  an  observatory  report  for  the  Millstone  Hill  Field  Station.  For 
reasons  of  expedience  in  getting  the  series  started,  the  first  report  contained 
only  descriptions  of  programs  being  carried  out  by  the  Radio  Physics  Division 
at  the  Millstone  Hill  Field  Station.  It  was  intended,  however,  that  these  re¬ 
ports  also  include  allied  work  performed  by  other  Lincoln  Laboratory  divi¬ 
sions  as  well  as  similar  work  by  M.I.T.  personnel  using  Lincoln  Laboratory's 
facilities.  In  that  respect  the  present  report,  which  contains  substantial  con¬ 
tributions  from  other  Lincoln  Laboratory  divisions  and  the  Research  Labora¬ 
tory  of  Electronics,  meets  these  original  objectives.  However,  the  Millstone 
Hill  Field  Station  is  the  Laboratory's  primary  research  facility  in  the  radio 
physics  and  astronomy  fields  and  for  this  reason  descriptions  of  its  programs 
still  account  for  the  major  portion  of  this  report. 
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I.  SUMMARY 


A.  INTRODUCTION 

Following  the  precedent  set  by  our  first  Radio  Physics  and  Astronomy  Report,  the  material 
in  the  present  report  is  organized  under  four  areas  of  research:  Radar  Astronomy,  Radio  As¬ 
tronomy,  Atmospheric  Research  and  Space  Surveillance  Techniques.  A  fifth  section,  Systems 
Development,  has  been  added.  Two  papers  which  do  not  clearly  fall  in  any  of  the  above  classi¬ 
fications  are  included  as  appendices.  One  of  these  papers  gives  a  detailed  description  of  eval¬ 
uation  tests  performed  several  years  ago  on  the  1000-wavelength  (28-ft  diameter)  antenna  at 
Lincoln  Laboratory.  The  other  paper  presents  an  extension  of  perturbation  theory  for  exam¬ 
ination  of  the  scattering  characteristics  of  irregularly  shaped  objects. 

B.  RADAR  ASTRONOMY 

The  23.2-cm  Millstone  radar  and  the  3.8-cm  Haystack  radar  continued  to  be  employed  in 
lunar*  and  planetary  studies.  In  the  lunar  experiments  the  scattering  laws  at  the  two  wave¬ 
lengths  were  determined  and  studies  of  the  polarizing  properties  of  the  lunar  surface  were  made 
at  23.2  cm. 

The  planetary  radar  studies  consisted  mainly  of  observations  on  the  planet  Venus.  Meas¬ 
urements  were  made  of  range,  cross  section  and  the  spectrum  of  the  reflected  signals  which 
yielded  information  regarding  the  orbit  and  the  scattering  laws  at  the  two  wavelengths. 

C.  RADIO  ASTRONOMY 

Radio  astronomy  continued  to  be  one  of  our  most  fruitful  areas  of  research.  Interesting 
observations  were  made  of  anomalous  spectral  line  emissions  from  OH  at  1665  MHz  near  W3. 
Absolute  flux  measurements  of  Cassiopeia  A  and  Taurus  A  were  made  at  3.64  and  1.94  cm.  The 
results  of  these  flux  measurements  enable  the  use  of  these  sources  for  measuring  the  gains  of 
large  antennas.  Observations  at  these  two  wavelengths  were  also  made  of  the  secular  variations 
in  the  radio  flux  from  3C273.  A  digital  mapping  technique  developed  by  personnel  of  Lincoln 
Laboratory  and  M.I.T.  Research  Laboratory  of  Electronics  was  used  to  map  the  center  of  our 
galaxy. 

The  Haystack  120-foot  antenna  is  the  basic  instrument  for  most  of  our  radio  astronomy 
programs  at  wavelengths  from  hydrogen  line  to  about  1.5  cm.  However,  the  84-foot  antenna  at 
Millstone  has  been  used  with  the  Haystack  antenna  to  form  a  2000-foot  base-line  radiometric 
interferometer.  First  use  of  the  interferometer  will  be  to  determine  the  positions  and  angular 
sizes  of  OH  emission  regions  in  the  1665-MHz  band. 

The  28 -foot  millimeter-wave  antenna  at  Lincoln  Laboratory  in  Lexington  was  used  in  radio- 
metric  studies  of  Venus,  Jupiter,  the  moon  and  the  terrestrial  atmosphere  at  wavelengths  near 
one  centimeter. 


*  The  lunar  studies  program  is  supported  by  NASA. 
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D.  ATMOSPHERIC  RESEARCH 

Studies  of  the  ionosphere  by  the  incoherent  backscatter  technique  continue  to  be  an  important 
area  of  research.  However,  because  of  the  lag  between  data  gathering  operations  and  the  data 
reduction  and  analysis  operations,  there  is  little  to  report  in  this  area  that  was  not  discussed  in 
Radio  Physics  and  Astronomy  Report  (1965:1). 

No  new  data  were  taken  in  the  auroral  studies  program  during  the  period  covered  by  this 
report.  All  the  effort  in  this  area  was  devoted  to  the  reduction  and  analysis  of  auroral  data  taken 
earlier. 

Because  of  the  increased  use  of  X-band  wavelengths  in  satellite  communication  systems, 
there  is  great  interest  in  obtaining  a  better  understanding  of  the  effects  of  precipitation  on  radio 
propagation  at  these  wavelengths.  The  Haystack  X-band  radiometer  has  been  used  by  the  Lincoln 
Laboratory  Communications  Division  in  conjunction  with  the  M.I.T.  weather  radar  research  facil¬ 
ity  to  examine  the  relation  between  radar  reflectivity  and  thermal  effects  for  a  number  of  showers 
observed  in  August  1965. 

E.  SPACE  SURVEILLANCE  TECHNIQUES 

The  space  surveillance  program  for  the  Millstone  L-band  radar  was  discussed  in  detail  in 
the  last  radio  physics  and  astronomy  report  (1965:1).  The  most  important  part  of  this  program 
is  the  development  of  improved  techniques  for  the  detection  and  tracking  of  distant  satellites. 

A  versatile  system  of  radar  data  handling  and  orbit  computation  programs  is  being  developed, 
partially  through  contractual  assistance,  for  the  Millstone  computer.  Two  new  approaches 
have  also  been  studied  relative  to  extending  the  range  of  the  radar  by  computer  integration  tech¬ 
niques  while  retaining  the  ability  to  autotrack. 

A  new  tracking  experiment  using  a  ruby  laser  radar  is  being  conducted  at  the  Millstone  Hill 
Field  Station  by  the  Optics  and  Infrared  Group  of  Lincoln  Laboratory's  Solid  State  Division.  This 
instrument  uses  a  precision  azimuth-elevation  mount  installed  on  a  concrete  silo  adjacent  to  the 
Millstone  radar.  The  mount  is  slaved  to  the  Millstone  L-band  radar,  which  is  used  to  bring  the 
laser  radar  beam  on  target. 

Another  interesting  development  at  the  Millstone  Hill  Field  Station  is  the  experimental  use 
of  the  Haystack  120-ft  dish  for  closed-loop  passive  tracking  of  cooperative  X-band  satellites. 

The  precision  angle  tracking  provided  by  this  facility  can  be  used  to  upgrade  the  orbital  elements 
of  satellites  in  very  high  orbits,  e.g.,  certain  communications  satellites,  which  are  beyond  the 
range  of  conventional  "skin-tracking"  satellite  radars. 

F.  SYSTEMS  DEVELOPMENT 

A  great  effort  is  being  made  to  complete  a  500-kw  average  power  X-band  radar  for  Hay¬ 
stack.  This  radar  will  be  used  initially  in  the  "fourth  test  of  general  relativity"  which  was  pro¬ 
posed  by  1. 1.  Shapiro  and  described  in  the  last  radio  physics  and  astronomy  report. 

A  millimeter-wave  research  facility  is  being  developed  around  the  28-ft  precision  surface 
antenna  at  Lincoln  Laboratory.  This  facility,  when  developed  according  to  present  plans,  will 
include  a  CW  1000-W,  35-GHz  radar. 
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II.  RADAR  ASTRONOMY 


A.  INTRODUCTION 

Observations  of  the  moon  were  made  as  part  of  a  program  sponsored  by  NASA.  An  extensive 

series  of  observations  at  Millstone  forms  the  basis  of  a  paper  (in  press  ).  This  work  consisted 

primarily  of  determining  the  angular  scattering  law  at  23 -cm  wavelength  and  comparing  this 

with  determinations  at  other  wavelengths  (e.g.,  the  68 -cm  results  reported  by  Evans  and 
2 

Pettengill  ).  Preliminary  results  of  new  experiments  in  which  the  depolarizing  properties  of 
the  lunar  surface  were  studied  in  greater  detail  than  hitherto  are  reported  here.  Measurements 
with  the  Haystack  radar  at  3.8  cm  have  also  recently  begun.  A  considerable  amount  of  effort  has 
been  devoted  to  the  programming  required  for  the  lunar  mapping  that  we  plan  to  carry  out  at  Hay¬ 
stack  Hill.  The  scattering  law  at  3.8  cm  has  also  been  redetermined  using  Haystack  and  is  found 
to  be  in  good  agreement  with  the  earlier  3.6-cm  determination  reported  by  Evans  and  Pettengill. 

The  Millstone  radar  was  employed  to  observe  Venus  during  the  past  year,  primarily  in  or¬ 
der  to  secure  precise  range  information  on  a  portion  of  the  orbit  not  previously  observed  by  ra¬ 
dar.  Range  determinations  were  made  for  a  period  of  six  months  centered  on  close  approach. 
During  the  central  four  months,  the  flight  time  was  measured  to  an  accuracy  of  approximately 
±10|isec  about  twice  a  week.  These  measurements  also  form  the  subject  of  a  separate  report. 

Perhaps  the  most  exciting  event  during  the  past  year  for  the  radar  astronomers  at  the 
Field  Station  was  the  availability,  for  the  first  time,  of  a  100-kw,  X-band  CW  radar  system 
at  Haystack  Hill.  A  few  observations  of  Mercury  have  been  made  with  this  radar,  but  most  of 
the  effort  so  far  has  been  devoted  to  studying  Venus.  A  particularly  intensive  series  of  obser¬ 
vations  was  made  during  the  January  1966  inferior  conjunction.  Since  these  measurements  will 
be  reported  separately  in  a  journal  article,  we  only  summarize  the  main  conclusions  here. 

B.  LUNAR  STUDIES 

1.  Polarization  of  Radar  Returns  from  Moon 

a.  Introduction 

Systematic  studies  of  the  state  of  polarization  of  radar  echoes  from  the  lunar  surface  for 
different  types  of  polarization  of  the  transmitted  wave  provide  information  which,  when  properly 
interpreted,  relates  directly  to  the  structure  of  the  lunar  surface.  Depolarization  experiments 
have  been  carried  out  at  A.  =  23  cm  with  the  Millstone  radar  and  have  been  started  at  A.  =  3.8  cm 
using  the  Haystack  radar.  Section  II-B-l-b  briefly  describes  the  equipment  which  provides  the 
necessary  flexibility  for  such  polarization  experiments.  Sections  II-B-l-c  and  d  describe  the 
experiments  and  give  the  results.  Section  II-B-l-e  discusses  their  physical  significance. 

b.  Polarization  Equipment 

Two  circularly  polarized  waves  of  opposite  sense  can  be  excited  in  the  Millstone  antenna 
by  feeding  a  signal  into  either  one  or  the  other  of  two  antenna  input  terminals.  By  providing  a 
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Fig.  11-1.  Circuit  to  control  transmitter  polarization. 


Fig.  11-2.  Polarized  and  depolarized  circular 
components  vs  delay. 


Fig.  11-3.  Percentage  circular  polarization  vs  delay. 
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network  which  can  arbitrarily  divide  the  power  between  the  two  antenna  ports  and  alter  their 
relative  phase,  it  is  possible  to  excite  a  wave  of  arbitrary  polarization.  As  shown  in  Fig.II-1, 
the  signal  from  the  transmitter  is  applied  to  the  90°  3-db  hybrid  HI.  In  one  of  the  two  signal 
paths  of  the  output  of  HI,  a  variable  phase  changer  (pi  is  inserted.  The  two  signals  are  then 
applied  to  another  hybrid  circuit  H2.  The  relative  power  at  the  output  terminals  A  and  B  can 
thus  be  varied  by  changing  the  phase  delay  in  (pi  and  the  relative  phase  of  the  two  signals  applied 
to  the  antenna  terminals  by  changing  the  phase  delay  in  (p 2.  This  combination  network  allows  an 
arbitrary  polarization  to  be  transmitted. 

The  combination  of  the  two  orthogonally  polarized  received  signals  at  the  two  input  terminals 
of  HI  in  a  second  network  similar  to  that  shown  in  Fig.II-1  allows  for  the  reception  of  a  speci¬ 
fied  pair  of  orthogonally  polarized  components  for  any  arbitrary  transmitted  polarization.  The 
synthesis  of  any  arbitrary  polarizations  for  the  received  components  could,  in  principle,  be 
carried  out  in  the  computer  since  quadrature  phase  detection  is  employed  and  detector  outputs 
are  digitized  for  both  polarizations.  This  procedure  was  not  followed,  however,  because  of 
inadequate  phase  and  gain  stability  in  the  receiver  chains. 

The  Haystack  system  operating  at  3.8  cm  is  at  present  capable  only  of  transmitting  circularly 
polarized  waves.  A  waveguide  network  similar  to  that  shown  in  Fig.II-1  is,  however,  available 
for  the  reception  of  an  arbitrary  pair  of  orthogonally  polarized  waves. 


c.  Experiments  and  Results  at  23  cm  at  Millstone 

The  moon  was  illuminated  with  circularly  polarized  waves  and  the  two  opposite  sense  cir¬ 
cularly  polarized  received  waves  were  observed.  The  ellipticity  of  the  transmitted  wave  was 
adjusted  by  means  of  a  near-field  test  antenna  to  give  a  ratio  between  major  and  minor  axes  of 
less  than  0.4  db.  The  ellipticity  of  the  two  receiver  channels  was  typically  0.6  - 0.8  db  and  a 
correction  was  applied  in  the  data  analysis  program  to  account  for  this.  Figure  II -2  shows  a 
plot  of  both  polarized  and  depolarized  components  as  a  function  of  delay  beyond  the  subradar 
point.  Range  resolution  was  lOOpsec  to  a  delay  of  4  msec  and  500psec  for  the  remainder  of  the 
curve.  Figure  II -3  shows  the  percentage  polarization  measured  as  a  function  of  delay  and  com¬ 
pares  the  results  with  previous  observations  at  68  cm  (Ref.  2).  The  polarization  is  defined  by 


P  = 


-  P 

TV 


L 

L 


where  P„,  PT  are  the  powers  in  the  right  and  left  circularly  polarized  receiver  channels, 

R  E 

respectively. 

A  second  experiment  consisted  in  illuminating  the  moon  with  linearly  polarized  radiation 
and  measuring  the  ratio  of  the  polarized  to  the  depolarized  components  of  the  received  signal. 
Because  of  uncertainty  in  the  amount  of  Faraday  rotation,  the  receiver  polarization  was  rotated 
between  runs.  The  relative  power  output  of  the  two  receiver  channels  will  then  vary  sinusoi¬ 
dally  with  the  angle  2 ft,  where  ft  is  the  feed  rotation  angle.  A  least-mean-square  sine  wave 
was  fitted  to  the  data  obtained  in  different  runs  and  the  depth  of  modulation  was  used  to  deter¬ 
mine  the  ratio  of  the  depolarized  to  the  polarized  components.  Figure  II-4  shows  the  two  com¬ 
ponents  plotted  as  a  function  of  delay.  The  spike  in  the  depolarized  component  near  zero  delay 
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Fig.  11-4.  Polarized  and  depolarized  linear 
components  vs  delay. 


Fig.  11-5.  Degree  of  linear  polarization  vs  delay. 


Fig.  11-6.  Angular  scattering  law  for  depolarized  circular  component. 
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is  at  present  believed  to  be  caused  by  insufficient  isolation  between  the  two  channels.  The 
degree  of  linear  polarization  plotted  as  a  function  of  delay  is  shown  in  Fig.  II-5.  The  definition 
of  polarization  is  similar  to  that  in  the  case  of  circular  polarization,  except  that  PDPT  have  to 
be  replaced  by  the  power  in  the  two  linearly  polarized  orthogonal  channels.  Note  that  the  linear 
depolarization  is  much  smaller  than  the  circular  depolarization. 

At  least  one  peculiarity  of  the  data  just  described  is  worth  mentioning.  When  the  depolarized 
circular  component  in  Fig.II-2  -  after  appropriate  correction  for  antenna  polar  diagram  effects  - 
is  plotted  as  a  function  of  cos  0,  where  0  is  the  angle  of  incidence,  the  power -range  relation¬ 
ship  very  closely  follows  a  cos  0  law,  as  shown  in  Fig.II-6.  This  corresponds  to  a  uniformly 
bright  disk.  When  the  depolarized  linear  component  in  Fig.II-4  after  similar  corrections  is 
plotted  in  the  same  way,  it  turns  out  that  the  scattering  law  goes  as  cos3//2  0,  as  shown  in 
Fig.  II  -7. 


Fig.  11-7.  Angular  scattering  law  for  depolarized  linear  components. 

In  the  second  of  the  two  experiments  just  described,  the  transmitted  linear  polarization 
will  lie  within  the  local  plane  of  incidence  at  two  points  on  a  constant  delay  ring  and  will  be  nor¬ 
mal  to  the  local  plane  of  incidence  at  two  other  points  90°  away.  At  intermediate  points  on  the 
range  ring  considered,  the  direction  of  polarization  of  the  illumination  with  respect  to  the  local 
plane  of  incidence  will  be  intermediate  between  0°  and  90°. 

In  order  to  examine  the  effect  of  the  orientation  of  the  transmitted  linear  polarization  with 
respect  to  the  local  plane  of  incidence,  the  linear  polarization  of  the  transmitter  was  programmed 
to  be  parallel  to  the  instantaneous  lunar  libration  axis.  The  resulting  measured  ratio  of  depolar- 
ized-to-polarized  components  is  shown  as  a  function  of  delay  in  Fig.  II- 8,  where  it  is  also  compared 


7 


Section  II 


INCOHERENT  DATA 

\ 

|3 -30-10248) 

.  • 

.SLX> _ °_L  < 

o°#**oo° 

0  0 

°l 

C 

O 

l 

/  ^  o  o  X 

P'0  0 

• 

x.  • 

COHERENT 

DATA  X 

•  •  N. 

MOON  ECHO,  X  =  23  cm 

TRANSMIT  LINEAR  II  LIBRATION  AXIS 

RECEIVE  ORTHOGONAL 

LINEARS 

5  APRIL  1966 

_ 1 _ 

1 

6  8 
DELAY  (msec) 


Fig.  11-8.  Ratio  of  depolarized-to-polarized  components  when 
surface  is  illuminated  with  a  wave  polarized  linearly  in  local 
plane  of  incidence. 


Fig.  11-9.  Square  root  of  ratio  of  backscattering  coefficients 
plotted  against  angle  of  incidence  for  X  =  23  and  3.7cm. 
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with  the  ratio  for  a  complete  range  ring  (obtained  from  incoherent  processing  of  the  data).  As 
explained  in  the  next  section,  the  data  shown  in  Fig.II-8  may  be  in  qualitative  agreement  only 
with  previous  work  and  the  measurements  will  be  repeated  to  reduce  noise  and,  if  possible,  to 
eliminate  any  systematic  errors  in  the  experiment. 

d.  Results  Obtained  at  Haystack  at  3.8  cm 

The  measurement  of  the  difference  of  backscattering  coefficients  for  waves  polarized  in  and 
across  the  local  plane  of  incidence  previously  carried  out  at  23  cm  (Ref.  3)  was  repeated  at  3.8  cm, 
using  a  slightly  different  technique.  Whereas  the  original  23-cm  experiment  used  frequency 
resolution  to  separate  small  areas  with  a  well-defined  plane  of  incidence,  the  3.8 -cm  experi¬ 
ment  depended  on  the  narrow  polar  diagram  of  the  antenna  beam  for  this  purpose.  As  before, 
the  antenna  was  made  to  transmit  a  circularly  polarized  wave.  Two  orthogonal  linearly  polar¬ 
ized  channels  were  synthesized  in  a  manner  similar  to  that  shown  in  Fig.  II -1.  One  channel 
was  aligned  for  maximum  sensitivity  to  correspond  with  E-field  parallel  to  the  disk-radial  di¬ 
rection  connecting  the  center  of  the  disk  and  the  point  of  observation.  The  variation  of  the  power 
ratio  in  the  two  received  components  with  angle  of  incidence  was  derived  by  moving  the  antenna 
beam  radially  from  the  center  of  the  disk  to  the  limb  with  a  fixed  polarization  of  the  two  re¬ 
ceiver  channels.  Figure  II -9  shows  the  square  root  of  the  ratio  of  the  two  received  components 
as  a  function  of  angle  of  incidence.  Also  shown  is  the  corresponding  curve  obtained  previously 

3 

at  23-cm  wavelength.  As  can  be  seen,  the  polarization  effect  is  significantly  smaller  at  X-band 
than  at  L-band. 

e.  Discussion 

The  depolarization  of  the  circularly  polarized  wave  may  be  explained  in  at  least  two  ways: 

(1)  a  difference  in  the  backscattering  coefficients  of  the  two  principal  linearly  polarized  com¬ 
ponents  but  no  coupling  from  one  into  the  other,  or  (2)  a  "spillover"  from  one  of  the  principal 
linear  components  into  the  other.  Table  II -1  shows  a  few  numerical  values  for  hypothesis  (1). 


TABLE 

ll-l 

PREDICTIONS  FOR  HYPOTHESIS  (1) 

If  Ratio  of 

Circular  Returns 
for  Circular 
Illumination  Is  (in  db): 

Ratio  of 
Backscattering 
Coefficients  for 
Linear  Components 
Should  Be  (in  db): 

2.0 

18.8 

4.0 

13.0 

6.0 

9.6 

8.0 

7.3 

10.0 

5.6 
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Since  the  ratio  of  the  two  orthogonal  backscattering  coefficients  for  linear  polarization  has 
previously  been  shown  to  be,  at  most,  about  3.0  db  (Ref.  3),  we  see  that  hypothesis  (1)  is  com¬ 
pletely  inadequate  to  account  for  the  observed  depolarization  effect,  which,  for  circular  trans¬ 
mitted  and  circular  received  waves,  corresponds  to  a  ratio  of  the  two  components  of  only  between 
2  and  3  db  (compare  Fig.  II- 2) . 

Hypothesis  (2)  appears  to  be  much  closer  to  the  actual  situation  than  (1)  in  accounting  for 

the  circular-circular  depolarization  effect  and  it  may  also  largely  explain  the  observed  linear- 

linear  depolarization  effect.  One  might  imagine  that  the  scattering  at  oblique  incidence  occurs 

from  a  discrete  structure  of  irregular  shape,  forcing  the  induced  currents  to  flow  in  directions 

which  differ  from  that  of  the  impressed  electric  field.  As  a  mathematical  model,  one  might 

imagine  the  scattering  to  occur  from  a  collection  of  randomly  oriented  dipoles  and  a  collection 

\ 

of  nondepolarizing  ideal  reflectors.  On  the  basis  of  such  a  naive  model,  one  can  work  out  the 
ratio  of  total  scattered  power  in  the  dipole  component  and  the  ideally  reflected  component  for 
various  ratios  of  depolarized-to-polarized  components  for  circular  and  for  linear  illumination. 

A  set  of  predictions  based  on  this  dipole-reflector  model  is  given  in  Table  II  —  2 .  Toward  grazing 
angles  of  incidence  on  the  surface,  it  appears  that  a  percentage  of  power  from  the  randomly 
oriented  dipoles  of  about  70  percent  will  explain  both  linear  and  circular  polarization.  This 
simple  model  does  not  explain  completely  the  variation  of  depolarization  with  angle  of  incidence 
because  the  linear/linear  depolarization  remains  virtually  constant  for  angles  of  incidence 
45  °  <  (p  <  90  °,  whereas  the  circular/circular  depolarization  increases  significantly  in  this 
range  of  angles. 

The  results  shown  in  Fig.  II -B  can  be  understood  at  least  qualitatively  on  the  basis  of  pre¬ 
vious  experiments  as  follows.  Denote  the  backscattering  coefficients  for  in-  and  across-plane 
of  incidence  polarization  by  avv  and  cthh.  By  reciprocity  the  "spillover”  coefficients  and 


TABLE  11-2 

PREDICTION  BASED  ON  DIPOLE-REFLECTOR  MODEL 

Function  of 

Ratio  of  Depolarized/ 

Ratio  of  Depolarized/ 

Power  Scattered 

Polarized,  Linear 

Polarized,  Circular 

From  Dipoles  (%) 

(linear  illumination)  (db) 

(circular  illumination)  (db) 

10 

15.9 

12.8 

20 

12.8 

9.5 

30 

10.9 

7.5 

40 

9.5 

6.0 

50 

8.4 

4.7 

60 

7.5 

3.6 

70 

6.6 

2.7 

80 

6.0 

1.8 

90 

5.2 

0.  85 

100 

4.  7 

0 
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o'  r  must  be  equal.  The  decrease  in  the  ratio  of  depolarized-to-polarized  components  with 
H  V 

respect  to  the  average  obtained  from  a  whole  range  ring  can  be  interpreted  in  terms  of  the 
systematic  variation  in  °VV  angle  of  incidence  as  previously  reported." 

The  less  pronounced  variation  of  the  ratio  of  the  two  principal  backscattering  coefficients 
with  angle  of  incidence  at  3.8  cm  than  at  23  cm  as  shown  in  Fig.  II-9  may  be  accounted  for  in 
several  different  ways.  It  may  be  that  the  scattering  from  within  the  tenuous  top  layer  occurs 
from  irregular  structure  actually  suspended  in  the  layer.  If  the  tenuous  layer  has  a  smaller 
dielectric  constant  at  the  surface  than  deeper  in  the  layer,  then  the  fact  that  the  ratio  tends 
toward  unity  with  decreasing  wavelength  may  be  taken  as  evidence  of  the  gradual  change  in  the 
material  properties  with  depth.  It  may  also  be  that  a  sufficient  number  of  small  rocks  or  pebbles 
lying  on  top  of  the  tenuous  layer  may  account  for  a  significant  amount  of  scattering  at  3.8  cm  but 
not  at  23  cm,  and  this  may  explain  the  observations. 

2.  Lunar  Scattering  Law  at  3.8 -cm  Wavelength 

The  observations  of  the  lunar  scattering  behavior  at  23 -cm  wavelength  as  described  so  far 
only  represent  a  small  part  of  the  work  required  by  our  contract  with  NASA.  Additional  meas¬ 
urements  are  in  progress  at  Haystack  to  examine  the  behavior  at  3. 6 -cm  wavelength. 

The  moon's  scattering  is  known  to  change  markedly  as  a  function  of  wavelength.  This  arises, 
it  is  thought,  partly  because  the  surface  is  covered  with  an  increasing  amount  of  small  scale 
structure  and  hence  appears  "rougher"  as  the  wavelength  is  reduced,  and  partly  because  the 
surface  appears  inhomogeneous  and  hence  different  wavelengths  penetrate  to  different  depths. 

4 

Previously  X-band  observations  were  made  of  the  moon  in  1962  by  Evans  using  the  Lincoln 
Laboratory  —  Camp  Parks  radar.  Better  accuracy  is  expected  from  the  Haystack  system  because 
of  the  better  pointing  accuracy  and  superior  signal-to-noise  ratio  provided  by  the  Haystack  radar, 
coupled  with  the  better  range  tracking  and  integration  techniques  available  with  the  Millstone 
radar  sequencer  and  SDS  9300  computer  combination. 

Observations  were  made  with  Haystack  each  morning  from  Wednesday,  16  March  through 
Saturday,  19  March  for  about  7  hours.  At  Millstone  the  antenna  beam  illuminates  essentially 
the  whole  moon  so  that  the  complete  law  for  the  variation  of  echo  power  with  delay  can  be  es¬ 
tablished  with  the  beam  directed  at  the  moon's  center.  In  practice  it  is  usually  necessary  to 
employ  a  variety  of  pulse  lengths  so  that  good  delay  resolution  can  be  obtained  for  the  region 
around  the  subradar  point,  and  yet  good  signal-to-noise  ratio  may  be  had  for  the  limbs.  At 
Haystack  the  pulse  length  must  be  varied  and,  in  addition,  the  beam  must  be  scanned  off  center. 

In  the  measurements  made  in  March,  a  pulse  of  10  psec  with  a  10-psec  associated  sample  spacing 
was  employed  to  explore  the  central  region.  Measurements  were  made  for  periods  of  5  minutes 
for  beam  positions  of  0  to  5  arc  minutes  with  respect  to  the  center  in  1 -arc-minute  steps.  The 
beam  was  scanned  away  from  the  center  in  four  mutually  perpendicular  directions.  The  region 
6  to  12  arc  minutes  from  the  center  was  later  explored  in  1 -arc-minute  steps,  using  40-psec 
pulses  and  an  associated  20-psec  sample  interval.  The  beam  was  scanned  away  from  the  center 
in  8  directions,  intersecting  at  45°  angles.  Finally,  the  limb  region  12  to  15  arc  minutes  from 
the  center  was  explored  in  the  same  8  directions,  using  a  100-psec  sample  interval. 
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TABLE  11-3 

RELATIVE  ECHO  POWER  VS  DELAY  FOR  THE  MOON 
OBSERVED  AT  LINCOLN  LABORATORY  M.l.T. 


Delay  (t) 
(nsec) 

9° 

X  =  3.  6  cm* 

23  cmt 

68  cmt 
(db) 

Delay  (t) 
(msec) 

9° 

X  =  3.  6  cm 

23  cm 

68  cm 
(db) 

10 

2.38 

0 

0 

0 

2.0 

34.  16 

-12.5 

-18.65 

-21.8 

20 

3.37 

-0.6 

-0.85 

-0.6 

2.25 

36.30 

-12.9 

-19.  1 

-22.3 

30 

4.  11 

-1.15 

-1.40 

-1.5 

2.50 

38.33 

-13.  35 

-19.  55 

-22.7 

40 

4.  77 

-1.55 

-1.9 

-2.2 

2.75 

40.28 

-13.75 

-19.  85 

-23.  1 

50 

5.31 

-1.9 

-2.35 

-2.8 

3.0 

42.  16 

-14.05 

-20.2 

-23.5 

60 

5.83 

-2.2 

-2 .75 

-3.3 

3.25 

43.97 

-14.4 

-20.5 

-23.8 

70 

6.30 

-2.5 

-3.  15 

-3.8 

3.50 

45.  72 

-14.7 

-20.  85 

-24.  1 

80 

6.73 

-2.75 

-3.55 

-4.3 

3.75 

47.  42 

-15.0 

-21.  15 

-24.3 

90 

7.  13 

-3.0 

-3.  95 

-4.8 

4.0 

49.08 

-15.3 

-21.4 

-24.5 

100 

7.  53 

-3.2 

-4.3 

-5.2 

4.25 

50.69 

-15.55 

-21. 7 

-24.  7 

125 

8.42 

-3.75 

-5. 05 

-6.2 

4.50 

52.27 

-15.8 

-21.95 

-24.9 

150 

9.22 

-4.2 

-5.8 

-7.0 

4.75 

53.  82 

-16.  1 

-22.35 

-25.  1 

175 

9.96 

-4.  55 

-6.45 

-7.7 

5.0 

55.  33 

-16.3 

-22.5 

-25.  35 

200 

10.65 

-4.9 

-7.0 

-8.4 

5.  25 

56.  82 

-16.6 

-22.  75 

-25.6 

225 

11.30 

-5.2 

-7.5 

-9.0 

5.50 

58.29 

-16.9 

-23.05 

-25.9 

250 

11.92 

-5.5 

-8.0 

-9.  7 

5.75 

59.  72 

-17.  15 

-23.3 

-26.  15 

275 

12.50 

-5.  85 

-8.  45 

-10.  15 

6.0 

60.95 

-17.5 

-23.6 

-26.  45 

300 

13.06 

-6.  1 

-8.  85 

-10.7 

6.25 

62.55 

-17.8 

-23.  85 

-26.8 

325 

13.60 

-6.4 

-9.2 

-11.1 

6.50 

63.93 

-18.  1 

-24.2 

-27.  15 

350 

14.  11 

-6.65 

-9.55 

-11.6 

6.75 

65.  36 

-18.5 

-24.55 

-27.5 

375 

14.60 

-6.9 

-9.9 

-11.9 

7.0 

66.65 

-18.9 

-24.  95 

-27.  85 

400 

15.09 

-7.  15 

-10.2 

-12.35 

7.  25 

67.  98 

-19.3 

-25.  35 

-28.  25 

425 

15.56 

-7.4 

-10.5 

-12.7 

7.50 

69.31 

-19.7 

-25.8 

-28.65 

450 

16.01 

-7.6 

-10.8 

-13.0 

7.75 

70.63 

-20.2 

-26.2 

-29.0 

475 

16.45 

-7.8 

-11.2 

-13.3 

8.0 

71.94 

-20.6 

-26.65 

-29.45 

500 

16.88 

-7.  95 

-11.35 

-13.6 

8.25 

73.23 

-21.05 

-27.  1 

-29.  95 

600 

18.51 

-8.5 

-12.3 

-14.6 

8.50 

74.  53 

-21.55 

-27.6 

-30.45 

700 

20.01 

-9.0 

-13.2 

-15.4 

8.75 

75.  79 

-22.05 

-28.  15 

-30.95 

800 

21.40 

-9.4 

-14.0 

-16.2 

9.0 

77.  06 

-22.55 

-28.  7 

-31.5 

900 

22.  72 

-9.8 

-14.6 

-16.9 

9.25 

78.33 

-23.  15 

-29.3 

-32.  1 

1000 

23.96 

-10.  1 

-15.  35 

-17.6 

9.  50 

79.59 

-23.  7 

-29.  95 

-32.  75 

1100 

25.  16 

-10.25 

-15.  95 

-18.  25 

9.75 

80.84 

-24.5 

-30.6 

-33.  35 

1200 

26.29 

-10.  75 

-16.4 

-18.9 

10.0 

82.09 

-25.2 

-31.35 

-34.05 

1300 

27.  39 

-11.0 

-16.8 

-19.5 

10.25 

83.  34 

-26.  1 

-32.2 

-34.  9 

1400 

28.44 

- 1 1 . 25 

-17.  15 

-20.0 

10.50 

84.  58 

-27.  1 

-33.25 

-35.9 

1500 

29.42 

-11.5 

- 1  7.  45 

-20.4 

10.  75 

85.82 

-28.5 

-34.35 

-36.9 

11.0 

87.  06 

-30.0 

-35.  85 

-38.  35 

11.25 

88.  29 

-31.2 

-37.7 

-40.  1 

11.50 

89.  53 

-32.0 

-40.35 

-42.  7 

*  Pulse  =  10  ysec;  resolution  at  the  receiver  =  10  ^isec. 

t  Pulse  =  10  ysec;  resolution  at  the  receiver  =  10  jjsec. 

t  Pulse  =  12  |jsec;  resolution  at  the  receiver  =  10  psec. 
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In  all  these  observations  the  antenna  offsets  from  the  moon's  center  were  inserted  manually. 
The  timing  of  the  Haystack  radar  was  locked  to  that  at  Millstone  so  that  the  Millstone  radar 
sequencer  (Craig/Fleck  machine)  could  be  used  to  sample  the  echoes  as  in  earlier  Venus  runs. 

One  change  for  the  moon  was  that  the  IF  signal  transmitted  to  Millstone  via  the  intersite  coupling 
was  at  2  MHz  instead  of  200  kHz,  so  that  adequate  bandwidth  (100  kHz)  could  be  obtained  for  the 
10-psec  pulse  measurements.  The  runs  were  processed  at  Millstone  using  the  SDS  9300  computer 
to  obtain  the  integrated  echo  power  vs  delay.  The  speed  of  this  program  is  such  that  the  data 
integration  is  completed  shortly  after  the  end  of  the  operation. 

Only  one  difficulty  was  encountered  during  these  operations.  This  was  the  failure  of  the 
power  amplifier  pulse  modulator  to  operate  continuously  above  approximately  30 -kv  applied 
voltage.  It  was  necessary  to  bring  up  the  voltage  slowly  during  the  first  hour  in  order  to  reach 
75-kw  peak  power  output.  Even  at  this  level,  however,  the  modulator  safety  circuits  would  trip 
intermittently. 

The  data  obtained  in  the  manner  described  above  have  been  analyzed  to  obtain  the  plot  of 
echo  power  vs  delay  shown  in  Fig.  II- 10.  In  order  to  obtain  this  plot,  it  was  necessary  to  average 
all  the  echo  power  vs  delay  plots  obtained  for  a  given  antenna  offset  at  different  angular  locations 
on  the  lunar  disk  and  then  to  correct  this  mean  curve  for  the  weighting  imposed  by  the  antenna 
beam.  In  this  way  portions  of  the  total  curve  were  built  up.  A  check  on  the  accuracy  of  the 
procedure  was  afforded  by  a  comparison  of  the  curves  in  the  regions  where  they  overlap.  A 
change  in  the  assumed  beam  shape  was  required  before  such  consistency  was  achieved.  The 
curve  shown  in  Fig.  11-10  is  in  reasonable  agreement  with  the  earlier  results,  but  provides 
better  resolution  for  the  delay  region  of  10  to  200  psec. 


Fig.  11-10.  Lunar  echo  power  vs  delay 
at  X  =  3.  8  cm. 


Table  II— 3  lists  the  scattering  laws  obtained  for  the  moon  at  X  =  68  cm,  X  =  23  cm  and 
X  =  3.8  cm,  using  Lincoln  Laboratory  radars.  Since  the  resolution  in  delay  was  the  same  in  all 
three  sets  of  measurements,  these  data  should  accurately  reflect  the  change  in  scattering  be¬ 
havior  with  wavelength.  j  v  Evans 

T.  Hagfors 
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C.  PLANETARY  STUDIES 

1.  Radar  Observations  of  Venus  at  Millstone 

Venus  was  observed  at  Millstone  (X  =  23  cm)  from  22  October  1965  to  7  April  1966.  The 

range  was  determined  to  an  accuracy  of  ±10  psec  from  27  December  to  21  March  at  approximately 

3-  to  4-day  intervals.  The  measurements  terminated  somewhat  earlier  than  anticipated  owing 

to  poor  signal-to-noise  ratio,  which  was  later  attributed  to  damage  caused  by  the  transmitter 

to  the  receiver  coaxial  line.  These  range  determinations  together  with  those  obtained  in  1964 

at  Millstone,  and  in  1964  and  1965  at  the  Arecibo  Ionospheric  Observatory  have  been  used  to 

5 

refine  the  elements  of  the  orbits  of  earth  and  Venus  and  redetermine  the  earth/moon  mass  ratio. 
A  preliminary  reduction  of  the  data  has  yielded  a  value  of  81.305  ±  0.005  for  this  latter  ratio 
(Smith,  et  al.^)  which  is  in  good  agreement  with  the  recently  adopted  value. 

The  scattering  law  for  Venus  was  determined  in  1964  (Evans,  et  al.,7),  using  a  combination 

g 

of  short  pulse  measurements.  Russian  workers  have  reported  a  change  in  the  scattering  law 
from  one  year  to  the  next  at  a  40-cm  wavelength.  In  our  view  this  seems  most  unlikely  because 
Venus  rotates  in  such  a  manner  that  approximately  the  same  face  is  presented  toward  the  earth 
at  successive  inferior  conjunctions.  Thus  the  change  in  the  scattering  behavior  reported  by  the 
Russian  workers  is,  we  feel,  more  readily  explained  as  an  instrumental  effect  arising  from  the 
much  improved  delay  resolution  that  they  achieved  in  1964.  Nevertheless,  during  the  1966 
conjunction  the  scattering  law  at  23-cm  wavelength  was  redetermined  by  using  methods  essen¬ 
tially  the  same  as  those  employed  in  1964.  Figure  II- 1 1  compares  the  new  and  old  results.  It 
can  be  seen  that  to  the  accuracy  of  these  measurements  there  is  no  difference  between  the  two 
sets  of  measurements. 


Fig.  11-11.  Echo  power  vs  delay  law  observed  in 
1966  is  here  compared  with  that  derived  by  Evans, 
et  o\7  during  1964.  Agreement  between  two 
curves  would  appear  to  refute  conclusion  of 
Kotel'nikov,  et  a] .8  that  law  changes  from  one 
inferior  conjunction  to  next. 
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The  depolarized  circular  component  of  the  echoes  was  observed  for  the  first  time  at  23-cm 

wavelength  during  this  series  of  measurements.  The  complete  scattering  law  for  this  component 

could  not  be  obtained  owing  to  the  weakness  of  the  echoes.  It  was  estimated  that  the  total  power 

in  this  component  was  15  db  lower  than  that  in  the  expected  sense.  Attempts  to  look  for  a  polariz- 

3 

ing  effect  attributable  to  a  surface  layer  similar  to  that  described  by  Hagfors,  et  al.  from  radar 
observations  of  the  moon,  failed  because  of  inadequate  signal  strength  at  a  delay  where  such  an 
effect  would  be  expected.  The  cross  section  of  Venus  was  determined  on  a  number  of  occasions 
and  the  average  was  of  the  order  of  18  percent.  This  is  somewhat  higher  than  the  mean  (15 
percent)  observed  in  1964,  though  not  significantly  so. 

2.  Radar  Observations  of  Venus  at  Haystack 

9 

Karp,  et  al.  made  observations  of  Venus  at  X-band  (3.6-cm  wavelength)  using  the  West  Ford 
bistatic  radar  system.  They  made  the  surprising  discovery  that  the  cross  section  was  at  least 
a  factor  of  ten  lower  than  that  observed  at  S-band  (12.5  cm)  and  longer  wavelengths.  The  most 
obvious  interpretation  of  this  result  would  be  that  attenuation  of  the  X-band  signals  occurs  in 
traversing  the  atmosphere  of  Venus.  However,  the  spectrum  of  the  signals  obtained  by  Karp, 
et  al.  was  wider  than  would  be  expected  if  the  signals  were  reflected  from  the  solid  surface,  so 
that  a  combination  of  atmospheric  absorption  together  with  a  scattering  agent  in  the  atmosphere 
(rain  or  ice)  seemed  to  be  required  to  explain  these  early  results. 

CW  radar  measurements  of  the  cross  section  of  Venus  have  been  made  at  3.8-cm  wavelength 
using  the  Haystack  radar  system  on  the  days  listed  in  Table  II-4.  The  results  shown  in  this 
table  amply  confirm  the  low  value  (~1  percent)  obtained  by  Karp,  et  al.,  although  some  fluctua¬ 
tion  in  the  cross  section  has  been  observed,  commencing  about  9  weeks  after  inferior  conjunc¬ 
tion.  The  depolarized  circular  component  of  the  echoes  (Table  II-4)  is  found  to  be  ~14db  weaker 
than  the  expected  component  —  in  fair  agreement  with  the  ratio  observed  at  longer  wavelengths. 

The  spectrum  of  the  reflected  signals  has  been  determined  by  processing  in  the  CDC  3200 
computer  signal  samples  taken  at  a  rate  of  500  Hz.  Figure  11-12  gives  an  example  of  such  a 
spectrum.  On  the  basis  of  these  results  it  can  be  stated  that  the  echoes  are  reflected  from  the 
solid  surface  for  the  following  reasons:  (a)  The  center  Doppler  shift  of  the  echo  corresponds 
precisely  to  that  estimated  from  simultaneous  measurements  at  Millstone.  It  is  inconceivable 
that  an  atmospheric  scattering  agent  could  be  supported  above  the  surface  in  such  a  way  that 
the  height  with  respect  to  the  surface  changes  at  less  than  ±5  cm/sec  (the  velocity  accuracy  of 
the  measurements),  (b)  The  spectral  width  of  the  signals  corresponds  in  all  cases  to  that  pre¬ 
dicted  on  the  assumption  that  the  echoes  are  reflected  only  from  the  surface.  This  predicted 
Doppler  width  (center  to  limb)  is  listed  in  Table  II-4  as  fmax  and  indicated  in  Fig.  11-12  by 
arrows,  (c)  Some  of  the  fine  structure  in  the  spectra  (Fig.  11-12)  can  be  attributed  to  locally 
rough  regions  on  the  surface  of  Venus.  Indeed,  it  is  possible  to  identify  the  same  features  in 

successive  spectra  shifted  in  position  owing  to  the  apparent  rotation  of  the  planet.  In  addition, 

10 

some  of  these  features  are  apparently  the  same  as  those  observed  in  19  64  by  Goldstein  and 
1 1 

Carpenter. 

As  a  conclusive  test  that  echoes  are  reflected  from  the  planetary  surface,  short-pulse 
range  determinations  were  made  nearly  simultaneously  at  Millstone  and  Haystack.  For  these 
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TABLE  II- 

-4 

CW  OBSERVATIONS  OF  VENUS 

Date  (1966) 

f  (Hz) 

max 

T«h„  «”K) 

b  (Hz) 

Correction 
for  Terrestrial 
Atmospheric 
Attenuation 
<*) 

Cross 
Section 
(%t ra^) 

(a)  Polarized  Component 

18  Jan. 

65.7 

400.0  ±0.4 

130.85 

0 

1.44 

1  Feb. 

61.0 

325.63  ±0.5 

121.08 

0 

1.10 

9  Feb. 

76.9 

213.55  ±  0. 3 

152.33 

0.2 

1.34 

15  Feb. 

91.6 

105.25  ±  0.51 

181.6 

1.2 

1.11 

21  March 

140.8 

8.89  ±0.54 

289. 08 

0 

1.22 

30  March 

146.5 

9.28  ±0.16 

289.08 

0 

2.08 

15  April 

154.3 

4.51  ±0.25 

304.69 

0 

2.92 

5  May 

161.4 

1.28  ±0.24 

304.69 

0 

1.33 

17  May 

164.5 

1.39  ±  0.21 

304.69 

0 

2.29 

28  May 

166.0 

1.02  ±  0.33 

304.69 

1.3 

3.28 

15  June 

165.6 

1.23  ±  .09 

117.19 

0 

1.56 

27  June 

163.2 

0.83  ±  .09 

117.19 

0 

1.29 

(b)  Depolarized  Component 

19  Jan. 

64.0 

15.94  ±0.41 

132.82 

0 

0.057 

4  Feb. 

65.9 

12.63  ±0.40 

132.82 

0 

0.052 
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FREQUENCY  RELATIVE  TO  PREDICTED  DOPPLER  SHIFT  (Hz) 

Fig.  11-12.  Venus  power  spectrum  observed  at  3.  8-cm  wavelength 
on  1  February  1966.  Expected  position  for  reflections  from  limbs 
have  been  indicated  by  arrows.  Note  bumps  in  spectrum  which  are 
thought  to  be  reflections  from  anomalously  rough  regions  on  surface 
of  Venus. 

the  same  phase  code  modulation  scheme  was  applied  to  the  transmitted  pulses  and  the  echoes 
were  "compressed,"  using  the  same  matched-filter  delay  line.  For  observations  on  24  January 
and  4  February  when  the  shortest  pulses  (40  psec)  were  employed,  the  residuals  with  respect 
to  a  common  ephemeris  agreed  to  within  5  psec,  although  the  accuracy  of  each  measurement 
was  only  approximately  ±10  psec.  Further,  the  shapes  of  the  echoes  at  the  two  wavelengths 
were  nearly  identical,  providing  separate  evidence  that  the  3. 8-cm  echo  is  from  the  planetary 
surface  (Fig.  11-13). 

Since  the  echoes  are  reflected  from  the  surface  we  must  attribute  the  low  cross  section  to 
the  effects  of  atmospheric  attenuation,  a  low  intrinsic  surface  reflectivity  or  a  combination  of 
the  two.  A  low  surface  reflectivity  would  require  a  layer  of  lossy  material  overlying  the  surface, 
which  is  not  more  than  about  one  wavelength  thick.  While  one  cannot  exclude  this  explanation, 
it  does  seem  artificial,  and  hence  a  search  was  conducted  for  the  limb  darkening  that  should  be 
present  if  atmospheric  attenuation  is  present.  Two  approaches  were  adopted.  In  the  first  the 
average  of  the  four  best  signal  spectra  was  obtained  and  transformed  (using  the  Bessel  tranforma- 

tion  routine  of  T.  Hagfors)  to  yield  the  angular  power  spectrum  P(<p)Y  of  the  X-band  signals. 

^12 

As  a  comparison  law  the  mean  S-band  spectrum  reported  by  Muhleman  was  similarly  trans¬ 
formed  to  yield  the  angular  scattering  law  P(<p)g  at  S-band.  It  was  then  assumed  that  in  the 
absence  of  atmospheric  attenuation  these  two  would  be  the  same,  i.e.,  that  the  scattering  law 
does  not  change  rapidly  with  wavelength.  In  the  case  of  the  moon  this  assumption  would  not 
hold,  but  Venus  appears  to  be  considerably  smoother  than  the  moon,  and  no  clear-cut  wavelength 
dependence  has  been  detected  in  the  scattering  law  thus  far.  If  the  one-way  zenithal  attenuation 
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Fig.  11-13.  Plots  of  echo  power  vs  delay  observed 
at  23-  and  3.  8-cm  wavelengths  on  24  January  1966 
when  two  radars  were  modulated  using  same  coded 
pulse  waveform,  and  reflected  signals  processed  in 
identical  fashion.  Absolute  delays  at  two  frequen¬ 
cies  were  found  to  agree  to  within  experimental 
error,  and  echo  shapes  are  strikingly  alike.  This 
confirms  that  the  surface  is  responsible  for  scattering 
at  both  wavelengths. 


in  the  atmosphere  of  Venus  is  A  db  we  expect  the  two  laws  to  differ  according  to 

10  l°g10P<^S  ~  10  lo^iOP{^X  =  ZAlseC(P  -  1]  db  .  (1) 

This  expression  assumes  rectilinear  propagation  and  that  the  attenuating  layer  has  a  height 
extent  that  is  very  small  compared  with  the  planetary  radius.  The  results  obtained  for  the  dif¬ 
ference  did  exhibit  a  sec  (p  dependence,  and  they  suggest  a  value  A  =  1.5  —  2.0  db. 

In  the  second  approach  the  angular  scattering  law  P(</>)^  was  determined  directly  from  a 
series  of  short-pulse  experiments.  Pulses  of  40,  360,  500  psec,  1.0  and  4.0  msec,  were  em¬ 
ployed  to  construct  the  mean  echo  power  vs  delay  curve  P(t)x.  This  was  then  compared  with 
the  function  obtained  at  Millstone  in  a  similar  manner.  The  uncertainty  in  P(t)^  is  somewhat 
larger  than  that  in  the  mean  X-band  spectrum,  because  of  the  reduced  average  power  achieved 
when  operating  the  Haystack  radar  in  a  pulse  mode.  A  comparison  of  the  X-  and  L-band  scat¬ 
tering  laws  in  the  limb  region  indicated  a  value  A  =  2.5  —  3.0  db.  The  difference  between  this 
result  and  that  obtained  from  an  examination  of  the  spectra  may  stem  from  a  difference  between 
P(0)  at  23-  and  at  12.5-cm  wavelength,  or  may  be  simply  due  to  experimental  error.  In  any  event 
we  conclude  that  atmospheric  attenuation  is  present  but  that  the  zenithal  absorption  is  only  of 
the  order  of  2  db  one  way,  and  not  5  db  as  anticipated  from  the  cross-section  measurements. 

It  should  be  pointed  out  that  the  above  tests  for  atmospheric  attenuation  depend  upon  the 
scattering  behavior  near  the  limbs,  whereas  the  largest  part  of  the  echo  energy  is  reflected 
from  the  center  of  the  disk.  Thus  a  variation  of  the  attenuation  over  the  face  of  the  planet  might 
account  for  the  difference  between  the  value  of  A  inferred  from  the  observed  limb  darkening  and 
the  low  cross  section. 
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The  variation  in  cross  section  observed  with  date  (Table  II-4)  supports  the  view  that  the 

attenuation  may  vary  over  the  face  of  the  disk.  It  would  then  seem  that  the  absorbing  agent  may 

be  a  constituent  of  the  atmosphere  (like  water  vapor)  whose  abundance  may  be  controlled  by 

1 3 

large  scale  weather  patterns.  Plummer  and  Strong  have  argued  for  such  a  model  on  the  basis 
of  a  number  of  lines  of  reasoning.  However,  it  may  still  be  possible  to  reconcile  a  dense 
COz-N2  model  atmosphere  with  the  observations  if  the  height  of  the  terrain  near  the  subradar 
point  with  respect  to  the  mean  surface  varies  by  several  kilometers,  since  the  absorption  is 
expected  to  vary  with  the  square  of  the  atmospheric  pressure. 

J.  V.  Evans 
R.  P.  Ingalls 


REFERENCES 


1.  J.  V.  Evans  and  T.  Hagfors,  "Radio  Echo  Studies  of  the  Moon  at  23  cm  Wavelength." 

J.  Geophys.  Res.  (in  press). 

2.  J.  V.  Evans  and  G.H.  Pettengill,  "The  Scattering  Behavior  of  the  Moon  at  Wavelengths 
of  3.6,  68  and  784  Centimeters,"  J.  Geophys.  Res.  68,  423-447  (1963). 

3.  T.  Hagfors,  R.A.  Brockelman,  H.H.  Danforth,  L.B.  Hanson  andG.M.  Hyde.  "Tenuous 
Surface  Layer  on  the  Moon:  Evidence  Derived  from  Radar  Observations,"  Science 
150,  1153-1156  (1965). 

4.  J.  V.  Evans,  "Radio-Echo  Observations  of  the  Moon  at  3.6-cm  Wavelength"  Technical 
Report  256,  Lincoln  Laboratory,  M.I.T.  (19  February  1962),  DDC  274669. 

5.  M.  E.  Ash,  W.  B.  Smith  and  1. 1.  Shapiro,  private  communication. 

6.  W.  B.  Smith,  1. 1.  Shapiro  and  M.  E.  Ash,  private  communication. 

7.  J.  V.  Evans,  et  al.,  Astron.  J.  70,  486  (1965). 

8.  V.  A.  Kotel’nikov,  et  al.,  Doklady  Akad,  Nauk  SSSR  163.50  [English  Transl.  Soviet 
Phys.  —  Doklady.  1965]. 

9.  D.  Karp,  W.  E.  Morrow  and  W.  B.  Smith,  Icarus  3,  473,  (1964),  DDC  613893. 

10.  R.  M.  Goldstein,  J.  Research,  Natl.  Bur.  Standards  69D,  1623  (1965). 

11.  R.  L.  Carpenter,  Astron.  J.  71,  142  (1966). 

12.  D.  O.  Muhleman,  J.  Research,  Natl.  Bur.  Standards  69D,  1630  (1965). 

13.  W.T.  Plummer  and  J.  Strong,  Astronautic  Acta.  LI,  365  (1965). 


19 


m.  RADIO  ASTRONOMY 


A.  DIGITAL  MAPPING  PROCEDURE  APPLIED  TO  GALACTIC  CENTER 

The  problem  of  producing  contour  maps  from  measurements  of  antenna  temperature  with  a 
radio  telescope  can,  in  principle,  be  handled  automatically  with  a  digital  computer.  The  pres¬ 
ent  program  is  designed  to  produce  contour  plots  from  radiometric  drift  scans  in  as  much 

A 

detail  as  the  antenna  resolution  permits.  The  problem  is  to  obtain  a  record  of  the  observa¬ 
tions  (time,  antenna  coordinates,  and  observed  antenna  temperatures)  in  a  form  which  can  con¬ 
veniently  be  read  into  a  computer.  The  Haystack  radiometric  system,  described  in  the  last 
report2  produces  these  records  on  magnetic  tape.  The  mapping  program  uses  these  particular 
tapes  as  input,  but  can  be  adapted  to  handle  similar  data  taken  on  other  radio  telescopes  in  the 
form  of  a  series  of  drift  scans  across  the  source,  spaced  one-half  beamwidth  apart.  An  IBM 
7094  computer  and  a  Stromberg-Carlson  4020  display  device  were  used  for  data  processing. 

The  region  mapped  in  the  first  application  of  this  computer  program  was  centered  on  the 
radio  source  Sagittarius  A,  believed  to  be  the  center  of  our  galaxy.  The  data  were  taken  jointly 

with  Dr.  Alan  Maxwell  and  Mr.  Dennis  Downes  of  Harvard  University,  and  manually  drawn  maps 

3 

of  the  central  region  have  been  previously  published.  The  present  map  was  prepared  by  computer 
from  the  data  taken  with  the  8-GHz  Haystack  radiometer.  This  radiometer  consists  of  two  chan¬ 
nels,  each  500  MHz  wide,  centered  at  7750  and  8250  MHz.  For  a  1-sec  integration  period,  the 
root-mean- square  fluctuations  are  about  0.1  °K  for  each  channel.  All  the  data  were  taken  in  drift 
scans  with  the  observed  antenna  temperatures  consisting  of  2-sec  integrals  of  the  radiometer 
output  with  the  recorded  antenna  pointing  in  right  ascension  a,  and  declination  <5,  corresponding 

to  the  center  of  each  integration  interval.  The  details  of  the  on-line  radiometric  data  system 

4 

are  described  in  a  previous  report.  The  half-power  beamwidth  of  the  Haystack  antenna  is  4.2 
arc  minutes  at  this  frequency. 

The  mapping  program  consists  of  three  parts:  Program  A  which  produces  smoothed  plots 
of  each  drift  scan,  showing  each  radiometer  channel  separately;  Program  B  which  compiles 
coincident  drift  scans  including  one  or  both  channels  as  specified;  and  Program  C  which  pro¬ 
duces  an  automatic  plot  of  the  contour  map.  These  programs  run  in  series,  giving  the  observer 
a  look  at  the  processed  data  after  each  step.  A  more  detailed  description  of  each  program  is 
given  below. 

Program  A:-  The  recorded  data  tapes  are  read  by  Program  A.  A  specified  length  of  rec¬ 
ord  at  the  beginning  and  at  the  end  of  each  scan  (intervals  greater  than  two  beamwidths  in  each 
case)  are  used  as  baseline  data  and  possible  linear  drifts  in  baseline  are  removed  from  the  cen¬ 
tral  portions  of  the  scan.  The  data  points  are  then  smoothed  and  interpolated  by  means  of  Sub¬ 
routine  Smooth,  which  fits  a  Fourier  sine  series  to  the  observed  data.  The  series  is  truncated 
at  the  highest  angular  frequency  which  the  antenna  resolution  permits.  The  number  of  Fourier 
coefficients  used  is  2©//3,  where  ©  =  (RA  interval  of  scan)  cos  6  and  0  is  the  half-power  beam- 

5 

width.  Subroutine  Smooth  produces  the  proper  smoothing  and  interpolates  the  data  to  a  stand¬ 
ard  set  of  equally  spaced  right  ascensions.  Thus  the  data  are  registered  in  preparation  for 
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Example  of  graphical  output  of  Program  A. 
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Programs  B  and  C.  Figure  III- 1  shows  an  example  of  the  graphical  output  of  Program  A.  The 
observed  data  points  are  plotted  as  well  as  the  smooth  curves  for  both  channels.  In  addition  to 
these  output  plots,  a  magnetic  tape  record  of  the  smoothed,  interpolated  data  is  also  written  for 
Program  B. 

Program  B:—  After  examination  of  the  individual  scans  above,  the  observer  may  omit  in¬ 
dividual  channel  records  if  interference  or  equipment  difficulties  are  apparent.  Program  B 
uses  the  tape  output  of  A  together  with  editing  information  on  punched  cards,  indicating  the 
scans  to  be  averaged  and  the  channel  to  be  omitted  on  specific  scans.  Program  B  makes  an 
atmosphere  attenuation  correction  to  each  scan,  averages  coincident  scans  together,  and  calls 
Subroutine  Smooth  for  the  averaged  data  for  scans  at  each  declination.  The  antenna  tempera¬ 
tures  are  corrected  for  atmospheric  attenuation  by  dividing  by  exp  [—  sec  0],  where  is  the 
zenith  attenuation  and  0  is  the  zenith  angle  of  the  antenna  beam.  The  zenith  attenuation  tq  was 
determined  from  radiometric  observations  of  the  center  of  the  moon;  a  value  tq  =  0.023  (or 
0.1  db)  was  used.  Figure  III-2  shows  an  example  of  the  graphical  output  from  this  program. 

Data  points  from  Program  A  are  also  plotted. 

Program  C:-  All  the  information  necessary  to  plot  a  contour  map  made  with  an  axially 
symmetric  antenna  of  resolution  0  is  contained  in  a  regular  array  of  antenna  temperature  values 
at  positions  spaced  at  invervals  0/2  in  right  ascension  and  declination.  Program  C  assembles 
such  an  array,  termed  the  data  matrix,  from  the  output  of  Program  B.  The  contour  curves, 
however,  are  computed  from  an  array  which  is  more  finely  subdivided  than  the  Data  Matrix. 

In  the  present  case  the  spacing  of  this  array,  the  Augmented  Matrix,  was  p/8.  The  program 
is  written  so  that  subsections  of  the  Data  Matrix  can  be  interpolated  to  obtain  Augmented  Matri¬ 
ces  which  do  not  exceed  the  memory  capacity  of  the  computer.  A  large  map  can  thus  be  as¬ 
sembled  in  sections.  Interpolation  from  the  Data  Matrix  to  the  Augmented  Matrix  is  done  by 
using  Subroutine  Smooth,  first  parallel  to  the  right  ascension  axis  and  then  parallel  to  the  decli¬ 
nation  axis.  The  contour  lines  are  determined  by  computing  the  intersection  points  of  a  given 
contour  level  with  a  rectangular  grid  which  has  intersections  at  the  coordinate  points  of  the 
Augmented  Matrix.  Linear  interpolation  is  used  to  locate  these  intersection  points  between 
pairs  of  adjacent  points  on  the  grid.  The  contour  lines  are  plotted  with  a  series  of  straight- 
line  segments  connecting  the  grid  intersection  points.  Figure  II1-3  shows  a  plot  of  the  contours 
for  the  galactic  center  region  exactly  as  plotted  on  the  SC4020  graphical  output  device. 

The  brightest  region  in  Fig.  III-3  is  the  radio  source  Sagittarius  A.  These  observations 
confirm  that  the  microwave  spectrum  of  Sagittarius  A  is  nonthermal  and  show  that  the  diam¬ 
eter  of  the  source  is  about  3.5  minutes  of  arc.  The  contours  north  of  Sagittarius  A  evidently 
correspond  to  thermal  sources  which  appear  to  be  still  unresolved  with  a  beamwidth  of 

4.2  minutes  of  arc.  M.  L.  Meeks 

J.  C.  Henry 

B.  REAL-TIME  PROGRAM  FOR  SPECTRAL  LINE  MEASUREMENTS  AT  HAYSTACK 

Microwave  spectral  analysis  is  performed  "on-line"  by  the  same  computer  which  guides 
the  Haystack  antenna,  in  addition  to  its  previous  duties  of  antenna-pointing  and  real-time  radio- 
metric  data  reduction,  the  Univac  490  computer  now  also  sums  100  channels  of  data  from  the 
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Fig.  111-3.  Contour  map  of  the  Sagittarius  A  region  at  3.7  cm  wavelength.  This  map 
was  prepared  by  the  SC4020  output  device  from  contours  computed  in  the  IBM  7094. 
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Output  of  the  spectral 


line  program  showing  antenna  temperature  vs  frequency. 
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Digital  Autocorrelator.  On  recognizing  the  end  of  the  preset  integration  period,  the  computer 
program  transforms  the  multichannel  sums,  to  produce  a  spectrum  which  is  immediately  printed 
out  on  the  high-speed  printer  as  well  as  on  magnetic  tape.  The  observer  can  select  many  op¬ 
tions,  some  via  console  switches,  others  via  typewriter  dialog. 

During  transforming  and  printing  (plotting),  the  computer  program  avoids  summing  of  fur¬ 
ther  data,  but  continues  the  antenna  pointing.  The  spectral  summations,  performed  in  a  burst 
of  200  every  l/5  second,  may  go  on  for  hours  of  integration.  Not  quite  20  seconds  are  needed 
for  the  three-page  high-speed  printout,  whereas  the  transformations  take  only  about  5  seconds 
in  spite  of  concurrent  requirements  for  pointing  the  antenna.  The  mathematical  basis  for  the 

g 

computation,  and  the  design  of  the  Digital  Autocorrelator  have  been  described  by  Weinreb. 

The  three  printed  pages  include  one  containing  the  Main  Data  Plot  (Fig.  III-4)  of  the  tem¬ 
perature  profile  T.  =  T(f).  Abscissas  are  always  displayed  both  as  frequencies  f.  and  as  equiv¬ 
alent  velocities  with  respect  to  the  local  standard  of  rest.  Ordinates,  however,  can  be  altered 
in  appearance  if  the  observer  exercises  his  option  to  vary  the  scale  and  zero  line;  he  may  also 
add  a  slope-correction  factor,  which  affects  the  plot  but  not  the  ordinate  axis  values.  Every 
other  plotted  point  is  omitted,  for  brevity,  but  all  161  points  (T. ,  f.)  are  tabulated  on  an  auxil¬ 
iary  page  of  the  printed  output.  At  the  top  of  each  page  is  printed  the  title  of  the  run,  along  with 
date  and  day  number.  The  next  3  lines  of  the  Main  Data  Plot  page  allow  comments,  the  inser¬ 
tion  of  which  is  among  the  options  available  to  the  user.  Other  options  include  presetting  some 
of  the  values  of  parameters  itemized  in  the  legend  to  the  right  of  the  Main  Data  Plot.  In  the 
legend,  CAL  GMT  gives  the  time  at  which  the  user  terminated  the  last  prior  Cal-Base  run 
sequence. 

Console  switches,  and  typewriter  dialog,  give  the  observer  control  over  the  program.  The 
console  switches  include  the  integration-duration  dials  (values  000  to  999),  the  4-position  mode 
switch  (STOP,  CAL,  BASE,  OBS)  with  which  the  user  performs  a  Cal-Base  sequence  or  starts 
a  spectral  Observe  run,  and  eight  toggle  switches  which  trigger  such  actions  as  replotting,  re¬ 
transformation  or  reinitializing  (including  starting  a  new  typewriter  dialog). 

The  typewriter  dialog  is  automatically  activated  when  the  AUTOCOR  program  is  first  called 
in  from  the  system  tape:  after  initializing  itself,  this  program  identifies  itself  and  types  out 
(Fig.  Ill—  5 )  first  a  summary  of  toggle-switch  meanings  and  then  a  list  of  some  1 2  options  which  can 
be  chosen  via  typewriter.  The  user  replies  by  typing  a  number  followed  by  a  carriage  return(*  ). 

Several  choices  result  in  further  dialog.  For  instance,  if  the  user  types  1^  f°r  "Comments," 
the  program  further  inquires  which  line  (1-3)  he  desires  to  replace,  and  on  receiving  his  reply 
asks  for  the  text  of  the  comment;  it  then  asks  for  another  line  choice,  etc.  (Line  choice  0_  ter¬ 
minates  the  "Comments”  dialog  and  returns  to  the  "main"  dialog.)  If  the  user  chooses  8^  for 
"Slope,"  the  program  types  out  the  current  slope-correction  factor  and  asks  for  the  new  value; 
similarly,  in  choice  11*  ("Plot-scale")  and  1 2"  ("Zero  level"),  which  have  been  mentioned  pre¬ 
viously  as  user-operations  which  change  the  appearance  of  the  Main  Data  Plot.  Replying  0^  to 
the  "main"  dialog  terminates  it,  whereupon  the  program  leaves  the  initialization  section  and 
enters  the  working  section,  allowing  the  observer  to  start  a  run. 

In  summary,  the  spectral  analysis  program  used  at  Haystack  is  characterized  by:  real¬ 
time  integration  of  100-channel  spectral  data;  transforming  and  plotting  in  quasi-realtime  (data 
integration  terminated,  but  antenna  pointing  continues);  quick-look  capability  implicit  in  the 
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Fig.  111-5.  Sample  of  typewriter  dialogue  for  control  of  spectral  line  program. 

immediate  plotting,  and  especially  in  unlimited  replotting  with  parameters  altered  at  the  will 

of  the  observer;  magnetic-tape  logging  of  both  multi-channel  sums  and  transformed  spectrum 

(plus  auxiliary  parameters  of  the  run),  for  later  processing  by  other  methods,  and  for  pen- 

plotter  automatic  graphical  display  for  publication.  ^  Conaut 

Consultant 
M.  L.  Meeks 

C.  STOKES  PARAMETERS  FOR  1665-MHz  OH  EMISSION  NEAR  W3 
1.  Introduction 

The  anomalous  1665-MHz  spectral  line  emission  from  OH  near  W3  has  been  observed  to  be 
partially  linearly  polarized  by  Weinreb,  Meeks,  Carter,  Barrett  and  Rogers  9  Subsequent  ob¬ 
servations  of  right  and  left  circular  polarization  from  this  source  by  Barrett  and  Rogers10  re¬ 
vealed  that  all  four  of  the  OH  transitions  at  1612,  1665,  1667  and  1720  MHz  produce  features  in 
the  emission  spectrum  with  predominantly  circular  polarization. 

Linear  and  circular  polarization  in  the  1665-MHz  emission  from  the  W3  source  have  also 

1 1 

been  observed  by  Davies,  de  Jager  and  Verschuur.  This  article  presents  a  measurement  of 
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the  Stokes  parameters  as  a  function  of  frequency  for  the  1665-MHz  OH  emission  near  W3.  The 
polarization  parameters  -  ellipticity,  ellipse  position  angle,  and  degree  of  polarization  -  are 
computed  and  their  rapid  variation  with  frequency  is  described  in  terms  of  a  simple  model. 


2.  Polarization  Parameters 

The  polarization  properties  of  an  electromagnetic  wave  are  completely  specified  by  the  four 
Stokes  parameters,  which  may  be  determined  from  measurements  of  the  intensities  of  the  wave 
with  four  independent  polarizations.  If  I (<p)  is  the  intensity  for  linear  polarization  at  a  position 
angle  <p,  and  if  I(R)  and  I(L)  are  the  intensities  for  right  and  left  circular  polarizations,*"  then 
the  Stokes  parameters  are  defined  as 

SQ  =  1(0°)  +  1(90°)  =  I(R)  +  I(L) 

Sd  =  1(0°)  -  1(90°) 

S2  +  1(45°)  -  1(135°) 

S3=I(R)-I(L)  .  (1) 


These  parameters  have  the  property  that  the  polarization  of  an  incoherent  superposition  of  waves 
is  described  by  Stokes  parameters  which  are  just  the  sums  of  respective  parameters  for  the  sep¬ 
arate  waves.  Other  parameters  describing  the  polarization  can  be  obtained  from  the  Stokes  pa¬ 
rameters  as  follows.  The  position  angle  </?  of  the  major  axis  of  the  polarization  ellipse  is  given 
by 

tan  (2 (p)  =  S2/S1  .  (2) 


The  amplitude  ratio  of  minor-to-major  axes  b/a  for  the  polarization  ellipse  is  expressed  as 


tan  x  =  ±b/a 


(3) 


where 


sin  2X  =  S3  Jsf  +  S22  +  S32  ,  (4) 

and  the  +  in  Eq.(3)  applies  to  right  elliptical  and  the  -  to  left  elliptical  polarization.  Finally  the 
degree  of  polarization  P  is  obtained  from 

P  =  Jsf  +  S22  +  S3/sQ  (5) 

It  is  convenient  to  work  with  the  Stokes  parameters  because  of  their  additive  property,  but  some 
what  better  physical  insight  can  be  obtained  from  the  ellipse  position  angle  <p,  the  ellipticity 
±b/a,  and  the  degree  of  polarization  P,  as  defined  above. 


3.  Spectral  Line  Measurements 

Measurements  of  the  polarization  parameters  for  the  1665-MHz  OH  emission  were  made 
with  the  1 2  0  -  ft  parabolic  antenna  of  the  Haystack  research  facility  of  Lincoln  Laboratory, 
Massachusetts  Institute  of  Technology,  during  the  period  28  November  through  9  December  1965. 

*  In  th,s  PaPer  right  circular  polarization  denotes  clockwise  rotation  of  the  electric  vector  when  viewed  along  the  direction  of 
propagation. 
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Fig.  111-6.  Polarization  parameters  for  1665-MHz  OH  emission  near  W3.  (a)  and  (b)  show  measured 

values  of  Stokes  parameters.  Points  plotted  in(c),  (d)  and  (e)  were  computed  from  Stokes  parameters, 
and  solid  lines  represent  values  derived  from  model  from  -46  to  -44  km/sec,  or  otherwise  best  esti¬ 
mates  of  parameters. 
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The  observation  of  various  linear  and  circular  polarizations  with  a  single  system  made  possible 
a  complete  determination  of  the  polarization  parameters  of  the  spectral  features  as  a  function  of 
frequency,  or  equivalently  as  a  function  of  radial  velocity  with  respect  to  the  local  standard  of 
rest.  At  this  OH-transition  frequency  a  radial  velocity  of  +1  km/sec  gives  a  Doppler  shift  of 
-  5.555  kHz. 

The  spectral  line  radiometer  used  in  these  measurements  included  a  100-channel  digital 
autocorrelator  and  a  room-temperature  parametric  amplifier.  The  autocorrelator  is  tied  di¬ 
rectly  to  a  real-time  computer  which  produces  plots  of  the  observed  spectra  immediately  after 
each  integration  period.  Most  of  the  observations  were  made  with  a  1-kHz  frequency  resolution, 
but  the  brightest  features  were  also  examined  with  200-Hz  resolution.  The  overall  system  tem¬ 
perature  was  320°K  for  linear  polarization  measurements  and  360°K  for  circular  polarization; 
the  increased  system  temperature  resulted  from  a  network  that  combined  horizontal  and  vertical 
linear  polarization  to  produce  the  circular  polarizations.  The  observed  antenna  temperatures 
for  Cassiopeia  A  continuum  emission  were  240°K  and  180°K,  for  linear  and  circular  polarization 
respectively.  The  isolation  of  right  and  left  circular  polarization  was  greater  than  15db,  and 
the  half-power  beamwidth  was  22  minutes  of  arc. 

Figure  1 1 1  -  6  (a )  and  (b)  shows  the  measured  Stokes  parameters  plotted  as  a  function  of  radial 
velocity  with  the  RMS  errors  and  frequency  resolution  specified.  In  Fig.  1 1 1  —  6  (b )  the  antenna  tem¬ 
perature  scale  has  been  expanded  by  a  factor  of  two  over  that  of  Fig.  Ill  -  6  (a ) .  It  is  clear  from 
these  curves  that  the  polarization  is  primarily  circular,  progressing  from  totally  right  circular 
on  the  right  side  of  Fig.  Ill  —  6 (a )  to  totally  left  circular  for  the  feature  just  to  the  left  of  the  cen¬ 
tral  maximum.  Figure  III- 6 (b )  shows  that  the  amount  of  linear  polarization  is  relatively  small 
by  comparison.  The  polarization  behavior  across  the  central  maximum  is  evidently  very  complex. 

4.  Interpretation  of  Polarization 

Once  the  Stokes  parameters  have  been  determined,  it  is  desirable  to  use  them  to  compute 
the  ellipticity,  the  ellipse  position  angle,  and  the  degree  of  polarization  for  the  emission  spec¬ 
trum.  These  last  three  parameters,  however,  are  not  proportional  to  the  received  signal  and 
can  be  determined  only  when  the  signal-to-noise  ratio  is  sufficiently  high. 

Equations  (3)  through  (5)  were  used  to  compute  the  ellipticity  ±b/a  and  the  degree  of  polariza¬ 
tion  P.  The  ellipse  position  angle  </?  was  determined  as  a  parameter  in  the  least-square  fitting 
procedure,  Eq.(8).  The  root-mean-square  (RMS)  errors  associated  with  these  three  parameters 
were  computed  by  estimating  the  RMS  errors  in  the  Stokes  parameters  and  assuming  linear  error- 
propagation.  The  polarization  parameters  and  the  RMS  error  bars  are  shown  in  Fig.  1 1 1- 6  (c ),  (d), 
and  (e).  The  plotted  points  here  are  separated  by  an  equivalent  filter  width  in  frequency;  so  they 
represent  independent  measurements.  The  large  errors  are  associated  with  low  signal-to-noise 
ratios,  and  no  points  are  plotted  when  the  errors  become  excessive. 

The  polarization  parameters  in  F ig.  Ill- 6  show  continuous  variation  around  the  central  maxi¬ 
mum,  but  outside  this  region  of  variability,  the  polarization  parameters  appear  to  remain  con¬ 
stant,  within  experimental  error,  across  the  individual  emission  features.  We  have  used  a 
simple  model  to  describe  the  variability  of  the  polarization  parameters  in  the  region  of  the  cen¬ 
tral  maximum.  This  model  assumes  that  there  are  three  incoherent  emission  features  in  the 
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radial  velocity  interval -46  to— 44  km/sec.  These  features  are  assumed  to  be  100-percent  po¬ 
larized  and  to  have  Gaussian  profiles.  Polarization  variability  then  results  from  incoherent 
superposition.  The  problem  of  fitting  this  model  to  the  observed  polarization  parameters  in 
the  region  of  the  central  maximum  is  not  a  simple  one.  Five  parameters  are  necessary  to  spec¬ 
ify  the  center  frequency,  width,  peak  intensity  and  polarization  of  each  feature.  These  param¬ 
eters  were  obtained  by  repeated  trial  and  error  with  a  computer  program  that  combined  trial 
values  of  the  Stokes  parameters  in  this  model  and  plotted  the  resulting  polarization  parameters 
with  the  observed  data  superimposed.  The  solid  curves  in  Fig.  Ill  —6 (c ) ,  (d)  and  (e)  show  the  best 
fit  that  was  obtained  by  this  procedure.  The  parameters  that  produced  this  fit  are  given  in 
Table  III-l  which  also  includes  the  corresponding  parameters  for  all  features  in  the  spectrum 
shown  in  Fig.  III-6.  The  overall  agreement  is  good  between  the  observed  data  and  the  model  in 
Fig.  III-6(c),  (d)  and  (e).  However,  there  are  small  discrepancies  that  nevertheless  may  be 
large  compared  with  the  RMS  error  bars  on  the  data  points.  We  believe  that  these  discrepancies 
result  from  the  following  causes:  (a)  the  features  are  not  exactly  Gaussian  in  shape;  (b)  the  res¬ 
olution  of  1  kHz  is  not  quite  sufficient  to  completely  resolve  the  spectrum  of  the  source;  (c)  the 
cut-and-try  fitting  procedure  did  not  produce  the  best-fit  parameters.  In  general,  this  simple 
model  does  indeed  appear  to  fit  the  observed  variations  in  polarization  parameters. 

Outside  the  interval -44  to-46km/sec  in  which  an  overlap  of  100-percent  polarized  features 
is  assumed  to  account  for  P  <  1,  only  two  of  the  remaining  five  features  are  significantly  less 
than  100-percent  polarized.  The  feature  at  —  49. 1  km/sec  is  51-percent  polarized  (±3  RMS  error) 

and  the  features  at  —  43.7  km/sec  is  86-percent  polarized  (±6  RMS  error). 

1 1 

It  has  been  proposed  by  Davies,  de  Jager  and  Verschuur  that  the  Zeeman  effect  might  ac¬ 
count  for  this  spectrum.  However,  the  identification  of  Zeeman  components  suggested  by  these 
authors  does  not  appear  consistent  with  our  observed  polarization  parameters  and  the  results  in 
Table  III- 1 .  It  is,  of  course,  impossible  to  say  on  the  basis  of  our  measurement  that  no  com¬ 
bination  of  simple  Zeeman  components  can  produce  the  observed  spectrum  for  the  1665-MHz 

transition.  Rut  the  circular  polarization  observations  of  the  other  three  OH  transitions  from 

1  0 

the  W3  source  by  Barrett  and  Rogers  appear  to  rule  out  a  simple  interpretation  of  this  1665- 
MHz  spectrum  in  terms  of  Zeeman  splitting. 

5.  Conclusions 

The  spectrum  of  the  1665-MHz  OH  transition  observed  near  W3  shows  the  bright,  narrow 
emission  features  that  are  characteristic  of  the  anomalous  emission  associated  with  interstellar 
gas  clouds  near  H  II  regions.  A  detailed  examination  of  the  polarization  parameters  shows  that 
the  spectrum  can  be  interpreted  as  an  incoherent  superposition  of  emission  features,  each  with 
polarization  parameters  that  are  constant  across  the  feature.  The  features  are  roughly  Gaussian 
in  shape,  but  some  significant  deviations  from  Gaussian  shape  are  observed.  Fight  features  ap¬ 
pear  to  be  present  in  the  spectrum,  but  some  additional  low-level  emission  may  be  present,  par¬ 
ticularly  in  the  interval— 49  to  — 47  km/sec.  Three  emission  features  are  observed  to  be  100- 
percent  polarized  within  the  observational  errors.  Three  other  features,  which  overlap  in  fre¬ 
quency,  are  consistent  with  a  model  in  which  the  features  themselves  are  100-percent  polarized 
and  the  observed  degree  of  polarization  is  a  consequence  of  an  incoherent  superposition.  The 
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TABLE  MI-2 

RADIOMETRIC  SYSTEM  PARAMETERS 


Quantity 

Center  Frequency  (GHz) 

Units 

8.25 

15.50 

Bandwidth 

500 

1000 

MHz 

Noise  temperature 

1300 

2150 

°K 

Sensitivity  for  1-sec  integration 

0.10 

0.16 

°K 

Fig.  111-7.  Histogram  of  percent  gain  change 
in  40-minute  time  interval  for  2-cm  radiometer. 
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two  remaining  features  are  not  obviously  overlapped,  yet  give  P  <  1.  However,  some  overlap 
with  low-level  emission  may  be  present  in  these  two  cases.  Generally  the  observations  are 
consistent  with  a  superposition  of  100-percent  polarized  features  that  are  mutually  incoherent. 
The  observed  polarization  varies  continuously  from  left  circular  through  linear  to  right  elliptical 
and  back  to  left  elliptical  between  radial  velocities  of —  46  and  —44  km/sec.  The  model  simu¬ 
lates  this  behavior  with  features  that  are  elliptical  with  ellipticity  |b/a|  >0.65  or  equivalently 
with  linear  polarization  amounting  to  41  percent  or  less. 

The  mechanism  by  which  these  features  are  produced  is  not  clearly  understood.  However, 
the  description  of  the  emission  in  terms  of  our  simple  model  makes  it  possible  to  conclude  that 
the  1665-MHz  emission  near  W3  is  not  qualitatively  different  from  that  of  the  other  three  transi¬ 
tions,  notwithstanding  the  observed  variations  in  the  polarization  parameters  as  a  function  of 

radial  velocity.  „ „  _ 

J  M.  L.  Meeks 

J.  A.  Ball 
J.  C.  Carter 

D.  ABSOLUTE  FLUX  MEASUREMENTS  OF  CASSIOPEIA  A  AND  TAURUS  A 
AT  3.64  AND  1.94  cm 

1.  Introduction 

Measurements  of  the  absolute  intensity  of  the  radiation  received  at  the  earth  from  a  celes¬ 
tial  radio  source  are  important  for  two  reasons:  such  measurements  at  different  wavelengths 
yield  information  about  the  physical  processes  occurring  in  the  object  itself,  and  they  provide 
a  method  of  measuring  the  gain  of  large  antennas  which  is  free  of  many  of  the  spurious  effects 
associated  with  the  use  of  pattern  transmitters. 

The  requirements  for  absolute  flux  measurements  at  centimeter  wavelengths  are:  sensitive, 
stable,  accurately  calibrated  radiometers,  and  an  antenna  of  known  aperture  efficiency. 

2.  Radiometric  System 

At  the  Haystack  Research  Facility,  the  first  requirement  previously  mentioned  has  been 
admirably  met  with  equipment  designed  by  Dr.  Sander  Weinreb.'  Briefly,  the  radiometers  at 
both  frequencies  are  simple  Dicke-switched  TRF  receivers,  employing  a  cascade  of  wide-band 
tunnel  diode  amplifiers  followed  by  a  traveling- wave  tube.  The  signal  at  the  output  of  the  final 
RF  amplifier  is  split  into  two  500-MHz  bands  (as  an  aid  in  identifying  interference)  and  delivered 
first  to  separate  square-law  detectors,  then  synchronously  demodulated  and  applied  to  linear 
integrators.  The  500-MHz  band  at  7.7 5  GHz  is  frequently  hampered  by  interference;  the  two 
channels  at  15.25  and  15.75  GHz  are  averaged  after  final  processing.  Table  I II- 2  lists  the  rel¬ 
evant  parameters  for  the  radiometric  system. 

The  great  stability  of  these  radiometers  is  exhibited  in  Fig.  Ill- 7 .  The  overall  system  has 
been  examined  for  linearity  by  means  of  standard  techniques  (gas  discharge  tube  and  precision 
attenuator)  and  correction  to  the  present  measurements  has  been  negligible  (0.6  percent  at 
8.25  GHz  and  0.1  percent  at  15.50  GHz).  The  calibration  noise  source  in  the  radiometer  has 
been  referred  to  the  temperature  difference  between  a  termination  at  room  temperature,  and 
a  similar  termination  immersed  in  liquid  nitrogen.  The  experimental  accuracy  of  this  measure¬ 
ment  is  better  than  1  percent. 

*  See  Ref.  2,  Sec.  III. 
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3.  Cornucopia  Reflector 

2 

A  horn  reflector  antenna  with  approximately  6-m  aperture  was  obtained  for  observations 

of  the  strong  radio  sources  Cassiopeia  A  and  Taurus  A.  Measurements  of  the  aperture  efficiency 

of  this  antenna  from  2.7  to  16.0 GHz  were  performed  by  A.  Sotiropoulos,  utilizing  the  substitution 

method  on  a  simulated  free  space  range.  The  data  have  been  analyzed  by  Sotiropoulos  and 

13 

Dr.  J.  Ruze,  who  have  published  these  results  in  a  technical  report.  The  conclusions  pertinent 
to  the  present  observations  are  presented  below. 


Frequency  (GHz)  r\ 


P 


8.25 

15.25 

15.75 


0.719  ±  0.028  633  ±  25 

0.583  ±  0.023  780  ±  31 

0.572  ±  0.023  795  ±  31 


where  rj  =  aperture  efficiency 


P 


2k/n  Ag 


is  the  number  that  converts  antenna 
temperature  from  an  unpolarized  point 
source  to  flux  in  units  of  10“26  watts. 
m2.  Hz*l  according  to  =  pr T  . 


(1) 


This  antenna  is  mounted  as  a  transit  instrument  movable  in  elevation  along  the  local  me¬ 
ridian,  and  connected  to  the  radiometer  through  a  rotary  joint  and  a  short  piece  of  flexible  wave¬ 
guide.  The  loss-VSWR  combination  of  this  additional  plumbing  has  been  measured  across  the 
bandwidth  of  the  radiometers  and  introduces  losses  of  4.9  percent  at  8.25 GHz  and  6.2  percent 
at  15.50GHz. 


4.  Observations 

Transit  observations  of  the  two  strong  radio  sources  Cassiopeia  A  and  Taurus  A  were  made 
during  December  1965  and  January  1966.  Drift  scans  on  a  source  were  taken  at  several  slightly 
different  elevations  to  remove  any  pointing  error.  Each  drift  scan  was  corrected  for  baseline 
drift,  atmospheric  attenuation  on  that  particular  day,  scaled  to  temperature,  and  fitted  by  a 
Gaussian  curve  with  a  least-squares  iteration  technique.  A  sample  record  is  shown  in  Fig.  III-8. 
The  resultant  peak  temperatures  at  different  elevation  offsets  were  again  fitted  by  a  Gaussian 
curve,  and  the  final  peak  temperature  obtained.  The  data  shown  in  Fig.  II I- 9  illustrated  the 
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Fig.  111-8.  Drift  scan  of  Taurus  A  at  2  cm  with  Cornucopia  horn  antenna  on  15  January  1966  (region  2). 
Recording  time  constant  is  10  sec  and  time  lapse  from  region  1  to  region  3  is  about  40  minutes. 


♦See  Kef.  2,  p.  24. 
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Fig.  111-9.  Peak  temperatures  of  all  drift  scans  on  Cassiopeia  A  at  15.75  GHz.  Solid  curve 
is  least-squares  Gaussian  curve  through  data  points.  One  unit  on  abscissa  corresponds  to  ap¬ 
proximately  2.2  arc  minutes  of  offset  in  elevation. 


results  of  drift  scans  at  different  elevations.  The  peak  temperatures  were  corrected  for  rotary 
joint  loss  and  source  size,  and  converted  to  flux  using  Eq.(l).  (The  data  for  Taurus  A  have  also 
been  corrected  for  polarization.)  The  tabulation  below  summarizes  the  results. 


Source 


Frequency 


Flux 


Cassiopeia  A  8.25 

15.50 


616  ±  4.5%  RMS 
380  ±  6%  RMS 


Taurus  A 


8.25  567  ±  4.5%  RMS 

15.50  461  ±  6%  RMS 


The  data  indicate  that  the  spectrum  of  Cassiopeia  A  is  well  represented  by  a  constant  spec¬ 
tral  index  at  least  to  15.5  GHz,  in  contradiction  to  the  results  of  Baars,  Mezger,  and  Wendker.14 
Figure  III- 1 0  shows  the  present  observations  for  Cassiopeia  A  plotted  along  with  data  obtained 

by  other  groups  listed  in  Ref.  14.  The  data  for  Taurus  A  agree  fairly  well  with  the  lower  fre- 

a 

quency  measurements,  although  a  fit  to  S^  =  k^  for  constant  ot  is  less  convincing  than  for 
Cassiopeia  A.  Figure  III  —  1 1  shows  the  spectrum  of  Taurus  A. 


R.J  Allen 
Prof.  A.  H.  Barrett 
(Research  Laboratory 
of  Electronics,  M.I.T.) 


E.  SECULAR  VARIATIONS  IN  CENTIMETER  RADIO  FLUX  FROM  3C273 

Observations  of  the  quasi-stellar  radio  source  3C27  3  have  been  made  at  8.25  and  15. 50  GHz 
with  the  Haystack  antenna.  Data  obtained  during  the  first  few  months  of  1966  have  been  com¬ 
pared  with  similar  observations  during  the  spring  and  summer  of  1965.  The  results  at  8.25  GHz 

1  5 

indicate  a  moderate  increase  in  radio  flux  of  16  percent  per  year,  in  agreement  with  Dent. 
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Fig.  111-10.  Radio  spectrum  of  Cassiopeia  A(a).  Present  measurements  are  shown  as  open  triangles. 
Open  circles  are  previously  published  measurements  of  less  than  10-percent  quoted  error. 
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Fig.  111-11.  Radio  spectrum  of  Taurus  A. 
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Fig.  111-12.  Increase  in  centimeter  radiation  received  from  quasi-stel lar  source  3C273. 

On  the  other  hand,  observations  at  15. 5  GHz  show  a  nonlinear  increase  in  flux  by  a  factor  of  al¬ 
most  two  in  the  past  year.  On  6  April  1966,  the  flux  at  8.25  GHz  was  35.0  flux  units,'  while  at 

15.50GHz  it  was  54.5  units.  Thus  at  that  time,  3C273  has  an  upturned  spectrum  similar  to  3C84 
1 6 

and  3C279.  The  flux  from  Virgo  A  has  remained  constant  at  both  frequencies;  42.7  units  at 
8.25GHz  and  27.5  units  at  15.5  GHz.  Figure  III- 1 2  shows  the  time  dependence  of  the  radio  emis¬ 
sion  at  the  two  frequencies.  ~  T 

^  R.  J.  Allen 

Prof.  A.  H.  Barrett 
(Research  Laboratory 
of  Electronics,  M.I.T.) 


F.  OBSERVATIONS  OF  HYDROGEN  RECOMBINATION  LINES 

In  April  1966,  the  Haystack  radar  receiving  system  was  modified  to  permit  radiometric 
observations  of  the  hydrogen  recombination  line  at  7.793  GHz.  The  modifications  included 
mounting  a  Clavin  feed  next  to  the  normal  horn  feed  of  the  radar  and  adjusting  the  two  local  os¬ 
cillators  to  appropriate  frequencies.  Doppler  frequency  shift  caused  by  the  motion  of  the  Hay¬ 
stack  site  with  respect  to  the  local  standard  of  rest  was  compensated  by  tuning  the  second  local 
oscillator  during  observations.  The  parametric  amplifiers  were  cooled  by  liquid  N^  to  give  a 
system  temperature  of  approximately  200°K.  The  spectral  line  equipment  normally  used  with 
the  radiometer  receivers  was  arranged  to  process  the  signal  from  the  radar  receiver. 

Observations  were  made  of  the  92cv  hydrogen  line  in  emission  in  the  Orion  Nebula,  Omega 
Nebula  and  W3.  The  frequency  resolution  was  100  kHz  over  a  spectral  range  of  4  MHz.  Table  III- 3 
lists  the  measured  temperatures  and  line  widths  for  two  of  the  sources.  The  parameters  T^,,  Tj  , 
and  refer  to  the  increase  in  antenna  temperature  due  to  the  continuum  and  line  emission,  and 

the  half-power  width  of  the  line,  respectively. 

*  I  flux  unit  =  10’26  w  m"2  (Hz)’1 


39 


Section  III 


TABLE  III -3 

TEMPERATURES  AND  LINEWIDTHS  FOR  TWO  SOURCES 

Source 

Tc 

tl 

a,l 

Orion  Nebula 

Omega  Nebula 

21.1°  (±0.6) 

14.5°  (±0.6) 

1.5°  (±0.2) 

0.95°  (±0.1) 

850  kHz  (±100) 

1100  kHz  (±50) 

The  velocity  of  the  hydrogen  recombination  line  source  in  W3  is  -  45km/sec/±2km,  the 
same  as  that  of  the  strong  OH-line  source  in  W3.  This  coincidence  in  velocity  space  suggests 
a  spatial  association  of  the  two  emission  regions. 


M.  L.  Meeks 
M.  A.  Gordon 


G.  HAYSTACK-MILLSTONE  INTERFEROMETER 

The  Haystack  120-ft  antenna  and  the  Millstone  84-ft  antenna  have  been  used  jointly  as  an 
interferometer  in  an  attempt  to  obtain  the  positions  and  effective  sizes  of  OH  emission  regions. 
The  separation  of  2257  feet  between  the  two  antennas  produces  a  minimum  fringe-lobe  separa¬ 
tion  of  54  seconds  of  arc  at  1665  MHz  as  illustrated  in  Fig.  Ill- 1 3.  Phase  measurements  rela¬ 
tive  to  a  radio  source  of  known  position  allow  position  measurements  to  about  20  seconds  of  arc 
if  the  phase  is  measured  to  within  a  radian.  While  instrumental  phase  changes  and  the  accuracy 
of  the  baseline  limit  the  precision  of  the  absolute  measurements,  the  measurements  of  position 
of  one  feature  relative  to  another  can  be  made  with  much  higher  precision,  since  the  features 
are  separated  by  only  a  few  kilohertz. 

A  block  diagram  of  the  interferometer  is  shown  in  Fig.  Ill- 1 4.  The  radiometers  are  image- 
rejected,  tunnel-diode  amplifiers  followed  by  mixers  and  30-MHz  IF  amplifiers.  A  common 
LO  signal  is  derived  from  oscillators  phase-locked  to  the  sum  of  a  harmonic  of  67  MHz  and  an 


Fig.  111-13.  Lobe  pattern  for  Haystack-Millstone  interferometer  at  1665  MHz. 
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adjustable  frequency  around  28  MHz.  The  phase  shift  of  the  67  and  28  MHz  signals  caused  by 
line-length  changes  is  eliminated  by  a  line  servo  system  that  maintains  a  constant  electrical 
length  along  the  line  between  radiometers.  Typically,  the  line  servo  has  to  adjust  for  a  length 
change  of  6  feet  between  day  and  night  temperatures. 

The  30-MHz  signals  are  alternately  added  and  subtracted  with  a  200-psec  period.  The  auto¬ 
correlation  functions  for  the  clipped  versions  of  the  sum  and  difference  signals  are  transferred 
from  the  output  of  a  digital  correlator  to  the  Univac  490  computer.  These  autocorrelation  func¬ 
tions  together  with  antenna  command  coordinates  are  written  onto  magnetic  tape.  A  post-real¬ 
time  program  then  analyzes  the  data  for  fringe  amplitude  and  phase  information  over  the  spectrum. 
A  reference  phase  is  computed  every  200psec  from  the  relation 

co 

“■ref  =  -f  1°  sin  Eh  sin  Eg>t  +  D  cos  EH  cos  Eg(As<t  -  AR) 

-  d  cos  E  .1 
s,tJ 

Here  D,  EtT,  ATT,  and  d  are  parameters  of  the  baseline  and  A  E  .  are  the  refraction- 
H  H  r  s,t  s,t 

corrected  azimuth  and  elevation  of  a  reference  position  in  the  sky. 

Multiplication  of  the  difference  of  autocorrelation  functions  of  the  added  and  subtracted 
signals  by  the  sine  and  cosine  of  the  reference  phase  allows  a  least- squares  estimate  of  the 
fringe  amplitude  and  phase  to  be  made. 

Plots  of  the  un-normalized  fringe  amplitude  and  phase  are  shown  in  Figs.  II I  —  1  5  and  III- 1 6 

as  typical  examples  of  the  data  output.  .  ^  ^ 

r  ^  A.  E.  E.  Rogers 

(Research  Laboratory 
of  Electronics,  M.I.T.) 

G.  M.  Hyde 

H.  SPECTRUM  MEASUREMENTS  NEAR  1-cm  WAVELENGTH  USING  THE  28-ft 
MILLIMETER-WAVE  ANTENNA* 

1.  Introduction 

Several  research  programs  have  been  conducted  since  April  1964  using  the  28-ft  millimeter 

antenna  in  conjunction  with  a  five-channel  microwave  radiometer  developed  by  the  M.I.T. 

17  18 

Research  Laboratory  of  Electronics  (Staelin,  Staelin  and  Barrett  ).  A  block  diagram  of  the 
radiometer  is  shown  in  Fig.III-17.  The  major  programs  include  spectrum  studies  of  the  fol¬ 
lowing:  the  planet  Venus,  the  planet  Jupiter,  the  moon,  the  sun,  and  the  terrestrial  atmosphere. 

2.  The  Planet  Venus 

Venus  was  observed  on  15  days  from  June  to  August  1964  at  wavelengths  of  0.92,  1.02, 

I. 18,  1.28,  1.37,  and  1.42  cm  and  was  observed  on  14  days  from  January  to  March  1966  at 
wavelengths  of  1.18,  1.28,  1.35,  1.43,  and  1.58cm.  The  planetary  spectrum  was  determined 
by  measuring  the  planet-moon  ratio  and  then  correcting  the  ratio  for  the  measured  antenna 
pattern.  An  extensive  set  of  atmospheric  absorption  measurements  was  used  to  correct  the 
observed  spectra  for  atmospheric  absorption. 

*  See  Sec.  YI-B  and  Appendix  B  for  description. 
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Fig.  111-15.  Un-normal ized  fringe  amplitude,  W3  OH-emission  source,  right  circular 
polarization.  Points  on  emission  features  are  connected  with  straight  lines. 
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Fig.  111-16.  Fringe  phase,  W3  OH-emission  source,  right  circular  polarization.  Points  repre¬ 
senting  phases  of  two  emission  features  in  previous  figure  are  connected  with  straight  lines. 
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Fig.  111-17.  Five-channel  microwave  radiometer. 

This  region  of  the  microwave  spectrum  of  Venus  is  important  because  it  includes  a  rather 
sharp  transition,  the  shape  of  which  may  distinguish  various  atmospheric  models.  Certain 
atmospheric  models  that  have  been  proposed  incorporate  molecules  such  as  HzO  which  show 
absorption  at  these  wavelengths  (Barrett  and  Staelin19).  The  first  series  of  observations  indi- 
cated  the  possible  presence  of  a  time- varying  absorption  feature  near  1.35-cm  wavelength,  a 
resonant  frequency  of  water  vapor  (Staelin  and  Barrett  ). 

More  recent  observations  of  1966  were  made  with  an  improved  antenna  feed  and  radiometer. 

These  observations,  not  fully  analyzed,  presently  indicate  no  spectral  feature  at  1.35-cm  wave- 

20 

length  (Staelin  and  Neal  ). 

3.  The  Planet  Jupiter 

Jupiter  was  observed  on  five  days  from  February  to  March  1966,  using  the  same  equipment 
and  procedures  as  for  the  simultaneous  Venus  observing  program.  The  preliminary  average 
spectrum  shows  no  sharp  feature  near  1.35-cm  wavelength  (Staelin  and  Neal  ).  Resonance 
features  of  NH^  may  also  occur  in  this  wavelength  region,  and  will  be  sought  later. 

4.  The  Moon 

2 

The  moon  was  observed  principally  during  the  summer  of  1964  (Moran,  Barrett  and  Staelin; 
Moran22).  Scans  through  the  center  of  the  moon  were  used  to  determine  the  lunar  brightness 
temperature  at  the  center  of  the  lunar  disk  as  a  function  of  phase.  By  measuring  the  asymmetry 
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DETERMINATION 

TABLE 

OF  LUNAR 

MI-4 

HARMONIC  COMPONENTS 

Frequency 

Assumed  T 

o 

Ti 

91 

T2 

02 

T3 

(GHz) 

(°K) 

(°K) 

(deg) 

(°K) 

(deg) 

(°K) 

21.9 

217 

21  ±2 

41  ±5 

5  ±  1 

38  ±  5 

1 

23.5 

217 

22 

48 

6 

31 

1 

25.5 

217 

25 

42 

6 

40 

1 

32.4 

217 

29 

41 

8 

38 

1 

in  the  lunar  brightness  temperature  as  a  function  of  phase,  the  harmonic  components  were  de¬ 
termined.  The  results  are  summarized  in  Table  III-4,  where  the  lunar  brightness  temperature 
Tg  is  given  as  a  function  of  lunar  phase  (p  by 

oo 

tB  S  Ti  cos  (*♦  -  ej) 
i=0 

Measurements  of  the  difference  in  temperature  observed  when  the  antenna  feed  is  linearly 
polarized  along  and  across  the  direction  of  scan  as  a  function  of  angular  offset  from  the  center, 
as  shown  in  Fig.  Ill- 1 8,  are  consistent  with  a  lunar  dielectric  constant  of  1.5  ±  0.1  for  a  smooth 
lunar  surface,  and  1.7  for  a  surface  with  an  RMS  slope  of  15°. 

5.  The  Terrestrial  Atmosphere 

Three  major  series  of  atmospheric  observations  have  been  made:  June- August  1964,  June- 
August  1965,  and  January- March  1966.  The  wavelengths  observed  were  0.92,  1.02,  1.18,  1.28, 
and  1.37  cm;  0.92,  1.02,  1.28,  1.  35  and  1.37  cm;  1.18,  1.28,  1.  35,  1.43  and  1.58  cm,  respectively 
The  observations  consisted  principally  of  atmospheric  absorption  measurements  made  by  means 


Fig.  111-18.  Polarization  dependence  of  lunar  brightness  temperature  at  35  GHz. 
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Fig.  111-19.  Absorption  spectra  of  terrestrial 
atmosphere. 


Uj — i _ i _ i _ i _ i _ i _ i _ i _ i  i  i 

20  25  30 


FREQUENCY  (GHz) 


of  solar-extinction  observations.  Over  40  of  these  were  accompanied  by  a  radiosonde  launched 

either  from  Hanscom  Field  or  from  the  M.I.T.  campus  in  Cambridge.  The  observed  spectra 
23  24  ^ 

(Staelin,  Gant  and  Staelin  )  all  exhibit  the  1.35-cm  water  vapor  resonance,  which  varied  in 

amplitude  and  shape  depending  upon  the  water  vapor  profile  of  the  atmosphere.  Figure  III- 19 

shows  typical  observed  spectra. 


Prof.  D.  H.  Staelin 
Prof.  A.  H.  Barrett 
N.  E.  Gant 
J.  J.  Moran 
J.  J.  G.  McCue 
J.  Ft.  Cogdell 


Research  Laboratory 
of  Electronics,  M.  I.  T. 
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IV.  ATMOSPHERIC  RESEARCH 


A.  INCOHERENT  IONOSPHERE  BACKSCATTER 

i 

In  the  last  radio  physics  and  astronomy  report  we  described  in  some  detail  the  program 
of  incoherent  backscatter  observations  of  the  earth’s  ionosphere.  No  major  changes  have  been 
made  to  the  equipment  nor  to  the  methods  of  data  gathering  described  in  that  report,  though 
some  changes  have  been  made  in  the  methods  of  data  reduction. 

Observations  using  the  440-MHz  radar  have  continued  throughout  1966  for  two  periods  of 
24  hours  each,  every  month.  Measurements  are  also  made  with  the  1295-MHz  steerable  radar 
once  a  month  for  a  12 -hour  period. 

The  backlog  of  data  reduction  continues  to  be  the  most  pressing  problem.  A  second  data 
technician  has  been  acquired  to  assist  with  the  hand  analysis.  In  addition,  some  computer  pro¬ 
gram  changes  have  been  made  which  should  eliminate  one  tedious  part  of  the  hand  analysis. 
These  improvements  are,  however,  too  recent  to  have  yet  seriously  affected  the  situation.  At 
present  the  complete  analysis  of  the  1964  results  is  nearly  finished.  Partial  analysis  (compu¬ 
ter  reduction)  of  the  1965-1966  results  has  been  carried  out,  but  much  long  and  tedious  hand 
analysis  remains.  The  remedy  ultimately  lies  in  providing  the  sophisticated  computer  pro¬ 
grams  required  to  permit  much  greater  dependence  upon  machine  processing. 

J.  V.  Evans  R.  F.  Julian 

W.A.  Reid  R.A.  Brockelman 


B.  LOWER  ATMOSPHERIC  STUDIES  AT  X-BAND 

A  program  to  investigate  the  effects  of  the  lower  atmosphere  on  earth  to  space  propagation 
paths  was  initiated  in  support  of  the  Laboratory's  Space  Communications  Program.  The  initial 
area  of  investigation  is  in  the  effects  of  the  atmosphere  on  X-band,  line-of-sight  systems.  In 
this  frequency  band,  the  predominant  propagation  problems  are  caused  by  rain:  signal  attenua¬ 
tion,  enhanced  scatter  interference  from  ground  sources,  and  increased  receiver  noise.  The 
current  X-band  propagation  study  is  directed  primarily  toward  the  analysis  of  these  effects. 

1.  Scattering  by  Rain  and  Clouds 

The  microwave  rain  scattering  phenomenon  has  been  known  for  many  years.  Since  World 
War  II,  this  phenomenon  has  been  used  by  radio  meteorologists  for  storm  detection  and  study. 
The  increased  use  of  microwave,  line-of-sight  paths  for  reliable  communications  has  revived 
a  study  of  the  details  of  this  phenomenon.  The  communication  systems  engineer  now  needs  sta¬ 
tistics  of  the  rain  effects:  attenuation,  bistatic  scattering  cross  sections,  and  increased  sky 
temperature.  Each  of  these  scattering  effects  may  be  computed  by  using  simplifying  assump¬ 
tions  about  the  rain  drops,  their  size,  shape,  dielectric  properties,  and  numbers,  about  the 
scattering  process,  by  neglecting  cooperative  effects  between  individual  drops,  and  by  invoking 
single  scattering  theory  to  describe  the  overall  effects  of  an  ensemble  of  drops.  With  these 
assumptions,  computations  of  the  attenuation  and  backscattering  cross  sections  have  been  made 
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2  3  4  5 

for  model,  homogeneous  rain-drop  distributions  (Ryde  and  Hyde  ;  Gunn  and  East  ;  Oguchi  '  ; 

Medhurst^).  Recent  computations  of  bistatic  scattering  cross  sections  have  been  made  by 

7  8 

Deirmendjian  and  of  sky  temperature  change  by  Hogg  and  Semplak. 

The  computations  and,  in  some  cases,  experimental  data  indicate  the  order  of  magnitude 
of  the  scattering  effects.  The  attenuation,  sky  temperature,  and  scattering  statistics  depend 
both  upon  the  per-unit -volume  scattering  statistics  and  on  the  way  these  statistics  are  organized 
in  rain  showers  and  storms.  Weather  radar  data  have  been  used  to  study  the  details  of  the  large 
scale  storm  processes.  These  data  show  that  the  rain  storm  structure  varies  from  storm  to 
storm  and  location  to  location.  Statistical  analysis  of  rain  storm  data  is,  however,  not  avail¬ 
able  for  the  estimation  of  the  rain  scattering  effects  of  interest.  The  only  crude  statistical 

9  10 

estimates  available  of  rain  attenuation  (Bussey  ;  Ugal  and  Kaneda  )  were  made  from  meas¬ 
urements  of  rain  rate  taken  at  a  point  on  the  surface. 

The  program  described  here  for  studying  rain-scattering  effects  consists  of  making  a  de¬ 
tailed  analysis  of  both  the  microscale  and  macroscale  scattering  effects  of  rain.  The  micro¬ 
scale  or  per-unit- volume  scattering  statistics  are  being  studied  by  using  a  large  sample  of 
measured  rain-drop  distributions  and  by  experimentally  checking  the  results  of  scatter  com¬ 
putations  on  artificial  rain-drop  distributions.  The  macroscale  scattering  analysis  consists  of 
a  study  of  the  details  of  storm  organization  through  comparisons  of  measured  sky  temperature 
with  computations  based  on  weather  radar  measurements.  The  goal  of  these  studies  is  a  method 
by  which  climatological  rain  data,  such  as  rain  gauge  statistics  and  weather  radar  statistics 
can  be  converted  into  attenuation,  sky  temperature,  and  interference  statistics. 

2.  Microscale  Studies 

The  interpretation  of  weather  radar  data  in  terms  of  rain  or  scattering  parameters  requires 

a  detailed  knowledge  of  the  integrated  scattering  and  rain  parameters  for  many  different  drop- 

size  distributions.  To  provide  a  statistical  transformation  from  weather  radar  measurements 

of  reflectivity  to  rain  rate,  liquid  water  content,  attenuation  coefficient  and  sky  temperature, 

computations  were  made  for  a  large  sample  of  drop-size  measurements  by  means  of  the  filter 

paper  method  at  M.I.T.,  in  Cambridge*  The  rain  drops  were  assumed  to  be  spherical,  to  be 

all  at  a  constant  drop  temperature,  and  to  have  a  parent  drop-size  distribution  given  by  the 

measured  drop-size  distribution.  The  scattering  effects  were  computed  using  Mie  Theory  as 

described  by  Van  de  Hulst.11  For  these  computations,  the  dielectric  constant  for  water  was 
1 2 

taken  from  Kerr  (Table  8.8,  page  675).  The  rain  rate  was  computed  using  the  terminal  ve- 

13 

locity  data  of  Gunn  and  Kinzer. 

The  effective  reflectivity  factor  (Ze^)  as  defined  by  Eq.  (1)  is  used  to  characterize  the  back- 
scattering  of  an  ensemble  of  scatterers. 

„  6/3  r  4 \ 

Zeff  =  ~ FTTZ  mm  'm  (1) 

7T  |k| 

P  =  backscatter  cross  section  per  unit  volume  in  millimeters 
squared  per  cubic  meter. 


*The  sample  consists  of  338  drop-size  measurements  made  over  a  two-year  period  by  the  M.I.T.  Weather  Radar  Research  Facility 
and  it  represents  many  different  types  of  storms. 
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X  =  wavelength  in  millimeters 


m  =  complex  refractive  index  of  water. 


Figure  IV-1  shows  the  relation  between  rain  rate  and  at  2.8  GHz  obtained  from  Eq.  (1) 
and  the  observed  drop-size  distributions.  This  figure  provides  the  connection  between  values 
of  obtained  from  weather  radar  measurements  and  the  usual  rain  parameter,  rain  rate. 
The  line  labeled 


1.6 


Z  ,,  =  200R 
eff 


where  R  =  rain  rate  in  millimeters  per  hour, 


(2) 


is  normally  used  by  radar  meteorologists  to  provide  this  connection.  The  Laws  and  Parsons 


1 4 

model  represents  the  averages  of  drop-size  measurements  made  over  a  three-year  period  in 


the  Washington,  D.C.  area.  This  model  is  used  for  most  of  the  attenuation  computations  ref¬ 
erenced  above.  The  spread  in  data  points  is  caused  in  part  by  the  natural  variability  of  parent 
drop-size  distributions  that  occur  in  the  many  different  storm  types  represented  by  the  New 
England  storm  sample.  A  second  cause  in  the  spread  in  data  points  is  instrumental,  being 
mostly  due  to  the  use  of  a  limited  sample  volume.  The  relative  importance  of  the  instrumental 
uncertainty  is  presently  unknown  and  is  the  subject  of  continuing  investigation. 

The  attenuation  coefficient  in  decibels  per  kilometer  as  a  function  of  rain  rate  and  of  re¬ 
flectivity  is  presented  in  Figs.IV-2  and  IV-3.  The  spread  in  data  points  about  the  Laws  and 
Parsons  line  is  smallest  at  low  rain  rates  or  at  high  reflectivities.  From  this,  it  is  seen  that 
for  the  very  low  rain  rates,  a  measurement  of  rainfall  rate  not  dependent  upon  weather  radar 
data  is  desirable  while  at  the  highest  rainfall  rates,  weather  radar  data  provides  a  better  meas¬ 
urement  of  attenuation  than  does  a  method  based  on  rainfall  rate  data.  For  this  data  sample, 
less  than  a  factor  of  two  uncertainty  in  the  attenuation  coefficient  is  evident  for  log  Z  values 
greater  than  4.0  or  rain  rates  above  10  to  15  mm/ hr. 

The  sky  temperature  for  a  homogeneous,  constant  temperature  slab  volume  of  rain  drops 
is  given  by 


A  =  total  attenuation  along  path 
0ext  =  extinction  cross  section  per  unit  volume 

Tjj  =  constant  drop  temperature 
d  =  distance  across  slab  volume 
s  =  distance  along  the  path. 
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Fig.  IV— 1  •  Reflectivity  vs  rain  rate  for  New  England  storm  sample. 
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Fig.  IV-2.  8.0-GHz  attenuation  coefficient  vs  rain  rate  for  New  England  storm  sample. 
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Fig.  IV— 3.  8.0-GHz  attenuation  coefficient  vs  reflectivity  for  New  England  storm  sample. 
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g 

This  integral  as  given  by  Hogg  and  Semplak  represents  a  low-order  scattering  solution  to  the 
radiative  transfer  equation.  Sky  temperature  computations  for  an  isolated  slab  at  0.0°  C  are 
given  in  Figs.IV-4  and  IV-5.  These  results  apply  reasonably  well  to  drop  temperatures  that 
may  occur  in  New  England  because,  at  8.0  GHz,  the  attenuation  coefficient  changes  only  slightly 
with  temperature  and  T^  varies  by,  at  most,  ten  percent. 

An  experimental  check  is  being  sought  for  the  above  calculations  by  computing  and  meas¬ 
uring  the  attenuation  and  phase  shift  caused  by  an  artificial  water-drop  distribution  existing  over 
a  distance  of  15  feet  in  the  center  of  a  1000-ft  line-of-sight  propagation  path,  facilities  are 
being  provided  for  measuring  both  the  instantaneous  and  average  values  of  attenuation  and  phase 
shift  for  a  time -stationary  parent  water-drop  distribution.  The  experimental  check  is  neces¬ 
sary  to  answer  some  of  the  questions  raised  by  the  lack  of  experimental  verification  of  scat¬ 
tering  computations  made  for  longer  path  lengths  as  reported  by ’Medhurst.^ 

The  final  part  of  the  microscale  investigation  is  concerned  with  the  bistatic  scattering  cross 
sections  per  unit  volume.  The  bistatic  scattering  cross  section  is  a  function  of  both  the  polar¬ 
ization  of  the  incident  radiation  and  the  scattering  angle.  The  polarization  dependence  can  be 
investigated  by  considering  two  linear,  orthogonal  polarizations,  one  in  the  plane  of  scattering 
and  the  other  perpendicular  to  this  plane.  For  spheres,  using  this  polarization  representation, 
there  are  no  cross -polarization  terms  generated  by  single  scattering  and  the  result  is  a  function 
only  of  the  scattering  angle,  the  angle  between  the  direction  of  propagation  of  the  incident  and 
the  scattered  radiation.  Two  bistatic  scattering  cross  sections  are  sufficient  to  describe  the 
scattering  process,  one  for  polarization  in  the  plane  of  scattering,  /?2(0)  and  the  other  perpen¬ 
dicular  to  this  plane  /L(0).  Figure  IV-6  shows  the  normalized  bistatic  scattering  cross  sections 

1  14 

per  unit  volume  as  a  function  of  scattering  angle  for  the  Laws  and  Parsons  model.  In  this  plot. 


/»!<©> 

I 3 

scat 


P2(0) 


/J2(0) 


scat 


^1’  ^2  =  bistatic  scattering  cross  section  per  unit  volume 


^scat^1^0’*^0*1  sin0d0 


=  total  scattering  cross  section  per  unit  volume. 


(4) 

(5) 


(6) 


The  lines  labeled  Rayleigh  scatterer  depict  the  results  for  Rayleigh  scattering  where  the  phase 
functions  p^(0),  p2(0)  are  given  by 


Pt(0)  =  3/2 

(7) 

p2(9)  =  3/2  cos2  ©  . 

(8) 

A  measure  of  the  departure  of  the  bistatic  scattering  cross  sections  for  rain  at  8  GHz  from 
those  for  Rayleigh  scatterers  is  given  by  the  ratio  of  p^O0)  to  p1(180°)  or  of  the  forward-to- 
backscattered  energy.  Figure  IV-7  gives  this  ratio  as  a  function  of  reflectivity.  A  second 
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Fig.  IV-4.  Antenna  temperature  at  8.0  GHz  vs  rain  rate  for  a  homogeneous, 
1 .0-km-thick  rain  slab. 
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Fig.  IV-5.  Antenna  temperature  at  8.0  GHz  vs  reflectivity  for  a  homogeneous, 
1 .0-km-thick  rain  slab. 
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Fig.  IV-6.  Normalized  bistatic  scattering  cross  sections  at  8.0GHz  vs 
scattering  angle  for  Laws  and  Parsons  model  data. 
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Fig.  IV-7.  Ratio  of  forward-to-backscattering  cross  section  at  8.0  GHz  vs 
reflectivity  for  New  England  storm  sample. 
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Fig.  IV-8.  Comparison  of  reflectivity  calculated  for  8.0  GHz  using  Mie  theory 
with  values  calculated  using  Rayleigh  theory. 
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measure  is  given  by  the  departure  of  the  Z  effective  from  Z  as  given  in  Fig.IV-8.  The  quantity 
Z  is  found  by  substituting  in  Eq.  (1)  the  backscatter  cross  section  calculated  for  an  ensemble  of 
rain  drops  using  Rayleigh  theory.  The  data  show  that,  for  scattering,  an  ensemble  of  rain  drops 
behaves  as  an  ensemble  of  Rayleigh  scatterers,  except  at  the  highest  value  of  reflectivity. 
Rayleigh  scattering  theory,  however,  cannot  be  used  to  compute  attenuation  or  sky  temperature 
at  8  GHz. 

3.  Macroscale  Studies 

The  microscale  studies  described  above  are  oriented  toward  providing  a  means  to  interpret 
weather  radar  data  with  the  least  possible  uncertainty  due  to  rain  storm  types  and  the  location 
of  the  scattering  volume  in  the  storm.  The  macroscale  study  is  directed  at  using  weather  radar 
measurements  of  storm  structure  to  compute  attenuation,  sky  temperature,  and  composite  bi¬ 
static  scattering  cross-section  statistics.  The  method  used  to  establish  the  validity  of  these 
computations  is  to  make  simultaneous  measurements  of  Z e ^  and  sky  temperature,  using  the 
10-cm  weather  radar  operated  by  the  M.I.T.  Weather  Radar  Research  Facility  at  Cambridge 
and  the  X-band  radiometer  facility  at  Haystack. 

The  mode  of  operation  is  to  observe  rain  storms  at  a  distance  from  the  Haystack  site.  Az¬ 
imuth  scans  are  taken  with  Haystack  through  a  rain  storm  at  a  fixed  elevation  angle.  Under 
these  conditions,  the  change  in  sky  temperature  is  due  to  the  emission  from  rain  and  clouds  and 
to  the  changes  in  the  energy  emitted  by  the  surface  of  the  earth  and  collected  by  the  antenna  side 
lobes.  The  contribution  of  the  ground  and  clouds  represents  uncertainties  in  the  measurements 
because  both  the  side -lobe  ground  structure  and  the  cloud  structure  are  not  readily  measurable. 
For  operation  at  elevation  angles  between  0.5  and  2.0  degrees,  measurements  of  the  ground  con¬ 
tribution  show  less  than  a  2.0°  K  fluctuation  after  correction  using  previously  measured  data  for 
each  azimuth.  The  uncertainties  in  measurement  due  to  clouds  is  related  to  the  uncertainty  in 
the  integrated  liquid  water  content  of  the  clouds  along  the  antenna  beam.  For  clouds  at  0.0°C, 
the  sky  temperature  for  an  isolated  1.0-km-thick  cloud  is  given  by 

T  =  4. 0 W  degrees  Kelvin  (9) 

where 

W  =  liquid  water  content  in  grams  per  cubic  meter. 

The  cloud  error  can  be  large,  being  as  much  as  20.0°  K  for  a  large,  nonraining  cumulus  cloud. 
These  errors  far  exceed  the  instrument  errors  caused  by  calibration  uncertainties  and  the  short, 
0.3-second  integration  time. 

The  results  of  a  measurement  made  on  10  August  1965  are  given  in  Figs.IV-9  through  IV-11. 
Figure  IV-9  is  a  weather  radar  map  which  shows  a  section  of  squall  line  visible  both  to  the  M.I.T. 
SCR-615B,  2.8  GHz  weather  radar  and  Haystack.  The  showers  in  the  wedge  bounded  by  60°, 

90°  and  1.0°  elevation  horizon  for  Haystack  are  numbered  1  to  8.  The  horizon  was  chosen  as 
the  distance  at  which  the  1°  elevation  ray  would  reach  5.0  km  height.  Per-unit  volume  reflec¬ 
tivities  computed  from  weather  radar  cross-section  data  are  shown  in  Fig.IV-10  as  reflectivity 
contours.  The  figure  also  includes  a  reflectivity  contour  table  which  gives  the  maximum  and 
minimum  values  of  reflectivity  which  the  volumes  between  successive  levels,  as  described  by 
the  outer  contours,  can  have. 
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Fig.  IV-9. 


Weather  radar  map  for  section  of  squall  line  east  of  the  Haystack  site. 
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.  IV-10.  Reflectivity  cross  section  through  section  of  squall  line  east  of  Haystack  site. 


60  70  80  90 
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Fig.  IV—  1 1 .  Comparison  of  computed  and  measured  antenna  temperature  for  azimuth  scan 
through  section  of  squall  line  east  of  Haystack  site. 
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The  composite  sky  temperature  for  the  eight  showers  is  shown  in  Fig.IV-11.  Showers  1 
and  4  contribute  less  than  15°K  each  and  are  lost  in  the  contributions  of  the  stronger  rain  cells. 
Reflectivity  values  in  excess  of  L5  are  necessary  to  generate  the  peaks  in  the  measured  and 
computed  temperature  data.  The  peaks  due  to  showers  2,  5  and  6  are  obvious.  The  measured 
sky  temperature  peaks  fall  within  the  expected  maximum  and  minimum  temperatures.  The  max¬ 
imum  and  minimum  expected  temperatures  are  generated  by  using  the  Laws  and  Parsons  curve 
in  Fig.IV-3.  The  attenuation  coefficient  for  the  maximum  and  minimum  Ze^  value  for  each 
level  was  multiplied  by  the  distance  along  the  ray  in  the  volume  corresponding  to  the  level  on 
the  weather  map,  then  converted  to  sky  temperature.  The  azimuths  of  the  measured  and  com¬ 
puted  temperature  peaks  2,  5  and  6  do  not  correspond  exactly.  A  range  error  of  the  order  of 
10km  is  possible  at  the  distances  greater  than  100km  from  M.I.T.  This  can  contribute  to  a 
2°  azimuth  error  relative  to  Haystack,  sufficient  to  reconcile  the  computations  and  measure¬ 
ments  for  showers  2,  5  and  6. 

The  results  for  azimuths  between  80°  and  90°  show  a  difference  between  the  measured  tem¬ 
perature,  100°  peak,  and  the  computed  data  113°- 160°,  narrow  peak.  The  source  of  attenua¬ 
tion  at  these  azimuths  is  the  strong  cell  in  shower  3  along  the  cross  section  given  in  Fig.IV-10. 
The  ray  from  Haystack  crosses  the  shower  at  about  the  same  height  as  the  center  of  the  beam 
for  the  weather  radar.  The  map  data  should  be  directly  applicable  to  sky  temperature  compu¬ 
tations  for  this  shower  position.  The  height  data  show  that  the  southern  edge  of  the  cell  is  in 
reasonable  agreement  with  the  map,  but  the  northern  edge  is  not.  An  error  of  as  much  as  4° 
may  be  evident  in  positioning  the  northern  edge  on  the  radar  map.  This  means  that  except  for 
the  sharp  peak  in  the  computed  curve,  the  shapes  of  the  observed  and  computed  curves  can  be 
reconciled.  The  sharp  peak  may  be  due  to  hail  which  can  cause  a  large  backscatter  cross  sec¬ 
tion  but  does  not  contribute  to  the  sky  temperature. 

Peaks  7  and  8  evident  in  the  measured  data  between  azimuths  75°  and  80°  appear  to  be 
caused  by  showers  7  and  8.  The  identification  of  these  peaks  is  difficult  and  the  lack  of  agree¬ 
ment  in  sky  temperature  is  obvious.  Lack  of  agreement  in  the  temperature  for  shower  3  is 
also  evident.  Interpretation  of  results  for  shower  3  is  difficult  because  it  is  in  the  near  field 
of  the  Haystack  antenna  system. 

4.  Summary 

The  rain  propagation  studies  are  progressing  toward  refining  the  estimates  of  attenuation 
and  sky  temperature  by  means  of  weather  radar.  At  present  a  sizable  statistical  spread  in  the 
per -unit -volume  scattering  statistics  and,  in  some  instances,  a  lack  of  agreement  in  the  sky 
temperature  and  weather  radar  measurement  comparisons  are  evident.  These  problems  are 
being  worked  on  by  improving  both  the  computation  procedures  and  the  measurement  techniques. 

R.  K.  Crane 


C.  AURORAL  OBSERVATIONS 
1.  Introduction 

Radar  investigation  of  the  aurora  borealis  using  the  Millstone  L-band  system  and  the  84 -ft 
parabolic  antenna,  which  was  begun  in  February  1965,  was  discontinued  in  December  1965 
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largely  because  radar  echoes  from  the  aurora  were  being  obtained  less  frequently  than  earlier 
in  the  year  and  were  usually  quite  weak.  Analysis  of  the  data  recorded  during  this  period  of 
observation  is  continuing  and  will  be  the  subject  of  a  future  report. 

On  the  majority  of  days  when  auroral  echoes  were  obtained,  they  were  of  relatively  short 
duration  —  from  a  few  minutes  to  almost  an  hour  —  and  were  received  only  at  the  lowest  eleva¬ 
tion  angles,  i.e.,  from  zero  to  about  3°.  On  a  few  days,  however,  the  echoes  persisted  for  pe¬ 
riods  up  to  several  hours  and  were  also  obtained  at  appreciably  higher  elevation  angles  -  on  one 
day  up  to  13°.  The  data  obtained  on  these  days  are  of  particular  interest  because  they  can  be 
used  to  determine  the  spatial  distribution  of  the  scattering  centers  which  are  the  source  of  the 
auroral  echoes.  Some  of  the  data  recorded  on  5  November  1965  have  been  chosen  for  analysis 
in  this  report,  primarily  because  the  echoes  on  this  day  were  the  strongest  obtained  in  nearly 
a  year  of  auroral  observation.  These  data  are  examined  to  determine  the  spatial  distribution 
of  the  auroral  echoes,  their  aspect  sensitivity  and  height  dependence  and  the  radar  cross  sec¬ 
tion  of  these  echoes  for  the  Millstone  L-band  radar  system. 

2.  Observation  Procedures 

The  method  of  recording  auroral  data  consisted  in  stepping  the  antenna  in  5°  increments  of 
azimuth  at  constant  elevation  angle  across  the  sector  from  which  the  auroral  echoes  were  re¬ 
ceived.  At  the  completion  of  each  pass  across  the  sector,  the  elevation  angle  was  changed  by 
0.5°  and  the  antenna  again  stepped  across  the  sector.  At  each  step  the  received  signal  was 
integrated  for  a  period  of  10  seconds.  During  the  time  when  the  data  examined  here  were  ob¬ 
tained,  the  region  producing  the  auroral  echoes  extended  from  about  320°  true  bearing  in  the 
west  to  beyond  15°  in  the  east  and  from  zero-degree  elevation  to  4.5°.  However,  the  region  to 
the  east  of  15°  was  not  investigated  because  of  the  possibility  that  the  Millstone  radar  might 
interfere  with  the  Haystack  site  at  a  bearing  of  about  19°.  The  total  time  required  to  probe  the 
volume  producing  the  auroral  echoes  was  approximately  thirty  minutes. 

3.  Results 

The  shape  of  a  typical  strong  auroral  echo  is  shown  in  Fig.IV-12.  Each  point  in  the  curve 
gives  the  cross-section  per  unit  volume  at  the  indicated  slant  ranges.  These  values  of  cross 
section  are  computed  for  a  volume  whose  diameter  is  determined  by  the  half-power  beamwidth 
of  the  antenna  (0.7°)  and  whose  length  is  determined  by  the  length  of  the  transmitted  pulse 
(500psec).  The  shape  of  this  curve  shows  only  approximately  the  distribution  of  the  echo- 
producing  sources.  It  is  distorted  by  the  width  of  the  pulse  over  the  entire  curve  and  at  the 
leading  and  trailing  edges  by  the  fact  that  at  these  low  elevation  angles  the  beamwidth  of  the  an¬ 
tenna  is  not  filled  over  the  first  and  last  60  to  75  km  of  slant  range.  Height  above  the  earth  for 
the  center  of  the  pulse  at  the  center  of  the  antenna  beam  has  been  computed  from  the  slant  range 
and  the  refraction-corrected  elevation  angle. 

Values  of  cross-section  per  unit  volume,  similar  to  those  in  Fig.IV-12,  were  obtained  for 
each  azimuth  step  of  the  antenna  at  each  elevation  angle.  The  values  at  all  azimuths  for  a  given 
elevation  angle  were  arranged  in  the  form  of  a  grid,  and  contours  which  are  1,  3,  and  6db 
down  from  the  maximum  cross  section  at  this  elevation  angle  were  drawn.  These  contours  for 
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Fig.  IV  —  1 2 .  Distribution  of  auroral  echo  strength  with  slant  ranges,  0.5°  elevation, 

0°  azimuth,  5  November  1965,  17:17:15  -  17:17:25  EST. 

elevation  angles  of  0.5°,  1.5°,  2.5°,  3.5°,  and  4.5°  are  shown  in  Figs.IV-13  through  IV-17, 
respectively.  The  heavy  dot  in  each  figure  gives  the  location  of  the  strongest  portion  of  the  au¬ 
roral  echo  at  the  given  elevation  angle.  The  symmetrical  curves  on  these  figures  show  the  lo¬ 
cation  of  equal  angles  between  the  direction  of  the  radar  ray  and  the  normal  to  the  earth's  mag¬ 
netic  field  lines.  At  zero  angles  the  rays  are  parallel  to  the  normal  and  orthogonal  to  the  field 
lines.  These  curves  are  based  on  computations  made  by  using  a  48-term  Legendre  polynomial 
to  describe  the  earth's  magnetic  field.  (The  computations  were  made  by  I)r.  Shou-Chin  Wang 
of  Lincoln  Laboratory.)  The  height  above  the  earth  was  calculated  from  the  slant  range  as  ex¬ 
plained  previously  for  Fig.IV-12. 

Examination  of  Figs.IV-13  through  IV-17  shows  that  the  strongest  portion  of  the  echo  at 

each  elevation  angle  is  located  very  close  to  the  orthogonality  curve  and  that  the  1-db  contour 

shows  a  very  strong  tendency  to  be  centered  about  this  line.  It  is  evident,  also,  that  most  of 

the  echoes  are  within  about  ±2°  of  orthogonality.  This  is  in  good  agreement  with  Leadabrand, 

15 

et  aL,  and  other  workers  at  lower  frequencies  who  have  found  that  auroral  echoes  are  strongly 
aspect-sensitive  and  that  this  sensitivity  increases  with  increase  in  the  transmitted  frequency. 

The  location  of  the  6-db  contours  gives  some  idea  of  the  heights  from  which  the  auroral 
echoes  are  received.  However,  the  heights  shown  here  are  only  approximate  in  that  they  are  for 
the  center  of  the  antenna  pattern,  whereas  in  reality  the  lower  limit  of  the  auroral  region  height 
is  measured  by  the  upper  edge  of  the  half-power  width  of  the  antenna  beam  and  the  upper  limit 
by  the  lower  edge  of  the  antenna  beam.  From  the  data  considered  in  this  report,  it  was  deter¬ 
mined  that  the  average  thickness  of  the  region  producing  the  auroral  echoes  was  30  km  and  that 
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Fig.  IV  — 1 3.  Spatial  distribution  of  auroral  echoes,  0.5°  elevation, 
5  November  1965,  17:16:30  -  17:19:25  EST. 


Fig.  IV— 1 4.  Spatial  distribution  of  auroral  echoes,  1.5°  elevation, 
5  November  1965,  17:23:30  -  17:26:25  EST. 
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Fig.  IV —15.  Spatial  distribution  of  auroral  echoes,  2.5°  elevation, 
5  November  1965,  1 7:30:30  -  17:33:25  EST. 


Fig.  IV-16.  Spatial  distribution  of  auroral  echoes,  3.5°  elevation, 
5  November  1965,  17:37:30  -  17:40:25  EST. 
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Fig.  IV — 1 7.  Spatial  distribution  of  auroral  echoes,  4.5°  elevation, 
5  November  1965,  17:44:30  -  17:47:25  EST. 


it  extended  from  about  90  to  120km  above  the  earth.  This  is  a  somewhat  greater  thickness  than 

16  17  15 

that  previously  reported  by  Unwin,  Stone,  et  al.,  and  Leadabrand,  et  al.,  but  it  should  be 
emphasized  that  this  value  is  for  a  strong  aurora  whereas  those  reported  by  other  workers  are 
the  average  of  many  auroras,  both  weak  and  strong.  Moreover,  our  data  do  show  that  the  larger 
values  of  layer  thickness  are  associated  with  the  stronger  echoes. 

In  the  last  report  it  was  erroneously  stated  that  the  radar  cross  section  of  L-band  auroras 
could  have  values  up  to  many  tens  of  square  meters  for  the  Millstone  radar.  It  should  have  been 

given  as  tenths  of  square  meters.  The  largest  value  of  radar  cross  section  measured  to  date 

2 

was  approximately  0.4m  . 

Future  work  on  L-band  radio  auroras  at  Millstone  is  indefinite  at  this  time.  However,  it  is 
hoped  that  since  we  are  now  in  the  period  of  increasing  solar  activity  and  the  accompanying  in¬ 
crease  in  the  rate  of  occurrence  of  aurora,  an  opportunity  will  be  found  to  resume  auroral  ob¬ 
servations  on  the  same  modest  scale  as  during  1965.  In  the  meantime,  the  analysis  of  auroral 
data  is  continuing  and  it  is  expected  that  a  report  on  this  work  will  be  out  shortly. 

W.G.  Abel 
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V.  SPACE  SURVEILLANCE  TECHNIQUES 


A.  MILLSTONE  RADAR 

The  Space  Surveillance  Techniques  (SST)  program  for  the  Millstone  radar  was  described  in 

$ 

considerable  detail  in  the  last  Radio  Physics  and  Astronomy  Report.  This  program  consists 
basically  of  two  parts:  (1)  routine  satellite  tracking  operations,  including  support  to  the  MITRE/ 
Millstone  Interferometer  experiment  and  (2)  a  sensor  improvement  program,  including  the  ac¬ 
quisition  of  an  orbit  upgrading  program  for  the  SDS-9300  computer  and  development  of  a  range 
extension  program  utilizing  both  coherent  and  incoherent  signal  integration  techniques.  Routine 
satellite  tracking  continues  on  about  the  same  level  as  described  in  the  earlier  report.  Signifi¬ 
cant  progress  has  been  made  in  the  sensor  improvement  program. 

An  orbit  upgrading  program  has  been  procured  from  TRW  Systems,  Inc.  This  program 
which  is  based  on  ESPOD  (Electronics  Systems  Precision  Orbit  Determination)  prpgrams  de¬ 
veloped  for  the  Air  Force  is  a  two-part  program.  The  first  part  is  a  non-real-time  program, 
called  NRTPOD,  which  makes  use  of  tracking  data  from  one  or  more  sensors  to  update  orbital 
elements.  The  second  part,  called  MHESPOD,  is  a  real-time  program  which  makes  use  of 
nominal  orbital  elements  corrected  by  data  that  are  currently  obtained  by  the  radar  to  upgrade 
pointing  instructions  for  the  remaining  portions  of  a  satellite  pass.  NRTPOD  and  MHESPOD 
have  been  delivered  and  are  presently  being  integrated  into  the  Millstone  Hill  system. 

Discussion  of  the  range  extension  program  in  the  previous  report  described  two  methods  of 
digitizing  the  radar  data  and  passing  it  to  the  computer  at  the  required  rate.  The  first,  called 
DIGIMON,  used  a  special-purpose  computer  and  involved  the  use  of  frequency  multiplexing  and 
a  pair  of  high-speed  analog -to-digital  converters.  The  second,  called  MACS,  uses  sample-and- 
hold  circuits  and  analog-to-digital  converters  to  encode  the  outputs  of  sine  and  cosine  phase  de¬ 
tectors.  DIGIMON  has  since  been  abandoned  because  of  the  inability  to  obtain  A/D  converters 
with  the  required  speed  and  accuracy.  An  expanded  version  of  MACS  has  been  designed  which 
normally  uses  six  10-bit  A/D  converters  to  encode  the  sine  and  cosine  outputs  of  one  sum  re¬ 
ceiver  and  the  two  difference  receivers  in  the  monopulse  system.  An  additional  pair  of  A/D 
converters  is  provided  which  can  be  used  when  needed  to  encode  the  output  of  a  fourth  receiver. 
The  new  MACS  assembles  a  20-bit  word  consisting  of  10  bits  each  from  two  data  channels  and 
passes  it  to  the  computer  in  approximately  the  1.75-psec  cycle  time  of  the  memory.  The  A/D 
converter  system  which  employs  integrated  circuitry  has  been  ordered  for  delivery  in  July  1966. 

The  computation  rates  required  in  the  range  extension  program  are  very  high.  The  computa¬ 
tion  problem  was  discussed  in  the  last  report  and  one  method  of  attacking  it  by  Fourier  analysis 
using  the  Cooley-Tukey  algorithm  was  mentioned.  A  seemingly  more  promising  approach  is 
being  investigated  which  involves  a  correlation  filter  (CFE-1)  made  by  Scientific  Data  Systems. 
The  CFE-1  is  a  special-purpose  high-speed  accessory  designed  for  use  with  the  SDS-9300  com¬ 
puter.  It  has  a  basic  word  length  of  12  bits  and  it  can  perform  10  complete  multiply-and-add 


*  Radio  Physics  and  Astronomy  Report,  Lincoln  Laboratory,  M.I.T.  (1965: 1),  DDC  629923. 
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routines  in  one  memory  cycle  (1.75  fisec)  of  the  SDS-9300.  It  can  obtain  an  autocorrelation  func¬ 
tion  directly,  and,  in  conjunction  with  appropriate  programs  in  the  9300,  can,  at  a  great  saving 
in  time,  perform  most  of  the  digital  filtering  that  is  usually  done  by  the  main  computer  alone. 
These  properties  of  the  CFE-1  make  it  a  useful  tool  for  many  other  experiments  requiring  very 

high-speed  signal  processing.  j.  C.  Henry  D.  P.  Hynek 

G.  M.  Hyde  V.  C.  Pineo 


B.  HAYSTACK  X-BAND  PASSIVE  AUTOTRACKER 

An  X-band  passive  tracking  capability  utilizing  the  120-ft  antenna  has  been  developed  at 
Haystack  Hill.  This  tracker  is  intended  to  angle  track  cooperative  satellites  such  as  Lincoln 
Experimental  Communications  Satellites  LES-2  and  LES-4.  These  satellites  are  in  eccentric 
orbits  with  great  apogee  heights;  consequently,  conventional  satellite  radar  trackers  such  as 
Millstone  can  observe  only  a  small  portion  of  the  orbit  near  perigee.  However,  the  Haystack 
large-aperture  precision  angle  tracker  can  obtain  good  time -angle  data  on  approximately  an 
eight -tenths  segment  of  a  LES-4  10-hour  orbit  in  20  percent  of  the  passes.  Data  obtained  on 
such  a  large  portion  of  the  orbit  enable  a  precise  determination  of  the  orbital  parameters.  The 
accuracy  of  this  method  of  orbit  determination  has  been  demonstrated  on  LES-4.  Haystack 
angle  data  on  this  satellite  taken  on  24  and  25  February  1966  used  in  an  orbit  fitting  program  had 
an  RMS  error  of  only  0.0061  deg. 

The  Haystack  passive  tracker  employs  solid  state  electronics  with  a  four-horn  monopulse 
type  feed  in  an  auto-track  system  in  which  the  tracking  loop  is  closed  through  the  Univac  U-490 
computer.  A  duplexed  single-channel  receiver  is  used  for  the  elevation  and  azimuth  error  sig¬ 
nals.  Error  voltages,  determined  with  phase  detectors,  are  digitized  and  fed  to  the  D-490.  An¬ 
tenna  commands  are  given  as  digital  numbers  to  the  digital  angle  servos.  Doppler  frequency 
tracking  is  accomplished  with  a  computer-directed  local  oscillator  and  a  phase-locked  loop. 

R.  R.  Silva 


C.  RUBY  LASER  RADAR 

A  program  for  studying  the  extension  of  radar  techniques  into  the  region  of  infrared  and  op¬ 
tical  wavelengths  is  in  progress.  High  peak  powers  are  available  from  ruby  lasers  operated  in 
the  MQ-switchedM  mode  where  single  pulses  of  250  megawatts  or  more,  with  about  20  nanosecond 
duration,  may  be  produced  at  up  to  one  pulse  per  second  rates.  With  such  high  peak  powers 
available  in  narrow  beams  (about  one  milliradian),  it  has  been  possible  to  obtain  strong  reflections 
from  orbiting  satellites  which  have  been  equipped  with  cube-corner  optical  reflectors.  Such  short 
pulse  durations  enable  precise  timing  of  the  return  echoes  such  that  accurate  satellite  ranges  may 
be  obtained.^  The  low  repetition  rate  and  the  narrow  beam  angle  of  the  laser  effectively  prevent 
acquisition  and  tracking  with  the  laser  radar  alone.  However,  by  interconnecting  the  laser  mount 
with  a  microwave  radar  capable  of  tracking  satellites  within  the  accuracy  required  by  the  laser 
beam,  laser  echoes  are  not  necessary  for  tracking  but  may  be  used  for  precise  ranging.  In 


*11.  M.  Plotkin,  Quantum  Electronics,  Proceedings  of  the  Third  International  Congress  (Columbia  University  Press,  New  York,  1964). 
t  P.  H.  Anderson,  C.  G.  Lehr,  L.  A.  Maestre  and  G.  L.  Snyder,  IEEE  J.  Quantum  Electron.  QE-2  (August  1966). 
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addition,  microwave  radar  tracking  allows  laser  ranging  during  any  period  satellites  are  over¬ 
head,  not  solely  when  they  are  sunlit  against  a  darkened  sky. 

The  Millstone  L-band  radar  is  able  to  track  small  satellites  to  within  one  milliradian,  hence 
can  provide  pointing  information  to  a  laser  mount.  A  Nike-Ajax  radar  mount  has  been  modified 
to  accept  a  ruby  laser  head  and  a  ten-inch  receiving  telescope.  Position  synchros  on  the  mount 
are  interconnected  to  position  synchros  on  the  Millstone  antenna  to  allow  slave  tracking  by  the 
laser  mount.  The  mount  has  been  installed  on  a  40-ft  tower,  Tower  3,  located  approximately 
3  50  ft  west  of  the  Millstone  84-ft  antenna.  A  photograph  of  the  mount  with  the  laser  radar  being 
pointed  by  Millstone  is  shown  in  Fig.  V-l.  Operating  controls  are  provided  in  a  van  located  be¬ 
side  the  base  of  the  tower,  with  supporting  laboratory  space  available  in  a  50-ft  trailer.  The 
pointing  mount  is  installed  on  a  hexagonal  wooden  platform  with  side  walls  and  enclosed  by  a 
folding  aluminum  frame  covered  with  weatherproofed  canvas.  In  order  to  allow  for  optical  align¬ 
ment  of  the  equipment  during  daylight  hours,  the  canvas  cover  has  a  black  interior  surface  and 
has  a  white  exterior  to  minimize  the  temperature  rise  inside. 


Fig.  V-l.  Ruby  laser  radar  (right)  and  Millstone  84-ft  tracking  antenna  (left). 


The  Millstone  antenna  beam-pointing  direction  is  sometimes  as  much  as  a  milliradian  from 

* 

the  synchro  indications  over  its  full  elevation  range,  because  of  mechanical  deflections.  There¬ 
fore,  the  interconnecting  synchro  line  has  been  provided  with  a  pair  of  position  correction  syn¬ 
chros  that  allow  both  elevation  and  azimuth  position  signals  to  be  modified  by  handwheel  controls. 
A  boresight  television  camera  allows  the  mount  operator  in  the  control  van  to  view  an  illuminated 
object  (such  as  the  Echo  balloons)  during  a  track.  By  observing  a  number  of  Echo  I  and  II  tracks 
on  the  boresight  monitor,  it  was  possible  to  determine  a  set  of  elevation  corrections  to  apply  to 


*  Data  from  the  Millstone  computer  are  corrected  for  this  sort  of  error  as  well  as  refraction.  The  optical  radar,  however,  since  it  is 
synchro-positioned  from  the  antenna  directly,  receives  uncorrected  angle  data. 
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Fig.  V-2.  Ruby  laser  radar  on  Nike-Ajax  mount. 


RANGE  vs  TIME 

Fig.  V-3.  Range  determination  for  Explorer  XXII. 
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the  Millstone  antenna  positions.  Using  these  corrections,  it  is  possible  to  hold  the  image  of  the 
Echo  balloons  well  within  one  milliradian  of  the  Nike  mount  axis. 

The  ruby  laser  employed  in  this  program  is  a  commercial  unit  (supplied  by  Applied  Lasers, 
Inc.,  Stoneham,  Mass.)  consisting  of  a  5/8  in.  diameter  by  7-1/4  in.  long  ruby  mounted  for  out¬ 
put  at  Brewster’s  angle  and  excited  by  two  linear  flashlamps  in  a  dual  elliptical  cavity.  The  op¬ 
tical  cavity  is  Q-switched  by  a  combination  of  a  rotating  quartz  rooftop  prism  and  a  saturable 
dye  cell.  Circulating  water  (in  winter  months,  and  alcohol-water  mixture)  maintains  the  temper¬ 
ature  of  the  lamps  and  ruby  at  a  constant  temperature.  The  operating  temperature  of  the  laser 
is  selected  so  that  the  output  wavelength  does  not  coincide  with  any  of  the  narrow  atmospheric 
absorption  lines.  The  20-kw  power  supply  and  the  water  circulation  system  are  located  inside 
the  tower  below  the  mount. 

Next  to  the  laser  head  is  a  10-inch  Cassegrainian  telescope  for  reception  of  the  satellite 

o 

echoes.  A  red-sensitive  photomultiplier  with  a  30-A  interference  filter  provides  skylight  dis¬ 
crimination,  and  the  pulse  output  is  fed  through  an  adjustable  height  discriminator.  This  passes 
pulses  which  exceed  a  predetermined  threshold  level  during  a  gate  period.  The  discriminator 
output  pulse  is  used  to  stop  a  time  interval  counter  which  has  been  initiated  by  the  outgoing 
laser  pulse.  Time  intervals  are  resolved  to  a  precision  about  one  nanosecond  with  the  unit  in 
use.  The  laser  transmitter  and  receiving  telescope  in  place  on  the  mount  are  shown  in  Fig.  V-2. 

Echoes  from  the  satellite  Explorer  XXII  have  been  obtained  for  ranges  between  580  and  1000 
nautical  miles.  Data  given  in  Fig.  V-3  show  the  laser  radar  range  measurements  as  solid  cir¬ 
cles  and  position  prediction  (from  month-old  orbital  parameters)  as  open  circles.  Time-of-day 
determinations  for  these  observations  are  good  only  to  about  3  seconds,  because  recording  was 
performed  manually. 

Modification  of  the  laser  radar  is  currently  in  progress.  The  10-in.  telescope  has  been 

o 

replaced  by  a  23-in.  telescope.  A  thermoelectr ically  cooled  photomultiplier  housing  and  a  5-A 
interference  filter  are  being  included  to  increase  threshold  sensitivity  of  the  receiving  system. 
Accuracy  in  both  the  time  of  observation  and  time  interval  measurements  is  being  improved  by 
referring  a  digital  clock  and  the  interval  unit  to  the  station  frequency  standard  located  at  Haystack. 
Automatic  data  recording  is  also  being  provided.  jj  ^  Bostic  k 

P.  R.  Longaker 
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VI.  SYSTEMS  DEVELOPMENT 


A.  EQUIPMENT  DEVELOPMENT  FOR  FOURTH  TEST  OF  GENERAL  RELATIVITY 

A  new  high-accuracy  high-sensitivity  ranging  radar  is  being  developed  at  Haystack  for  the 
measurements  required  for  the  fourth  test  of  general  relativity.  The  detection  ability  of  the 
radar  will  be  improved  by  about  lOdb  over  that  provided  in  the  system  currently  in  use  for  plan¬ 
etary  observation.  System  parameters  for  a  new  plug-in  box  for  planetary  radar  astronomy  are 
given  below. 

Frequency  7.840  GHz 

Transmitter  power  500-kw  CW 

Antenna  gain 

(including  radome  effect)  66  db  over  isotropic 

System  temperature  75 °K 

The  increase  in  performance  will  be  achieved  by  the  use  of  a  newly  developed  klystron  power 
amplifier  and  a  traveling-wave  maser.  Varian  Associates  has  completed  initial  testing  of  a 
klystron  at  200  kw.  The  objective  of  the  development  program  is  a  250-kw  klystron.  Two  klys¬ 
trons  will  be  operated  in  parallel  to  produce  from  400  to  500  kw.  The  transmitter  power  output 
will  be  limited  to  400  kw  average  until  the  safe  operating  margins  on  the  waveguide  components 
have  been  established.  A  TW  maser  has  been  completed  and  delivered  to  Lincoln  Laboratory. 

The  maser  contributes  10°K  to  the  system  noise  temperature  while  operating  with  gain  in  excess 
of  35  db  and  with  a  bandwidth  of  25  MHz. 

The  system  block  diagram  shown  in  Fig.  VI-I  gives  the  configuration  of  the  microwave  trans¬ 
mitting  and  receiving  equipment  installed  in  the  planetary  radar  box  and  the  frequency  translation. 


Fig.  VI-1.  System  block  diagram  for  fourth  test  experiment. 
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control  and  digital  data  processing  equipment  located  in  the  Haystack  control  room.  The  antenna 
will  be  illuminated  with  a  circularly  polarized  feed  horn.  Initially  a  single  receiver  channel  po¬ 
larized  orthogonally  with  respect  to  the  transmitted  signal  will  be  used  for  receiving.  A  mechan¬ 
ically  operated  switch  is  provided  to  prevent  overloading  the  maser  during  the  transmitting  in¬ 
terval.  Coherence  of  the  transmitted  and  received  signals  is  established  by  a  phase-locked  stable 
LO.  An  appropriately  phase-coded  IF  signal  is  mixed  with  the  LO  output  to  provide  the  excitation 
for  the  transmitter.  The  received  signal  is  heterodyned  to  an  IF  of  130  MHz,  after  being  ampli¬ 
fied  by  the  TW  maser.  This  IF  signal  is  further  processed  by  equipment  located  in  the  control 
room  where  the  Doppler  frequency  shift  is  removed  and  the  signal  data  digitally  encoded  for  com¬ 
puter  processing.  Improved  ground  level  receiving,  data  processing  and  control  instrumentation 
are  required  for  the  fourth  test  experiment.  A  description  of  this  new  instrumentation  follows. 

1.  Receiver  Ground  Equipment 

The  Doppler  shift  of  the  received  signal  is  removed  by  offsetting  the  frequency  of  the  100-MHz 
second  receiver  LO  using  a  coherent  frequency  synthesizer.  A  slewed  1-MHz  clock  signal  is  de¬ 
rived  from  the  Doppler  synthesizer  output  which,  when  appropriately  divided  in  frequency,  is 
used  by  the  radar  timing  equipment  to  track  the  target  in  range.  The  frequency  synthesizer  is 
externally  programmed  to  track  the  predicted  Doppler  frequency.  The  existing  equipment  steps 
the  synthesizer  smoothly  in  0.1 -Hz  steps,  with  provision  for  extending  this  to  0.01 -Hz  increments 
for  greater  precision.  The  frequency  synthesizer  is  controlled  either  by  the  Univac  490  pointing 
computer,  using  a  computed  Doppler  ephemeris,  or  by  a  device  which  is  programmed  to  produce 
a  linear  rate  of  change  in  frequency  with  time  at  a  rate  which  best  matches  the  predicted  Doppler 
frequency  rate. 

2.  Terminal  Equipment 

A  new  set  of  terminal  equipment  operating  at  a  2-MHz  IF  frequency  has  been  designed  for 
the  experiment  and  is  currently  under  construction.  This  terminal  equipment  provides  phase  de¬ 
tection  (zero  frequency  conversion)  of  the  radar  signals  to  generate  the  two  orthogonal  components 
to  be  analyzed  by  computer  processing.  Matched  filtering  of  these  signals  to  the  pulse  width  or 
baud  length  of  the  modulation  is  used  to  obtain  the  highest  signal-to-noise  ratio. 

The  filtered  signals  are  converted  to  digital  form  for  computer  processing  with  an  analog 
sample-and-hold  multiplexer  and  a  high-speed  analog-to-digital  converter.  A  versatile  inter¬ 
face  subsystem  has  been  constructed  and  partially  checked.  This  subsystem  transfers  the 
digitally  encoded  data  from  the  converter  to  the  Control  Data  Corporation  3200  data  processing 
computer.  This  interface  also  reads  time  from  the  station  clock  for  use  in  the  data  processing 
procedure,  and  provides  a  channel  for  the  computer  program  to  control  events  or  devices  ex¬ 
ternal  to  the  computer  or  to  ascertain  the  status  of  the  experiment  (i.e.,  whether  transmitting 
or  receiving). 

3.  Radar  Timing  Sequencer 

The  radar  sequencer  equipment  presently  under  construction  will  control  the  system  timing 
for  the  NASA  Lunar  Mapping  Program  and  for  the  instrumentation  of  Shapiro's  fourth  test  of  the 
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general  theory  of  relativity.  Subsystem  timing  (pulse  length,  repetition  rate,  etc.)  will  be  con¬ 
trolled  by  existing  equipment  under  supervision  of  the  radar  sequencer. 

The  sequencer  unit  will  provide  for  operator  control  of  the  following  functions: 

(a)  System  start 

(b)  Phase  code  modulation  of  transmitted  signal 

(c)  Radar  return  processing 

(d)  System  stop. 


Various  status  signals  will  be  available  to  drive  external  equipment  and  to  present  a  visual  in¬ 
dication  of  the  overall  system  operation.  In  addition,  several  synchronization  tests  will  be  im¬ 
plemented  and  the  results  displayed  on  the  operator  control  panel.  Facilities  will  be  provided 
for  remote  sensing  of  all  status  and  test  signals  in  anticipation  of  future  requirements. 

The  radar  sequencer  is  presently  under  construction  with  a  1  July  1966  target  date  to  start 
checkout.  Unit  checkout  will  require  a  minimum  of  2  to  3  weeks,  at  which  time  the  system  tests 
are  scheduled. 

M.  L.  Stone 
R.  P.  Ingalls 


B.  THE  MILLIMETER-WAVE  FACILITY  AT  LEXINGTON 

1.  Introduction 

Lincoln  Laboratory  procured  a  28 -ft-diameter  antenna  in  1961  as  part  of  a  millimeter-wave 
system  development  program.  The  reflector,  which  was  required  to  operate  at  wavelengths  in 
the  region  of  8  mm,  was  formed  by  a  technique  known  as  "spincasting"  by  the  Kennedy  Antenna 
Division  of  the  Electronic  Specialty  Corporation.  This  technique  for  constructing  precision 
paraboloidal  reflectors  with  a  large  aperture-to-wavelength  ratio  has  been  described  by  Dawson.2 
A  more  complete  description  of  the  antenna  and  its  electrical  performance  is  given  in  Appen¬ 
dix  B  of  this  report.  Evaluation  tests  of  the  28-ft  reflector  were  performed  using  an  El-Az 
mount  installed  on  the  roof  of  Lincoln  Laboratory  Building  D,  and  a  test  transmitter  at  a  dis¬ 
tance  of  about  6  miles.  The  results  of  these  tests  were  so  favorable  that  the  antenna  was  later 
used  as  the  basis  of  a  radar  operating  at  8.6  mm  and  for  millimeter-wave  radiometry. 

The  first  data  on  the  reflectivity  of  the  moon  at  millimeter  wavelengths  were  obtained  with 
this  radar  and  were  reported  by  Lynn,  et  al.3  Use  of  the  radar  was  discontinued  following  the 
moon  measurements  because  of  lack  of  funds  required  to  accomplish  a  significant  improvement 
in  performance.  However,  the  antenna  has  since  been  used  for  radiometry  at  8.5mm  (Ref.  4) 
and  1.5  cm.  An  account  of  the  radiometric  measurements  at  1.5  cm  is  given  in  Sec.  Ill- A  of  this 
report. 

Recently,  as  a  result  of  renewed  interest  in  radar  studies  of  the  lunar  surface,  the  radar 
has  been  reassembled  and  a  program  to  upgrade  its  performance  has  been  undertaken. 

2.  The  8. 5 -mm  Radar 

The  present  radar  is  essentially  the  same  as  that  used  in  the  moon  experiments  reported 
by  Lynn,  et  al.  The  transmitter  consists  basically  of  a  10-watt  floating-drift-tube  klystron 
oscillator,  phase-locked  to  the  required  frequency  which  is  derived  from  a  crystal-controlled 
frequency  multiplier.  The  klystron  is  an  Elliott  Brothers  (England)  type  8FK1. 
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The  front  end  of  the  receiver  is  a  balanced  mixer  using  1N53CMR  crystal  diodes.  The  noise 
figure  is  about  14 db.  Triple  conversion  reduces  the  signal  to  2.5  kHz,  after  which  it  passes 
through  a  bandpass  filter  into  a  detector  and  a  pen  recorder.  Equipment  now  being  prepared  will 
integrate  the  detected  signal  and  register  the  received  energy  on  an  IBM  card.  The  present  band¬ 
width  of  the  receiver  is  170 Hz,  somewhat  more  than  the  greatest  Doppler  spread  caused  by  the 
moon's  libration.  Doppler  shift  introduced  by  rotation  of  the  earth  and  by  ellipticity  of  the  moon's 
orbit  is  compensated  by  programming  the  second  LO. 

3.  Millimeter  Radar  Improvement  Program 

The  first  improvement  will  be  made  in  the  transmitter  with  the  replacement  of  the  10-watt 
8FK1  klystron  by  a  50-watt  8FK15  during  the  latter  half  of  1966.  The  50-watt  klystron  oscillator 
is  similar  to  the  10-watt  tube  presently  in  use  and  phase-locking  to  a  frequency  standard  will  be 
necessary  in  order  to  attain  the  frequency  stability  required  in  moon  experiments.  A  1000-watt 
klystron  amplifier  is  being  developed  by  Varian  Associates  for  Lincoln  Laboratory  in  connection 
with  the  NASA  lunar  studies  program. 

Another  important  improvement  will  be  the  introduction  of  polarization  diversity  to  permit 
three  modes  of  operation: 

(a)  Transmit  circular  polarization  and  receive  two  circular  polarizations  of 
opposite  sense. 

(b)  Transmit  circular  polarization  and  receive  orthogonal  plane  polarizations, 
oriented  at  will. 

(c)  Transmit  plane  polarization  oriented  at  will  and  receive  in  that  plane  and 
in  the  orthogonal  plane. 

It  is  expected  that  at  millimeter  wavelengths  there  will  be  a  very  appreciable  amount  of  depolari¬ 
zation  by  most  or  all  parts  of  the  moon,  and  that  a  useful  signal-to-noise  ratio  will  be  found  in 
both  channels  of  the  dual-polarization  receiver  even  when  the  transmitted  power  is  only  50  watts. 
Later,  when  the  1 -kw  tube  is  delivered,  it  may  be  possible  to  split  the  beam  into  narrow  strips 
by  Doppler  filtering,  and  thus  to  examine  specific  features  of  the  lunar  surface  which  have  ex¬ 
hibited  anomalous  scattering  behavior  at  longer  wavelengths  (Thompson  ’  ). 

4.  Radiometry  with  Millimeter-Wave  Facility 

Heretofore,  the  usefulness  of  the  antenna  in  radiometry  has  been  limited  by  the  coarseness 
of  the  pointing  mechanism;  five  meshes  of  low  precision  gears  have  intervened  between  the  in¬ 
dicating  equipment  and  the  azimuth  or  elevation  axis  of  the  dish.  Sixteen-bit  angular  position 
encoders  are  now  being  placed  directly  on  the  azimuth  and  elevation  axes.  When  binary-to- 
decimal  round-off  errors  are  allowed  for,  these  encoders  will  indicate  position  with  maximum 
errors  of  +23  or  -32  seconds  of  arc.  Elasticity  of  the  mount  and  thermal  displacement  of  the 
supporting  columns  will  affect  the  fiducial  positions,  but  these  factors  can  be  measured  and 
considered.  It  seems  realistic  to  hope  for  pointing  indication  with  error  not  exceeding  one  min¬ 
ute  of  arc,  which  is  a  quarter  of  the  3-db  beamwidth. 

A  second  limitation  on  the  use  of  this  facility  for  radiometry  has  been  the  sensitivity  of 
the  available  receivers.  At  35 GHz,  the  noise  temperature  of  the  radiometer  receiver  has  been 
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more  than  3300°K.  Parametric  amplifiers  available  at  35 GHz  have  not  been  attractive  for  ra- 
diometry,  because  the  few  that  have  been  built  have  had  narrow  pass  bands.  A  new  development 
in  varactors  now  gives  promise  of  a  3  5-GHz  parametric  amplifier  with  a  noise  temperature  of 
only  about  20°K  and  a  bandwidth  of  at  least  1  GHz.  The  effective  temperature  of  the  system,  in¬ 
cluding  waveguide  and  atmosphere,  would  be  about  330°K.  The  possibility  of  procuring  a  pair 
of  these  amplifiers  to  be  used  with  the  dual-polarization  plumbing  is  being  considered. 

The  existing  parameters  and  planned  improvements  for  the  microwave  facility  are  summa- 


rized  below. 

Antenna: 

Half-power  beamwidth  (one 
way)  at  3 5  GHz 

4.0' 

Gain,  including  1.4-db  losses 

65.9  db 

Polarization  (f”®tUlg 

Horizontal 

Complete  flexibility 

Transmitter: 

Existing  —  locked  oscillator 

10  watts 

Near  future  —  locked  oscillator 

50  watts 

Goal  -  coherent  amplifier 

1000  watts 

Receivers: 

Existing  —  single  channel 

3300°K 

Goal  —  dual  channel 

330°  K 

J.J.G.  McCue 
A.  F.  Standing 
E.  A.  Crocker 
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APPENDIX  A 

EXPERIMENTAL  EVALUATION  OF  1000-WAVELENGTH  ANTENNA 

W.D.  Fitzgerald,  V.  L.  Lynn,  K.  J.  Keeping 


A.  INTRODUCTION 

The  recent  demands  of  space  surveillance  and  space  communications  have  greatly  increased 
the  need  for  very-high-gain,  narrow-beam  antennas.  Until  recently,  this  need  has  been  fulfilled 
by  increasing  the  physical  aperture  of  reflectors  operating  at  wavelengths  of  3  cm  or  longer. 
However,  the  higher  frequencies  of  the  millimeter  band  can  improve  gain  and  beamwidth  with 
relatively  small  structures,  while  simultaneously  permitting  extremely  wide  bandwidths,  opera¬ 
tion  in  a  relatively  unused  portion  of  the  spectrum,  small  size  of  waveguide  components  and 
reduced  interference. 

Following  this  reasoning,  Lincoln  Laboratory  initiated  a  program  of  millimeter  systems 
development.  In  conjunction  with  this  effort,  a  28-ft-diameter  antenna,  operating  at  wavelengths 
of  8.5  mm,  has  been  installed  and  tested.  The  reflector  was  formed  by  a  technique  known  as 
"spincasting"  by  the  Kennedy  Antenna  Division  of  the  Electronic  Specialty  Corporation,  with 
engineering  and  technical  assistance  from  Lincoln  Laboratory. 

This  report  describes  the  evaluation  of  the  test  range  and  the  electrical  testing  of  the  an¬ 
tenna  at  a  wavelength  of  8.5  mm. 

B.  TEST  OBJECTIVES 

The  test  series  described  in  this  report  was  initiated  to  determine  the  electrical  charac¬ 
teristics  more  precisely  than  had  been  done  in  the  initial  tests  performed  shortly  after  installa¬ 
tion,  and  to  evaluate  the  effect  of  environmental  conditions  on  these  characteristics.  Gain, 
pattern,  and  boresight  shift  measurements  were  taken  to  determine  the  effects  of  different 
ambient  temperatures,  nonuniform  solar  heating  and  wind  conditions.  On  a  typical  daily  run, 
measurements  were  taken  at  one  to  one-and-one-half  hour  intervals  beginning  before  dawn  and 
continuing  until  the  late  afternoon.  Approximately  15  complete  daily  runs  were  made  during 
various  weather  conditions  which  included  winds  up  to  25-30  mph  and  ambient  temperatures 
ranging  from  the  mid-30's  (°F)  to  the  high  80' s. 

The  properties  of  a  test  range  capable  of  measuring  an  antenna  with  a  diameter  of  1000  wave¬ 
lengths  pose  special  problems/1  Some  effects  of  terrain  and  propagation  phenomena  which  are 
negligible  in  the  measurement  of  antennas  with  smaller  electrical  apertures  become  serious 
obstacles  with  a  1 000-wavelength  aperture.  For  the  range  in  question,  the  directivity  of  the 
28-ft  antenna  at  these  frequencies  is  sufficient  to  virtually  eliminate  the  influence  of  ground 
reflections  from  the  measured  patterns.  However,  these  ground  reflections  can  still  cause 
errors  in  the  gain  measurement  (the  method  used  is  comparison  with  a  calibrated  standard), 
and  propagation  phenomena  can  cause  errors  in  all  the  major  characteristics  of  the  antenna 
(gain,  patterns  and  boresight).  The  propagation  phenomena  of  interest  (refraction  and  signal 
scintillations),  although  only  vaguely  associated  with  site  geometry,  are  considered  as  properties 
of  the  range. 
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Fig.  A-l.  28-ft  millimeter  antenna 
(front  view). 


Fig.  A-2.  28-ft  millimeter  antenna 

(rear  view). 


Fig.  A-3.  Cross  section  of  antenna  structure  —  basic 
geometry  of  panels  and  backup. 
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An  investigation  of  the  range  to  determine  the  dependence  of  terrain  and  propagation  effects 
on  wind,  sun,  temperature,  etc.,  was  undertaken  with  a  view  toward  separating  these  "range 
errors"  from  actual  changes  in  antenna  performance  caused  by  the  same  changing  environment. 

C.  MECHANICAL  DESCRIPTION 

Photographs  of  the  antenna  mounted  on  the  roof  of  one  of  the  Lincoln  Laboratory  buildings 
are  shown  in  Fig.A-1  and  Fig.A-2. 

The  reflector  is  supported  by  an  aluminum  plate-girder-tangential-beam  structure  which  is 
shown  in  Fig.A-2.  Also  shown  in  Fig.A-2  is  the  cab  between  the  counterweights  which  houses 
all  the  microwave  equipment.  This  backup  structure  is  attached  to  the  mount  adapter- stiffener 
at  6  hexagonally  distributed  points.  The  reflecting  surface  is  supported  by  a  total  of  24  panels 
individually  attached  to  the  backup  structure.  Each  panel  is  6  in.  thick  and  constructed  of  fiber¬ 
glass  over  paper-core  honeycomb.  The  precise  spuncast  paraboloid  is  formed  over  these  panels. 
The  actual  reflecting  surface  is  a  thin  flame- sprayed  coating  of  zinc  which  is  then  weather  pro¬ 
tected  by  a  thin  coating  of  plastic.  Figure  A-3  shows  a  cross  section  of  the  reflector  construction. 

The  feed  employed  is  a  "rear  feed"  type  supported  by  an  aluminum  column  attached  to  the 
vertex  and  guyed  at  the  feed  end.  The  feed,  which  is  described  more  fully  later,  consists  of  an 
open-ended  circular  waveguide  operating  in  the  TE^  mode  and  a  circular  splash  plate  housed 
in  a  hemispherical  radome.  The  radome,  feed  support  column  and  guy  wires  are  shown  in 
Fig.A-1.  The  cylindrical  object  seen  attached  to  the  feed  support  is  an  auxiliary  receiver  which 
was  removed  after  initial  testing  of  the  dish. 

To  permit  focusing  of  the  antenna  without  the  use  of  flexible  waveguide  or  movable  joints, 
all  the  microwave  components  of  the  radar  and  pattern  receiver  are  mounted  on  a  single  large 
plate  in  the  equipment  housing.  Adjustment  is  accomplished  by  movement  of  the  entire  plate 
and  feed. 

The  mount  is  a  manually  controlled  elevation-over-azimuth  type.  Both  the  azimuth  and 
elevation  control  points  are  provided  with  telescopes  mounted  on  the  elevation  axis  and  directed 
through  small  apertures  in  the  reflecting  surface.  One  telescope  is  a  4-inch  Cassegrainian  type 
used  for  fine  control.  Rotation  of  the  mount  is  limited  to  345°  in  azimuth  and  from  -2°  to  +95° 
in  elevation.  A  two-speed  synchro  system  provides  information  for  remote  readout  and  for  the 
pattern  recorder. 

A  more  detailed  mechanical  description  and  an  account  of  the  "spincasting"  technique  has 

2  3 

been  given  by  Lynn,  etal.  and  Dawson. 

D.  PATTERN  RANGE 

Figure  A -4  shows  a  vertical  cross  section  of  the  test  range  used  in  the  evaluation  program. 
The  transmitting  antenna  (Billerica)  is  a  4-ft  paraboloid  mounted  on  top  of  a  1 3  0  -  ft  water  tower 
located  approximately  31,500  feet  from  the  antenna  under  test.  The  power  source  is  a  reflex 
klystron  (Varian  VA-97)  with  20-  to  25-mw  CW  output,  remotely  controlled  with  telephone-line 
relay  circuits.  The  polarization  of  the  transmitting  antenna  was  horizontal  throughout  the  test 
program. 

The  intensity  with  which  the  transmitter  illuminates  the  terrain  is  as  indicated  in  Fig.  A-4. 
The  100-  to  1 5  0  -  ft  valley  between  the  two  sites  minimizes  the  interference  pattern  observed  at 
the  receiving  site.  This  interference  pattern  is  of  consequence  only  to  the  relatively  low-gain 
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Fig.  A-4.  Height  profile  of  pattern  range. 


13-30-10293 1 

8-in. D I  AMETER  STANDARD  GA! N  HORN 

PARABOLOID  ( 2-7/t6  in  x  2 -7/8  in.  aperture) 


Fig.  A-5.  Comparison  of  interference  patterns. 
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standard  horn  used  in  the  gain  measurements.  The  angular  discrimination  of  the  main  antenna 
(3-db  beamwidth  »  0.06°)  is  such  as  to  almost  completely  discriminate  against  ground  reflections. 
This  is  borne  out  experimentally  by  the  excellent  symmetry  observed  in  vertical  plane  patterns 
even  when  the  angle  of  declination  is  as  great  as  1°. 

The  distance  separating  the  two  sites  is  equivalent  to  1.15  D  /%  and  consequently  a  slight 
amount  of  axial  defocusing  («0.2\)  is  required  in  the  main  dish  to  simulate  the  far-field  patterns.4 

Early  in  the  program  the  interference  pattern  at  the  receiving  site  was  investigated  with  a 
probe  tower  38-ft  high  and  a  standard-gain  horn  (nominally  25  db)  traversing  the  length  of  this 
tower.  Numerous  runs  were  taken  for  incremental  changes  in  the  elevation  angle  of  the  trans¬ 
mitting  antenna,  and  the  envelope  of  the  interference  pattern  was  adjusted  to  illuminate  properly 
the  28 -ft  main  antenna.  The  final  angular  position  of  the  transmitting  antenna  was  such  that  the 
illumination  of  the  28-ft  main  antenna  was  slightly  greater  toward  the  top  (the  receiving  aperture 
is  located  on  the  lower  side  of  the  transmitted  beam).  The  amplitude  taper  is  approximately 
l/2  db,  which  has  a  negligible  effect  on  the  pattern  measurements.  The  gain  standard  is  positioned 
at  the  midpoint  of  the  vertical  aperture  so  that  the  taper  will  not  result  in  an  error  in  the  gain 
measurement. 

Since  one  of  our  main  objectives  is  to  evaluate  the  effect  of  changing  environmental  condi¬ 
tions  on  the  antenna  characteristics,  the  measurement  errors  caused  by  an  imperfect  range  take 
on  added  importance.  The  effects  of  changing  propagation  conditions  over  the  test  site  with 
changing  environment  must  be  separated  from  actual  changes  in  antenna  performance.  For  pur¬ 
poses  of  analysis,  the  properties  of  the  range  can  be  separated  into  three  categories  which  can 
be  related  to  specific  errors  in  the  antenna  measurements.  These  are: 

(1)  The  interference  pattern  resulting  from  ground  reflections. 

(2)  Refractive  effects  (angle-of-arrival  variations)  caused  by  the  changing 
refractive  index  with  height  relationship. 

(3)  Signal  scintillations  associated  with  atmospheric  turbulence,  etc. 

E.  THE  INTERFERENCE  PATTERN 

The  measurement  of  gain  by  comparing  the  signal  received  from  the  antenna  under  test  with 
that  of  a  calibrated  standard  is  based  on  the  assumption  that  both  antennas  are  placed  in  the  same 
uniform  field.  Because  of  its  high  angular  resolution,  the  main  antenna  receives  only  direct 
radiation  from  the  transmitter  (the  free-space  field).  The  signal  received  by  the  standard  of 
comparison,  however,  because  of  ground  reflections,  depends  on  its  directivity  and  its  vertical 
position. 

A  comparison  of  the  interference  patterns  seen  by  the  standard  horn  and  an  8-in. -diameter 
paraboloid  is  shown  in  Fig.  A-5.  The  half-power  beamwidth  of  the  standard  horn  is  approximately 
8°  as  compared  to  a  3°  half-power  beamwidth  obtained  with  the  8-in.  dish.  It  is  seen  that  the 
interference  pattern  of  the  8-in.  dish  shows  substantially  larger  angular  separations  between 
peaks  and  nulls  and  also  much  larger  intensity  variations. 

It  is  well  known  that  microwave  propagation  over  rough  terrain  is  characterized  by  poor 
specular  reflection  and  results  in  a  signal  intensity  near  the  free-space  level.  The  integrated 
effect  of  a  large  number  of  reflected  signals  with  a  random  phase  relationship  will  alternately 
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reinforce  and  oppose  the  direct  (free-space)  field  as  the  height  is  varied.  It  follows  that  the  free- 
space  level  in  relation  to  the  interference  pattern  seen  by  the  standard  horn  is  the  average  of 
the  loci  of  maxima  and  minima  of  the  interference  pattern.  The  free-space  level  in  relation  to 
the  interference  pattern  seen  by  the  8-in.  dish  is  doubtful,  and,  even  if  it  were  known,  the  error 
in  gain  measurement  would  be  2  db  or  more  depending  on  the  vertical  position  of  the  antenna. 

The  point  of  the  discussion  here  is  that  it  does  not  necessarily  follow  that  the  error  in  the  gain 
measurement  due  to  ground  reflections  is  decreased  by  increasing  the  directivity  of  the  standard. 
In  fact,  unless  the  gain  standard  can  be  made  directive  enough  to  completely  exclude  ground 
reflections,  the  best  approach  is  to  employ  a  very  broad-beamed  standard.  This  would  include 
a  large  area  of  diffuse  reflections  and,  hence,  approach  the  desired  free  space  condition.  The 
25 -db  standard  horn  was  used  in  all  the  gain  measurements:  The  error  in  the  gain  measure¬ 
ments  as  a  result  of  interference  may  be  as  much  as  ±0.5  db,  although  it  is  probably  only  a 
fraction  of  this. 

The  interference  pattern  is  surprisingly  constant.  Probe  measurements  taken  intermittently 
over  a  period  of  three  months  in  a  variety  of  environmental  conditions  showed  but  negligible 
changes  in  the  interference  pattern. 

F.  REFRACTION 

Although  we  were  unable  to  measure  changes  in  the  angle  of  arrival  of  the  incident  signal 
or  at  least  unable  to  separate  this  effect  from  the  beam  scan  and  bore  sight  errors  of  the  system, 
published  accounts  of  angle-of-arrival  measurements  suggest  that  it  may  give  rise  to  a  signifi¬ 
cant  error  in  the  boresight  measurements. 

Sharpless^  describes  angle-of-arrival  measurements  taken  at  a  wavelength  of  1.25  cm.  The 
test  range  is  a  12.6-mile  overland  path  with  an  average  height  of  155  ft  above  ground  level.  The 
results  indicate  that  the  angle-of-arrival  changes  in  the  horizontal  plane  are  slight,  the  maximum 
observed  being  1.8  angular  minutes.  In  the  vertical  plane,  however,  the  variation  is  from  -2.4 
to  +6.6  minutes  relative  to  the  line  of  sight. 

The  present  test  range  is  about  half  as  long  (»6  miles)  but  otherwise  similar  to  the  one  de¬ 
scribed.  Using  the  reported  results  as  a  rough  guide,  we  are  led  to  conclude  that  anomalous 
propagation  conditions  might  cause  a  number  of  minutes’  variation  in  the  RF  angle  of  arrival. 

The  reported  9-minute  variation  represents  the  maximum  observed  over  a  2-month  period.  Also, 
the  range  used  was  on  the  New  Jersey  coast,  much  closer  to  the  ocean  than  the  present  test  site, 
which  may  result  in  quite  different  propagation  conditions. 

In  any  event,  the  beamwidth  of  the  antenna  under  test  is  only  4.2  minutes;  hence,  1  to  2 
minutes  of  boresight  error  is  serious.  The  possibility  of  errors  this  large  or  larger,  as  a  re¬ 
sult  of  refraction,  cannot  be  overlooked. 

Optical  refraction  affecting  the  boresight  optics  used  in  the  evaluation  program  is  the  cause 
of  still  more  uncertainty  in  the  boresight  measurements.  Under  "normal"  conditions, optical  and 
electrical  refraction  on  passage  through  the  atmosphere  are  approximately  the  same  at  these 
frequencies.7  Over  a  short  horizontal  transmission  path,  the  relationship  between  the  variation 
of  optical  and  electrical  refraction  is  unknown. 
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G.  SIGNAL  SCINTILLATION 

Probe  tower  measurements  taken  with  the  8-in.  dish  on  a  typical  daily  run  are  shown  in 
Fig.  A-6.  The  central  section  of  these  traces  (marked  descending)  corresponds  to  the  field 
variation  over  the  height  of  the  tower  which  is  the  superposition  of  the  interference  pattern  and 
the  signal  scintillation.  The  regions  marked  top  and  bottom  show  the  signal  scintillation  with 
the  probe  antenna  held  fixed  at  these  positions.  The  amplitude  of  these  scintillations  varies 
from  a  negligible  amount  (before  dawn)  to  approximately  ±0.6  db  (early  afternoon)  with  rates  on 
the  order  of  1  Hz.  Signal  fluctuations  of  approximately  the  same  rate  and  amplitude  are  observed 
on  both  the  28 -ft  main  antenna  and  the  broader-beamed  standard-gain  horn.  This  indicates  that 
the  scintillations  are  independent  of  ground  reflections. 
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Fig.  A-6.  Probe  tower  measurements. 


This  effect  has  been  reported  elsewhere  and  is  generally  attributed  to  scattering  due  to 
atmospheric  turbulence.  Numerous  runs  of  this  type  were  taken,  and  it  is  apparent  that  the 
scintillations  are  negligible  when  fog  is  present  and  that  they  are  generally  much  more  severe 
in  the  daylight  hours  (especially  early  afternoon)  than  in  the  predawn  hours. 

The  period  of  these  fluctuations  is  comparable  to  the  time  required  to  manually  swing  the 
28-ft  antenna  through  the  beam  when  taking  a  pattern.  As  a  result,  a  slight  lack  of  repeatability 
is  noticed  when  patterns  are  taken  at  times  of  large  signal  scintillation.  Measured  side  lobes 
can  be  in  error  by  as  much  as  1  db  and  beamwidths  by  some  0.2  minute  of  arc. 

These  scintillations  have  a  negligible  effect  on  the  gain  measurement  as  the  signal  is  effec¬ 
tively  "averaged"  out  by  rapid  switching. 

The  errors  in  the  antenna  measurements  which  may  be  attributed  to  the  range  are  sum¬ 
marized  as  follows: 
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1.  Interference  Pattern:—  Because  of  the  extreme  angular  discrimination  of  the  antenna 
under  test,  pattern  and  boresight  measurements  are  unaffected  by  ground  reflections.  The 
interference  pattern  seen  by  the  relatively  broad-beamed  standard-gain  horn  can  cause  errors 
only  in  the  gain  measurement.  The  error  will  not  exceed  ±0.5  db  and  is  probably  less. 

2.  Refraction:-  Changes  in  the  angle  of  arrival  are  of  consequence  only  in  the  boresight 
measurements.  The  amount  of  error  here  is  not  known,  although  during  periods  of  anomalous 
propagation  conditions  it  may  be  large. 

3.  Signal  Scintillations:-  These  intensity  fluctuations  affect  only  the  pattern  characteristics. 
At  times  of  high  scintillation  (usually  early  afternoon)  the  measured  beamwidths  and  side  lobes 
may  be  in  error  by  0.2  minute  and  1  db,  respectively. 

It  may  seem  like  an  oversimplification  to  associate  the  interference  pattern  exclusively  with 
the  gain  measurement,  refraction  exclusively  with  boresight  measurements,  and  scintillations 
exclusively  with  pattern  information,  but  the  extreme  directivity  of  the  antenna  and,  to  a  lesser 
extent,  measurement  techniques  indicate  that  this  breakdown  is  quite  accurate.  In  fact,  this 
analysis  of  measurement  errors  caused  by  the  test  site  appears  to  be  generally  applicable  to 
the  measurement  of  electrically  large  apertures. 

These  discrepancies  in  the  range  are  not  independent  of  each  other;  for  example,  conditions 
which  might  cause  large  refractive  effects  are  certainly  related  to  scintillation.  However,  the 
interdependence  is  of  little  significance  when  identifying  each  of  the  three  propagation  phenomena 
with  a  specific  antenna  characteristic. 

Sharpless6  and  Crawford  and  Jakes9  have  reported  multiple-path  transmission  caused  by 
refraction  rather  than  reflection.  Two,  three,  and  sometimes  more  signal  components  were 
found  to  arrive  at  differing  angles  in  the  vertical  plane.  If  such  an  effect  were  present  on  our 
test  site,  it  would  cause  serious  errors  in  pattern,  gain,  and  boresight  measurements.  For¬ 
tunately,  no  evidence  of  multiple-path  transmission  has  ever  been  observed. 

H.  GAIN 

The  gain  measurement  is  accomplished  by  the  method  of  direct  comparison  with  a  calibrated 
gain  standard.  The  standard  used  is  a  nominal  25-db  Microwave  Associates  pyramidal  horn. 

The  standard  horn  is  mounted  on  one  side  of  the  main  reflector  and  at  the  mid-point  of  the 
vertical  aperture.  The  feed  line  to  the  standard  is  a  long  waveguide  run  (~26ft)  which  is  con¬ 
nected  to  one  port  of  a  waveguide  switch  mounted  in  the  equipment  cab  just  inside  the  vertex  of 
the  antenna.  The  main  feed  line  consists  of  a  waveguide  run  of  approximately  15  ft  and  two 
calibrated  attenuators.  This  line  is  attached  to  the  second  port  of  the  waveguide  switch.  The 
function  of  the  attenuators  is  to  reduce  the  signal  in  the  main  line  to  the  level  of  the  signal  in 
the  standard  horn.  The  common  port  of  the  switch  is  connected  to  the  receiver  and  a  level  meter. 

This  arrangement  enables  rapid  switching  from  the  main  antenna  to  the  standard.  Measure¬ 
ment  errors  as  a  result  of  the  aforementioned  scintillations  are  effectively  averaged  out. 

I.  LINE  LOSSES 

A  most  important  part  of  the  gain  measurement  is  the  accurate  determination  of  the  losses 
in  the  main  and  standard  horn  feed  lines.  It  was  attempted  initially  to  use  a  reflex  klystron  as 
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a  source  for  this  measurement.  However,  a  millimeter  tube  is  subject  to  output  power  variations 
with  physical  handling,  ambient  temperatures,  wind,  etc.  (The  measurement  must  ultimately 
be  performed  outdoors.)  This  approach  was  found  to  be  incapable  of  sufficient  accuracy.  A 
"Tucor"  gas-discharge  noise  course  was  used  as  the  signal  generator.  The  noise  source  proved 
to  be  unaffected  by  the  above  disturbances. 

The  experimental  setup  for  the  loss  measurement  consists  of  the  noise  source  and  an  isolator 
connected  to  the  extreme  end  of  the  line  to  be  measured  (i.e.,  substituted  for  the  standard  horn 
or  the  main  feed,  as  the  case  may  be).  The  rear  end  of  the  line  in  question  is  attached  to  the 
receiver.  The  receiver  output  is  directed  through  a  precise  IF  attenuator  to  a  level  meter. 

With  the  noise  source  on,  the  IF  attenuator  is  adjusted  to  equate  this  level  to  that  which  is 
obtained  with  the  noise  source  off.  The  expression  is: 


T  +  T  =  IT  (1  —  L)  +  T  L  +  T  I  A 
o  m  1  s  o  m 1 


(A-l) 


where 

L  =  Fractional  loss  of  the  line  length  in  question. 

A  =  Fractional  transmission  of  IF  difference  attenuation 
which  is  required  to  equate  levels. 

T0  =  Ambient  temperature  of  line  being  measured. 

Tm  =  Effective  receiver  noise  temperature. 

Tg  =  Effective  source  noise  temperature. 

Solving  for  L, 


This  can  be  shown  to  be 


(A-2) 


(A-3) 


The  quantity  N  /N  ^  determined  by  connecting  the  source  directly  to  the  receiver  input 
terminals.  It  is  then  the  ratio  of  the  noise  power  output  with  the  source  on  and  off  and  can  be 
read  directly  from  the  precision  IF  attenuator. 

The  result  of  these  measurements  shows  the  standard  horn  line  loss  to  be  4.7  db  and  the 
main  line  loss  to  be  2.9  db.  The  main  line  loss  includes  insertion  loss  of  the  two  calibrated  RF 
attenuators  and  the  waveguide  of  the  pattern  receiver.  Repeated  measurements  check  to  within 
0.1  db. 

An  alternate  method  was  used  to  measure  the  relative  attenuation  of  the  two  lines.  Again, 
the  noise  source  was  connected  successively  in  place  of  the  standard  horn  and  the  feed.  The 
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calibrated  attenuator  in  the  main  feed  line  was  adjusted  to  give  equal  levels  of  output  from  the 
receiver.  The  results  of  this  method  checked  with  those  of  the  former  to  within  ±0.1  db. 

The  standard  horn,  as  per  the  manufacturer's  curves,  has  a  gain  of  25.05  db  at  8.5  mm. 

The  calculated  value  is  24.61,  which  should  be  accurate  to  ±0.2  db  (Ref.  10).  Using  the  more  con¬ 
servative  figure  of  24.6  db  for  the  standard  horn,  the  measured  gain  of  the  main  antenna  is 
67.3  db.  This  represents  an  efficiency  of  «53  percent. 

Considering  all  the  possible  causes  of  error,  namely,  the  interference  pattern  seen  by  the 
horn,  line  losses,  calibration  of  the  RF  attenuators,  beam  alignment,  etc.,  the  error  could  be 
1  db  or  more.  It  is  believed,  however,  that  the  above  figure  is  accurate  to  within  ±0.5  db. 

The  gain  figure  of  67.3  db  does  not  include  the  losses  of  the  main  feed  line.  This  figure  is 
the  gain  from  the  feed  and  is  quoted  as  being  an  indication  of  the  quality  of  the  reflecting  surface. 
Operationally,  the  main  line  losses  can  be  reduced  to  1.4  db  (measured);  hence,  the  true  or 
operational  gain  of  the  antenna  is  65.9  db. 

Gain  measurements  were  taken  during  the  daily  runs  to  evaluate  environmental  effects. 

The  previously  described  extremes  of  wind,  temperature,  etc.,  were  encountered.  The  varia¬ 
tion  in  measured  gain  was  within  ±0.25  db  for  90  percent  of  all  the  readings  taken.  There  was 
no  recognizable  trend  in  this  variation  which  could  be  identified  with  time  of  day,  solar  heating, 
wind  conditions,  etc. 

J.  FEED 

11 

The  feed  developed  for  this  antenna  is  a  "rear  feed"  type  consisting  of  an  open-ended 
circular  waveguide  operating  in  its  dominant  mode  (TE^)  with  a  circular  reflecting  disk  posi¬ 
tioned  approximately  3-1/2  wavelengths  in  front  of  the  open  guide.  A  small  hemispherical 
radome  serves  to  hold  the  reflecting  disk  in  position  and  also  provides  weather  protection.  Fig¬ 
ure  A -7  shows  the  feed  in  a  4 -ft  paraboloid.  This  feed  has  proven  highly  satisfactory  in  that  it 


Fig.  A- 7.  Modified  Cassegrainian  feed  in  4— ft  paraboloid. 
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provides  dual  polarization  capability,  wide  bandwidths,  and  nearly  equal  beamwidths  in  the  prin¬ 
cipal  planes.  However,  it  does  have  the  small  disadvantage  of  not  offering  a  point  center  of 
phase  and  consequently  there  is  a  slight  difference  between  the  effective  phase  centers  in  the 
two  principal  planes. 

A  similar  type  of  feed  with  the  splash  plate  in  close  proximity  to  the  open-ended  guide  is 
9  12 

described  in  the  literature.  *  This  feed  behaves  more  like  a  ring  source  than  a  point  source, 
and  consequently  it  is  considered  not  suitable  for  use  with  a  reflector  having  a  point  focus. 

The  present  feed  differs  from  this  "ring  focus"  feed  in  that  the  reflecting  disk  is  well  into 
the  far  field  of  the  open-ended  waveguide.  The  result  is  a  configuration  which  is  Cassegrainian 
in  behavior,  with  the  hyperboloid  sub  reflector  having  degenerated  into  a  flat  plane.  The  image 
of  the  open-ended  guide  is  placed  at  the  focal  point  of  the  paraboloid. 


Fig.  A-8.  Patterns  of  modified  Cassegrainian  feed  in  4— ft  paraboloid, 

H-plane  focused. 

The  phase  error  with  this  arrangement  is  reduced  to  within  tolerable  limits.  Figure  A-8 
indicates  the  principal  plane  patterns  taken  in  a  4-ft  paraboloid  with  the  axial  setting  of  the 
feed  such  as  to  optimize  the  H-plane  patterns.  The  difference  between  the  effective  phase  cen¬ 
ters  in  the  two  principal  planes  is  approximately  0.3  wavelength.  The  compromise  feed  setting 
results  in  similar  principal  plane  patterns  with  the  side  lobes  clearly  distinguishable  from  the 
main  beam  (i.e.,  the  side  lobes  are  not  integrated  into  the  main  beam  as  in  the  E-plane  pattern 
of  Fig.  A -8).  This  slight  defocusing  has  no  appreciable  effect  on  gain  or  beamwidths. 

A  more  detailed  description  of  the  feed  can  be  found  in  Ref.  11. 

K.  PATTERNS 

Typical  azimuth  and  elevation  plane  patterns  of  the  28-ft  antenna  with  the  feed  positioned  to 
focus  the  elevation  plane  patterns  are  shown  in  Fig.A-9.  There  is  a  close  similarity  between 
these  patterns  and  those  taken  in  the  proven  4-ft  paraboloid.  The  amount  of  defocusing  in  the 
E-plane  patterns  of  both  sets  is  comparable.  With  the  exception  of  the  side-lobe  structure  in 
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Fig.  A-9.  Typical  azimuth  and  elevation  patterns,  elevation 
plane  focused. 


Fig.  A-10.  Azimuth  and  elevation  plane  patterns  for  28— ft 
paraboloid  with  feed  in  compromise  position. 
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the  E -plane,  the  beamwidths  and  side  lobes  for  both  antennas  check  exceptionally  well  when 
expressed  in  terms  of  A/D.  This  is  a  good  indication  of  the  excellent  performance  of  the  28-ft 
antenna  as  well  as  the  quality  of  the  test  site.  In  both  cases,  the  difference  in  axial  setting 
corresponding  to  the  different  focal  points  in  the  principal  planes  is  approximately  0.3  X. 

It  is  apparent  that  the  inability  to  obtain  perfectly  focused  patterns  in  both  planes  for  a  single 
axial  setting  of  the  feed  is  a  characteristic  of  the  feed  itself  and  not  the  result  of  a  distorted 
reflecting  surface. 

Figure  A-10  shows  azimuth  and  elevation  plane  patterns  with  the  feed  in  the  compromise 
position.  The  main  beams  of  these  patterns  are  almost  identical;  the  only  difference  in  the  two 
is  in  the  side-lobe  structure.  Patterns  showing  the  far-out  side-lobe  structure  are  shown  in 
Fig.  A-l  1 . 


Fig.  A— 1 1 .  Azimuth  and  elevation  plane  patterns  for  28-ft  paraboloid 
showing  far-out  side  lobes. 

In  the  pattern  measurements  described  below,  which  are  primarily  to  determine  the  effects 
of  solar  heating,  the  axial  position  of  the  feed  was  set  at  the  elevation  plane  focal  point.  This 
was  done  because  the  effects  of  solar  heating  are  more  pronounced  in  the  elevation  plane  (the 
vertical  temperature  gradients  are  larger  than  the  horizontal)  and,  hence,  will  be  more  easily 
distinguishable  if  the  elevation  plane  is  focused. 

Thermocouples  were  attached  to  the  aluminum  backup  structure  to  measure  temperature 
gradients.  The  results  indicate  a  maximum  vertical  temperature  gradient  (in  the  early  afternoon) 
of  approximately  30°  F.  The  horizontal  gradient  is  about  half  this  figure.  The  antenna  is  directed 
almost  due  north  to  receive  a  signal  from  the  transmitter  and,  consequently,  the  late  summer 
sun  swings  across  the  back  of  the  reflector. 

Figure  A -12  shows  the  effect  of  solar  heating  on  the  elevation  plane  pattern.  The  high  asym¬ 
metrical  side  lobe  varies  from  a  level  of  17-1/2  db  to  a  condition  of  symmetry  with  both  side 
lobes  at  24  db.  The  17-l/2-db  side  lobe  is  the  maximum  experienced  in  all  the  patterns  taken 
to  evaluate  the  solar  heating  effect.  On  the  particular  day  in  question,  the  sun  became  obscured 
in  the  early  afternoon  and,  as  seen  in  the  figure,  the  pattern  reverted  back  to  the  predawn  con¬ 
dition.  A  close  correlation  between  asymmetry  and  solar  heating  was  consistently  observed. 
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Fig.  A-12.  Elevation  pattern  (H-plane) 
for  the  28— ft  paraboloid  showing  effect 
of  solar  heating. 


Fig.  A-13.  Azimuth  pattern  (E-plane) 
showing  effect  of  solar  heating. 


MINUTES  OF  ARC 
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The  guy  wires  which  support  the  feed  column  are  attached  to  the  ends  of  the  aluminum  girder 
panels.  One  result  of  the  higher  temperatures  on  the  top  of  the  backup  structure  is  believed  to 
be  a  relaxing  of  tension  in  the  top  guy  wire.  This  causes  a  downward  motion  of  the  feed  column 
and  an  upward  scan  of  the  beam.  The  high  asymmetrical  side  lobe  appears  on  the  proper  side 
of  the  scanned  beam  to  be  identified  as  the  "coma"  lobe  associated  with  beam  scan.  As  the  beam 
is  scanned  off-axis,  due  to  a  lateral  motion  of  the  feed,  the  "coma"  lobe  appears  on  the  axial 
side  of  the  main  beam.  The  side-lobe  variation  in  Fig.  A -12  shows  the  asymmetrical  side  lobe 
high  in  the  predawn  pattern  and  decreasing  with  solar  heating.  This  "reversed  scan"  is  the 
result  of  positioning  the  feed  for  best  elevation  plane  focus  at  a  time  of  maximum  thermal  gra¬ 
dients  (afternoon).  The  up-and-down  designations  in  the  figure  refer  to  the  position  of  the  trans¬ 
mitter  relative  to  the  peak  of  the  elevation  plane  pattern. 

A  relaxing  of  tension  in  the  top  guy  wire  must  certainly  be  accompanied  by  some  thermal 
distortion  of  the  upper  portion  of  the  reflecting  surface  which  would  also  result  in  asymmetry. 

The  obvious  approach  to  determine  how  much  pattern  asymmetry  is  the  result  of  motion  of  the 
feed  column  and  how  much  is  the  result  of  reflector  distortion  would  be  to  measure  the  beam 
shift.  The  pattern  asymmetry  which  would  result  from  reflector  distortions  would  be  accom¬ 
panied  by  a  negligible  amount  of  beam  scan.  Unfortunately,  the  accuracy  of  the  boresight  shift 
measurements  (described  below)  is  not  sufficient  for  the  amount  of  beam  scan  involved.  Pattern 
measurements,  however,  do  indicate  that  the  asymmetry  due  to  distortion  of  the  reflecting  sur¬ 
face  is  negligible  in  comparison  to  the  "coma"  lobe. 

Consider  the  case  (Fig.  A-12)  where  the  predawn  position  of  the  beam  is  presumed  to  be 
downward  and  is  scanned  upward  by  the  effects  of  solar  heating  to  a  condition  of  symmetry.  In 
view  of  the  excellent  patterns  obtained  with  the  maximum  temperature  gradient  (symmetrical; 
well-focused,  etc.)  associating  the  change  in  side-lobe  structure  with  thermal  distortion  of  the 
reflecting  surface  leads  to  the  highly  unlikely  situation  of  the  surface  becoming  more  paraboloidal 
as  the  temperature  gradient  increases.  Further,  patterns  were  taken  with  the  feed  position 
adjusted  such  that  the  predawn  position  of  the  beam  is  on  axis  and  is  scanned  upward  with  solar 
heating  to  an  asymmetrical  condition.  Comparison  of  the  predawn  asymmetrical  pattern  (neg¬ 
ligible  temperature  gradient)  of  the  former  with  the  early  afternoon  asymmetrical  pattern  (maxi¬ 
mum  temperature  gradient)  of  the  latter  shows  them  to  be  mirror  images.  This  is  the  case  with 
patterns  representing  equal  amounts  of  scan  in  opposite  directions  and  can  be  true  only  if  the 
asymmetry  caused  by  thermal  distortion  of  the  reflector  is  negligible. 

During  the  initial  focusing  period,  it  was  found  that  a  lateral  displacement  of  the  feed  amount¬ 
ing  to  approximately  1  / 1 6  in.  changed  the  pattern  from  a  condition  of  symmetry  to  a  typical  off- 
axis  pattern  with  the  "coma"  lobe  at  the  17-  to  18-db  level.  The  l/16-in.  motion  of  the  feed, 
which  has  about  the  same  effect  as  the  maximum  temperature  gradient,  corresponds  to  a  beam 
scan  of  approximately  1.4  angular  minutes. 

Figure  A -13  shows  the  effect  of  solar  heating  on  the  azimuth  patterns.  It  is  seen  that  the 
asymmetry  is  slight.  Presumably  this  is  due  to: 

(1)  The  thermal  gradient  from  left  to  right  (15°  —  20 °F)  is  on  the  order 
of  half  the  top-to-bottom  gradient. 

(2)  The  slight  defocusing  in  the  azimuth  plane  serves  to  mask  asymmetrical 
effects  caused  by  temperature  gradients. 
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The  only  consequence,  then,  of  nonuniform  solar  heating  on  the  radiation  patterns  appears 
to  be  the  beam  scan  in  the  vertical  plane  amounting  to  approximately  1.4  angular  minutes.  This 
amount  of  scan  corresponds  to  the  maximum  asymmetry  observed  during  the  test  period  although 
on  clear  days  the  ’'coma”  lobe  consistently  reaches  the  17-  to  18-db  level.  Thermal  distortion 
of  the  reflecting  surface  seems  To  have  a  negligible  effect  on  pattern  asymmetry  relative  to  the 
effects  of  beam  scan. 

L.  AMBIENT  TEMPERATURES 

Measurements  were  taken  in  ambient  temperatures  which  ranged  from  the  mid  30' s  (°K)  to 
the  high  80' s.  No  noticeable  changes  of  the  electrical  characteristics  were  observed. 

It  might  be  mentioned  here  that  the  antenna  was  left  unprotected  on  the  roof  of  one  of  the 
Lincoln  Laboratory  buildings  for  the  winter  prior  to  the  present  series.  Pattern  and  gain 
measurements  taken  before  the  winter  are  in  close  agreement  with  the  present  results. 

M.  WIND 

Measurements  were  taken  during  periods  when  the  wind  velocity  reached  25  mph  with  gusts 
of  approximately  30  mph.  The  wind  direction  varied  from  the  front  to  the  back  of  the  dish.  No 
first-order  effects  were  observed  in  gain  or  pattern  measurements.  The  only  noticeable  effect 
of  the  wind  was  a  slight  oscillation  in  the  level  of  the  received  signal,  indicating  some  vibration 
in  the  antenna  or  more  probably  in  the  feed  column.  This  amounted  to  not  more  than  0.2  to  0.3db 
superimposed  on  the  "wind-free”  pattern. 

N.  BORESIGHT  MEASUREMENTS 

The  variation  in  the  alignment  of  the  boresight  optics  and  the  position  of  the  RF  beam  as  a 
function  of  solar  heating  and  wind  conditions  proved  to  be  most  difficult  to  evaluate.  The  most 
significant  contribution  to  this  difficulty  is  believed  to  be  refractive  effects,  with  both  electrical 
and  optical  angle-of-arrival  being  effected  relative  to  each  other  in  an  unknown  manner. 

Seven  consecutive  daily  runs  were  taken,  during  which  time  no  adjustments  were  made  on 
either  the  feed  position  or  the  synchros.  The  boresight  measurements  taken  over  this  period, 
although  not  as  accurate  as  might  be  desired,  give  some  indication  of  the  magnitude  of  both  the 
electrical  and  optical  refraction  encountered. 

A  total  of  about  50  measurements  were  made  during  these  seven  daily  runs.  Each  measure¬ 
ment  consisted  of  optical  line-of-sight  and  RF  beam  position  on  both  azimuth  and  elevation 
synchro  readouts.  The  range  of  variation  of  these  quantities  experienced  during  this  measure¬ 
ment  period  is  as  follows: 


Azimuth: 

Optical 

Electrical 

Elevation: 

Optical 

Electrical 


1/2  to  1  minute 
1  to  1  - 1  / 2  minute 

3  to  3-1/2  minutes 
3  to  3-1/2  minutes 
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There  is  no  clear  relationship  between  the  relative  optical  and  electrical  variations  in  either 
plane.  Initially,  the  optical  and  electrical  boresights  were  aligned  parallel  and  at  times  during 
the  measurement  period  they  are  separated  by  » 1-1/2  minutes  in  the  azimuth  plane  and  more 
than  2  minutes  in  the  elevation  plane.  Possible  errors  in  these  measurements  are:  the  synchro 
readout  (as  much  as  1/2  minute),  human  errors  in  optical  alignment,  and  errors  in  aligning  the 
RF  beam.  Hence,  the  azimuth-optical  variation  of  l/2  to  1  minute  is  about  as  small  an  angle  as 
the  system  is  capable  of  measuring. 

Clearly,  there  is  more  variation  of  both  the  optical  and  electrical  boresights  in  the  elevation 
plane.  There  is  no  reason  to  suspect  that  any  of  the  above  sources  of  error  would  be  more 
pronounced  in  the  elevation  plane  over  the  azimuth  plane.  The  variation  of  the  elevation  plane 
quantities  may  simply  be  positional  shifts  of  the  entire  mount  because  of  flexible  mount  supports 
or  lack  of  rigidity  of  the  base.  A  second  possibility  is  the  presence  of  both  optical  and  RF  re¬ 
fraction.  References  6  and  12  indicate  that  angle-of-arrival  variations  of  a  number  of  minutes 
might  be  expected.  If  we  assume  that  the  basic  accuracy  of  the  mount  is  adequate,  the  amount 
of  optical  refraction  for  the  period  in  question  would  appear  to  be  in  the  vicinity  of  2  to  3  minutes. 
The  RF  refraction  is  at  least  the  same,  probably  more.  It  is  thought  that  an  appreciable  part 
of  the  variation  in  the  elevation  plane  can  be  attributed  to  refraction,  if  it  is  not  the  prominent 
cause. 

Another  important  characteristic  to  be  evaluated  is  the  effect  of  elevation  angle  on  boresight 
and  patterns.  Changes  in  gravitational  loading  could  give  rise  to  beam  scan  and  pattern  distortion. 
Although  actual  measurements  have  not  been  made,  there  are  indications  that  boresight  and 
patterns  are  preserved  with  elevation  angle.  The  reflector  is  fabricated  while  pointing  at  zenith, 
and  hence,  the  horizontal  position  probably  represents  the  "worst  case"  of  gravitational  loading. 

In  that  the  recorded  patterns  were  taken  in  the  horizontal  position  and  show  no  evidence  of 
pattern  distortion,  it  can  be  assumed  that  pattern  distortion  at  the  higher  elevation  angles  is 
slight.  The  beam  can  be  held  to  the  optical  reference  to  at  least  a  beamwidth,  probably  better. 

0.  CONCLUSIONS 

It  is  clear  that  the  antenna  performs  quite  satisfactorily  at  a  wavelength  of  8.5  mm.  The 
exceptional  aperture  efficiency  («53  percent)  and  patterns  obtained  indicate  that  the  reflecting 
surface  is  paraboloidal  to  a  high  degree  of  accuracy  and  retains  its  shape  in  the  environmental 
conditions  encountered. 

The  only  measured  change  in  electrical  characteristics  resulting  from  the  wind,  nonuniform 
solar  heating  and  ambient  temperatures  experienced  is  a  lateral  motion  of  the  feed  column  with 
nonuniform  heating,  causing  the  beam  to  scan.  The  amount  of  scan  is  less  than  1.5  angular 
minutes.  There  is  no  evidence  of  any  significant  distortions  of  the  reflecting  surface  as  a  re¬ 
sult  of  these  environmental  conditions. 

It  was  shown  that  with  proper  precautions,  the  test  range  does  not  introduce  appreciable 
errors  into  the  gain  and  pattern  measurements.  The  influence  of  refraction,  however,  is  be¬ 
lieved  to  appreciably  affect  the  boresight  measurement s,  although  the  amount  which  might  be  ex¬ 
pected  during  various  weather  conditions  is  not  clear. 

While  there  is  still  some  uncertainty  in  regard  to  performance  at  high  angles,  in  addition  to 
the  uncertainty  in  the  boresight  characteristics,  recent  radiometric  experiments  involving  optical 
tracking  of  celestial  sources  indicate  that  the  optical  control  of  the  beam  is  quite  satisfactory. 
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EXTENSION  OF  BOUNDARY  PERTURBATION  THEORY 


The  classical  application  of  boundary  perturbation  theory  to  the  solution  of  electromagnetic 
scattering  problems  is  limited  to  boundaries  which  are  close  to  those  few  boundary  shapes  for 
which  exact  solutions  are  known.  By  continuing  the  general  approach  to  boundary  perturbation 
theory  introduced  by  Mitzner,*  it  has  been  possible  to  write  the  first-order  perturbation  con¬ 
tribution  to  the  scattered  field  in  a  particularly  simple  form  (Eq.(B-l)]  which  is  easily  applied 
to  any  unperturbed  boundary  shape  for  which  the  surface  field  quantities  can  be  computed.  Thus 
it  becomes  possible  to  deal  with  problems  which  are  amenable  neither  to  physical  optics  methods 
alone  nor  to  perturbation  theory  alone,  by  first  using  physical  optics  to  obtain  the  field  quantities 
at  the  surface  of  the  large  "bumps,"  and  then  regarding  the  small  "bumps"  as  perturbations  of 
the  larger  "bumps,"  completing  the  solution  using  perturbation  theory. 

The  first-order  contribution  6* E  to  the  field  scattered  by  an  arbitrary  boundary  between 
two  arbitrary  isotropic  linear  media  has  been  derived  as 

fi1!:  •  p'  =  r  {AE  •  D'-  AD  .  E'-  AH  •  B'  +  AB  •  H'}  hdS  .  (B-l) 

Here  AA  is  the  discontinuity  in  the  unperturbed  field  A  across  the  unperturbed  boundary,  A'  is 
the  field  that  would  exist  at  the  unperturbed  boundary  were  an  electric  dipole  source  p'  to  be 
placed  at  the  point  where  6  E  is  to  be  evaluated,  h  is  the  height  function  of  the  perturbations 
and  the  integration  is  carried  out  over  the  unperturbed  surface.  The  sign  convention  is  that 
AA  =  A^  —  A^  where  A^  [A^]  is  the  field  in  the  region  toward  which  perturbations  are  reckoned 
as  positive  [negative].  No  ambiguity  arises  in  the  value  to  be  taken  by  the  primed  fields,  be¬ 
cause,  for  example,  AE  is  a  vector  normal  to  the  surface  and  the  normal  component  of  D'  is 
continuous  at  the  surface.  Thus  only  the  continuous  components  of  the  primed  fields  are  involved. 

The  precise  conditions  to  be  imposed  on  the  form  of  the  unperturbed  surface  and  on  h  are 
unknown.  One  supposes  that  the  unperturbed  surface  can  be  perturbed  only  in  regions  free  of 
field  singularities  and  that  both  h/\  and  Vh  should  be  small  compared  to  unity,  where  \  is  the 
smaller  electromagnetic  wavelength  of  the  two  regions. 

The  Lincoln  Calibration  Sphere  (LCS)  program^  provided  an  opportunity  to  apply  Eq.(B-l)  to 
compute  the  far-field  backscattered  by  a  perfectly  reflecting  body.  In  this  case,  since  tangential 
E'  and  normal  B'  are  zero  at  the  surface,  Eq.(B-l)  reduces  to 


6^E 


E’  = 
o 


k2  eikz 

47 re  z 
o 


I 


(eoE 


+  (i  H 
ro 


H')  hdS 


( B  —  2) 


where  the  superfluous  jump  operators  have  been  omitted  and  the  dipole  p'  has  replaced  by  the 
equivalent  plane  wave  E^  incident  on  the  sphere.  Thus  E,  H  and  E',  H'  are  the  total  fields  at 


*  k.  M.  Mitzner,  “Effect  of  Small  Irregularities  on  Electromagnetic  Scattering  from  an  Interface  of  Arbitrary  Shape,”  J.  Math. 
Phys.  5,  1776-1786  (1964). 

t  R.  T.  Prosser,  ‘‘The  Lincoln  Calibration  Sphere,”  IEEE  Proc.  (Correspondence)  53,  1672  (1965). 
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the  surface  of  the  perfect  sphere  when  it  is  illuminated  by  the  plane  waves  exp{— ikz}  and 
E'  exp  {—ikz},  respectively. 

°  2 

The  LCS  is  a  hollow  aluminum  sphere  nominally  1  m  in  cross-sectional  area.  Since  it  was 

machined  on  a  vertical  boring  mill  after  assembly,  its  shape  is  essentially  that  of  a  body  of 
revolution.  Thus  the  height  function  h  of  the  surface  perturbations  is  independent  of  the  azi¬ 
muthal  angle,  and  can  be  expanded  as 

oo 

h(G)  =  a  l  hnPn(cos0)  ,  (B-3) 

n=2 

where  a  is  the  radius  of  the  mean  sphere  and  P  (cose)  is  the  Legendre  polynomial  of  degree  n. 

No  generality  is  lost  by  assuming  that  the  natural  axis  of  the  "sphere"  lies  in  the  plane  y  =  0 
and  is  inclined  at  an  angle  0  to  the  z-axis.  Then,  by  choosing  Eq  and  E^  each  to  be  either  in 
the  x  or  y  direction,  the  following  dyadic  form  for  the  total  backscattered  field  E  (to  first 
order)  can  be  derived 

E  «  6 'E  +  6*E  =  (xaxxx  +  yayyy)  •  6°E  ,  (B-4) 

where 


yy 


=  1  -i2ka  Z  hn{bnPn(COS0)  *  Pn  (^rl{l  Pn(cos0,) 


n=2 


(B-5) 


and  6°E  is  the  zeroth-order  backscattered  field.  The  coefficients  b^  and  p^  are  given  by  certain 
integrals  of  combinations  of  field  components  over  the  surface  of  the  sphere.  The  field  com¬ 
ponents  can  be  computed  exactly  using  their  Mie  series  solutions  or  can  be  approximated  over 
the  lit  region  of  the  sphere  by  using  physical  optics.  The  former  method  is  then  a  direct  appli¬ 
cation  of  boundary  perturbation  theory;  the  latter  is  an  application  of  the  dual  approach  in  which 
the  small  "bumps"  are  regarded  as  perturbations  of  the  larger  "bump"  which  is  the  sphere 
itself. 

A  comparison  of  the  two  methods,  for  ka  =  30,  is  illustrated  in  Fig.  B-l,  which  is  a  plot  of 

2  2 

|b^|  and  (n  —  2)L  |  p^  |  /(n  +  2)'.  It  can  be  seen  that  there  is  very  good  agreement  for  all  n  <  1.5ka, 
and  the  same  is  true  when  a  phase  comparison  is  made.  That  the  agreement  is  destroyed  at 
n  ~  2ka  is  to  be  expected,  because  the  phase  variation  of  the  field  products  in  the  integrand  in 
Eq.(B-2)  matches  that  of  the  spherical  harmonic  of  degree  2ka  only  in  the  vicinity  of  the  shadow 
boundary.  Thus  the  maximum  contribution  to  b^  and  p^,  for  n  close  to  2ka,  comes  from  a 
region  of  the  sphere  where  physical  optics  is  no  longer  valid,  and  a  more  detailed  analysis 
based  on  the  methods  of  Fock  should  be  used. 

The  expressions  for  the  b^  and  p^,  valid  for  large  ka  and  n/2ka  <1—6,  computed  using 
physical  optics  to  specify  the  zeroth-order  fields  and  used  for  the  dashed  curves  in  Fig.  B-l  are 


fc>n  «  i‘n+1  h^2)(2ka)/h|2)(2ka)  ,  n=  0,1,2... 


P 


n 


(n  —  l)n(n+l)(n  +  2) 
(2n  -  1)  (2n  +  1)  (2n  +  3) 


{( 2n  +  3)  bn_2  -  2(2n  +  1)  b^  +  (2n  -  1)  t>n+2) 
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n  /Zka 

Fig.  B-l.  The  b  and  p  coefficients  as  a  function  of  n/2ka. 

n  n  ' 

(2) 

where  hn  (y)  is  the  spherical  Hankel  function  of  the  second  kind  of  order  n.  Simpler  expres¬ 
sions,  which  also  lose  accuracy  for  very  low  values  of  n,  are 

_  i  /? 

t>n«(cos9o)  '  exp  [i2ka(  1  —  cos  0q)  —  i(n  +  1/2)  9^) 

2  2 

p  ~  —  n  sin  0  b 
rn  on 

where  9q  =  arcsin[(n  +  l/2)/2ka]. 

The  shape  of  the  and  p^  curves  shown  in  Fig.  B-l  is  roughly  the  same  for  all  large  ka 

except  that  the  height  of  the  maximum  increases  with  increasing  ka.  A  rough  analysis  indicates 

l/3 

that  the  height  of  the  peak  is  proportional  to  (ka)  '  . 

Careful  profile  measurements  of  one  of  the  spheres  produced  during  the  LCS  program  were 

made  and  from  them  the  hn  coefficients  for  the  expansion  in  Eq.(B-3)  were  computed.  Then, 

using  Eqs.  ( B-4)  and  (B-5)  the  variation  in  the  backscattering  cross  section  was  computed  for 

the  frequency  of  9.4  GHz  (ka  =  111),  at  which  a  measurement  of  the  sphere  on  a  radar  range  was 

made.  The  results  are  shown  in  Fig.  B-2,  where  the  range  measurement  is  compared  with  the 

cross  section  computed  from  both  the  exact  and  the  approximate  forms  for  the  b  and  p  .  The 

n  n 

0-db  level  for  the  range  measurement  is  set  arbitrarily,  because  the  sphere  was  better  than 
any  calibration  target  available.  (The  RMS  surface  deviation  of  this  sphere  was  3.2  x  lo"3  in.) 
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Fig.  B-2.  Comparison  of  measured  cross-section  variation  with  that  computed 
using  both  approximate  and  exact  bn  and  pn  coefficients. 

Fair  agreement  appears  to  exist  between  the  measured  and  computed  cross-section  varia¬ 
tion.  The  background  level  of  the  range  varied  between  -30  to  -40dbm  and  so  could  account 
for  discrepancies  as  large  as  ±0.15  db.  The  two  computed  curves  are  practically  indistinguishable 
from  one  another,  owing  to  the  fact  that  all  the  surface  variation  of  the  sphere  was  attributable 
to  spherical  harmonics  of  relatively  low  degree,  for  which  the  approximate  bn  and  pn  coefficients 
have  good  accuracy. 

A  further  theoretical  result  of  interest  connects  the  variance  of  a  particular  vector  com¬ 
ponent  61  E  of  6i:E  to  the  power  spectral  components  wn  of  statistically  isotropic  surface  per¬ 
turbations  on  a  perfectly  conducting  sphere.  Thus  if 

oo 

<h(e,  (P)  h<e\  <p')>  =  a2  Yj  wnPn(cosy) 

n=2 

where  (  y  denotes  the  expected  value  and  cosy  =  cos0  cosG'  +  sin  0  sin0'  cos  (<p  —  </?'),  then 

oo 

<I«4E|2>  =  |6°E|24k2a2  £  wngn  , 

n-  2 

where  gn  is  |bj2  +  (n  -  2)'.  |pn|2/[2(n  +  2)'.]  and  (n  -  2)'.  |  pn  |  2/[2(n  +  2)'.]  for  the  like  polarized 
and  orthogonally  polarized  components,  respectively,  when  the  incident  plane  wave  is  linearly 
polarized.  The  corresponding  values  for  gn  when  the  incident  wave  is  circularly  polarized  are, 
respectively,  (n  -  2)1  |pn|2/(n-  2)1  and  |t>n|2.  This  explains  the  choice  of  the  quantities  plotted 
in  Fig.  B-l . 
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2 

A  related  formula  for  the  variance  <(6ct)  y  of  the  backscattering  cross  section  is 

<(da)2>  =  a2l6k2a2  V  w  f  , 
n 

where  a  is  the  nominal  (unperturbed)  backscattering  cross  section  and  f  is  |lm{b  }  I  + 
o  2  2  n  n 

(n  -  2)1  I Im  {pn}  |  /[2(n  +  2)1]  or  |lm{bn}  |  depending  on  whether  the  measurement  is  made  with 

linear  or  circular  polarization,  respectively. 

M.  L.  Burrows 
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